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In Re Application of: 
AKERS, Douglas, W. 
Serial No. 10/788,743 
Filing Date: February 25, 2004 



For: 



METHOD FOR ON-LINE 
EVALUATION OF MATERIALS 
USING PROMPT GAMMA RAY 
ANALYSIS 



Atty Dkt: B-200 



Examiner: Palabrica, R.J. 
Group Art Unit: 3663 
Confirmation No.: 6111 



APPEAL BRIEF 

To: Commissioner of Patents 
P.O. Box 1450 
Alexandria, VA 22313-1450 

Sir: 

This Appeal Brief is submitted in response to the non-final rejections of the claims 
the office action, mailed January 16, 2007. A notice of appeal was filed on January 26, 
2007. 



REAL PARTY-IN-INTEREST 

The assignee of the entire right, title, and interest in the patent application is Bechtel 
BWXT Idaho, LLC. 
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RELATED APPEALS AND INTERFERENCES 

There are three related appeals of other United States patent applications that may 
directly affect, or be directly affected by, or have a bearing on, the Board's decision. 
Divisional application, serial no. 10/269,807, filed October 10, 2003, came before the Board 
in Appeal No. 2005-0855. The Board mailed its decision to the appellant on May 31, 2005. 
United States Patent No. 7,058,153 subsequently issued on June 6, 2006. Divisional 
application, serial no. 10/383,096, filed March 5, 2003, is currently pending on appeal in the 
briefing stage. In parent application, serial no. 09/932,531, filed August 17, 2001, an appeal 
brief was filed on June 21, 2006. Thereafter, on February 23, 2007, a notice of allowance was 
issued, allowing the claims at issue in the appeal. The current application is a continuation-in- 
part of parent application, serial no. 09/932,531, and a continuation-in-part of divisional 
application, serial no. 10/383,096. 

There are currently no related interferences known to appellant, appellant's legal 
representative, or the assignee which will directly affect, or be directly affected by, or have a 
bearing on, the Board's decision. 
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STATUS OF THE CLAIMS 

Claims 1-10, 12-19 and 21-23 are pending in the application. Claims 1 1, 20 and 24 are 
cancelled. Claims 1-10, 12-19 and 21-23 currently stand rejected. The rejections of claims 1- 
10, 12-19 and 21-23 are appealed. 

In the present application, the examiner provisionally rejected claim 2 under the 
judicially-created doctrine of double patenting over claim 3 of co-pending divisional 
application, serial no. 10/383,096. The appellant did not traverse this provisional rejection 
and has agreed to file the appropriate terminal disclaimer in the appropriate application at the 
appropriate time (i.e., upon the indication of allowance of either claim 2 in this application or 
claim 3 of the co-pending application). 
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STATUS OF AMENDMENTS 

No amendments were filed or entered subsequent to the non-final office action mailed 
January 16, 2007. 



SUMMARY OF CLAIMED SUBJECT MATTER 

The present invention is directed to a method for evaluating a material specimen by 
bombarding it with neutrons to create prompt gamma rays. The invention as claimed is 
summarized below with reference to the independent claims and claims separately argued. 
Claims 1 and 12 are independent claims. Claims 2, 13-19, 21 and 23 are dependent claims 
argued separately. The claims contain reference numerals and reference to the specification 
and drawings. All references are shown in the application at least where indicated herein. 
1. A method (226; Fig. 9; |72) for evaluating a material specimen (12, 1 12, 212; Figs. 1, 
7-9; U123-31, 33-34, 37-40, 43, 47, 50 : 53, 055-57, 59-63, 66, 68-69, 71-75, 77), comprising: 

mounting (228, Fig. 9; H1J37, 59-61, 72-74, 77) a neutron source (14, 54, 1 14, 214, 
254; Figs. 1, 7-9; ffl[23-30, 32-34, 36-37, 50, 52, 56-63, 72-74, 77) adjacent the material 
specimen (12, 112,212); 

mounting (Fig. 9; 38, 40, 66-69) a detector (16, 30, 32, 1 16, 130, 216, 230; Figs. 1, 
4, 6-9; HTI23-24, 36, 38-42, 44-46, 49, 52-54, 56, 60, 65-72, 74-77) adjacent the material 
specimen (12, 112,212); 

bombarding (232; Fig. 9; 1H[23-24, 30-31, 34, 37, 50, 56, 59-61,73) the material 
specimen (12, 112, 212) with neutrons (Figs. 1, 7-9; ffif30-34, 37, 50, 52, 56-63, 73) from the 
neutron source (14, 1 14, 214) to create prompt gamma rays (Figs. 1, 7-9; ffl|23-24, 27, 30, 34, 
37-39, 41, 43-45, 49-53, 57, 60-61, 63, 66-67) within the material specimen (12, 1 12, 212), 
some of the prompt gamma rays being emitted from the material specimen (12, 1 12, 212), 
some of the prompt gamma rays resulting in the formation of positrons (Figs. 1, 7-9; 1ffl23-25, 
27, 30-31, 34, 37, 43, 50-52, 56-57, 60-61, 63) within the material specimen (12, 1 12, 212) by 
pair production (18, 118, 218; Figs. 1, 7-9; ffl[23, 28, 30-31, 37, 50, 52, 57, 60-61, 63); 
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collecting (70, 82, 234; Figs. 3, 6, 9; 175) positron annihilation data (22, 122, 222; 
Figs. 1, 3, 4, 6-9; H 38, 42-45, 49, 53-54, 67, 70-71, 75-77) by detecting with the detector 
(16, 30, 32, 1 16, 130, 216, 230) a plurality of emitted annihilation gamma rays (Figs. 1, 7-9); 

1123- 24, 30, 38-39, 41, 43, 45, 51-53, 57, 66-67) resulting from the annihilation (H24, 30, 
43-44, 57) of positrons, the detector (16, 30, 32, 116, 130, 216, 230) producing the positron 
annihilation data (22, 122, 222); 

processing (84, 236; Figs. 6, 9; H26, 42, 47, 49, 53-54, 67, 70, 75, 77) collected 
positron annihilation data (22, 122, 222) in accordance with a Doppler-broadening algorithm 
(40; Figs. 2, 6; 1126-27, 38, 42, 47-49, 51, 70, 75); and 

continuing (238; Fig. 9; 175-76) to collect and process positron annihilation data (22, 
122, 222) to measure an accumulation of lattice damage (1124-26, 43, 47, 60-61, 67-69, 71, 
74, 77) over time. 

2. The method of claim 1, further comprising: 

collecting (234; Fig. 9, 175-76) prompt gamma ray data (20, 120, 220; Figs. 1, 4, 7-9, 

1124- 25, 38, 42-45, 49, 53-54, 67, 70-71, 75-77) by detecting with the detector (16, 30, 32, 

1 16, 130, 216, 230) a plurality of emitted prompt gamma rays (Figs. 1, 7-9; 1123-24, 27, 30, 
34, 37-39, 41, 43-45, 49-53, 57, 60-61, 63, 66-67), the detector (16, 30, 32, 1 16, 130, 216, 
230) producing the prompt gamma ray data (20, 120, 220); 

calculating (72; Fig. 3; H25, 42-46, 51, 70) positron lifetime data (H25, 27, 42-46, 
51, 53, 70) from the positron annihilation data (22, 122, 222) and the prompt gamma ray data 
(20, 120, 220); and 
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continuing (238; Fig. 9; ^75-76) to collect positron annihilation data (22, 122, 222) 
and prompt gamma ray data (20, 120, 220) and calculate (72) positron lifetime data to 
measure an accumulation of lattice damage over time. 

12. A method for evaluating a material specimen (12, 112, 212; Figs. 1, 7-9; ffi|23-31, 33- 
34, 37-40, 43, 47, 50-53, 055-57, 59-63, 66, 68-69, 71-75, 77), comprising: 

mounting (228, Fig. 9; ffl[37, 59-61, 72-74, 77) a neutron source (14, 54, 1 14, 214, 
254; Figs. 1, 7-9; ffl[23-30, 32-34, 36-37, 50, 52, 56-63, 72-74, 77) adjacent the material 
specimen (12, 112, 212); 

mounting (Fig. 9; ffl[ 38, 40, 66-69) a detector (16, 30, 32, 1 16, 130, 216, 230; Figs. 1, 
4, 6-9; 1ffl23-24, 36, 38-42, 44-46, 49, 52-54, 56, 60, 65-72, 74-77) adjacent the material 
specimen (12, 112,212); 

bombarding (232; Fig. 9; U1|23-24, 30-31, 34, 37, 50, 56, 59-61,73) the material 
specimen (12, 1 12, 212) with neutrons (Figs. 1, 7-9; ffl[30-34, 37, 50, 52, 56-63, 73) from the 
neutron source (14, 54, 1 14, 214, 254) to create prompt gamma rays (Figs. 1, 7-9; f1f23-24, 
27, 30, 34, 37-39, 41, 43-45, 49-53, 57, 60-61, 63, 66-67) within the material specimen (12, 
1 12, 212), some of the prompt gamma rays being emitted from the material specimen (12, 
1 12, 212), some of the prompt gamma rays resulting in the formation of positrons (Figs. 1, 7- 
9; H123-25, 27, 30-31, 34, 37, 43, 50-52, 56-57, 60-61, 63) within the material specimen (12, 
1 12, 212) by pair production (18, 1 18, 218; Figs. 1, 7-9; 1Jf23, 28, 30-31, 37, 50, 52, 57, 60- 
61,63); 

collecting (70, 82, 234; Figs. 3, 6, 9; f75) positron annihilation data (22, 122, 222; 
Figs. 1, 3, 4, 6-9; fflj 38, 42-45, 49, 53-54, 67, 70-71, 75-77) by detecting with the detector 
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(16; 30, 32, 1 16, 130, 216, 230) a plurality of emitted annihilation gamma rays (Figs. 1, 7-9); 
1HI23-24, 30, 38-39, 41, 43, 45, 51-53, 57, 66-67) resulting from the annihilation flfl[24, 30, 
43-44, 57) of positrons, the detector (16, 30, 32, 116, 130, 216, 230) producing the positron 
annihilation data (22, 122, 222); 

storing flffi 71, 76) the positron annihilation data (22, 122, 222) on a data storage 
system (225; Fig. 8; fflf 71, 76) for later retrieval and processing; and 

continuing to collect and store positron annihilation data (22, 122, 222), the continued 
collected and stored positron annihilation data (22, 122, 222) being indicative of an 
accumulation of lattice damage flffl24-26, 43, 47, 60-61, 67-69, 71, 74, 77) over time. 

13. The method of claim 12, further comprising: 

collecting prompt gamma ray data (20, 120, 220) by detecting with the detector (16, 
30, 32, 116, 130, 216, 230) a plurality of emitted prompt gamma rays, the detector (16, 30, 32, 
116, 130, 216, 230) producing the prompt gamma ray data (20, 120, 220) ; 

storing prompt gamma ray data (20, 120, 220) on the data storage system (225) for 
later retrieval and processing; and 

continuing to collect and store prompt gamma ray data (20, 120, 220) , the continued 
collected and stored prompt gamma ray data (20, 120, 220) being indicative of an 
accumulation of lattice damage over time. 

14. The method of claim 12, wherein said mounting (228) a neutron source (14, 54, 1 14, 
214, 254) adjacent the material specimen (12, 1 12, 212) comprises mounting (228) the 
neutron source (14, 54, 1 14, 214, 254) to the material specimen (12, 1 12, 212). 
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15. The method of claim 14, wherein said mounting a detector (16, 30, 32, 116, 130, 216, 
230) adjacent the material specimen (12, 112, 212) comprises mounting the detector (16, 30, 
32, 116, 130, 216, 230) to the material specimen (12, 1 12, 212). 

16. The method of claim 15, further comprising positioning a shield (56, 256; Figs. 1, 8; 
ffl|33, 62) adjacent the neutron source (14, 54, 1 14, 214, 254) to absorb stray neutrons. 

17. The method of claim 16, further comprising positioning a moderator (60, 260; Figs. 1, 
8; 1fl|34-36, 63-65) between the neutron source (14, 54, 1 14, 214, 254) and the material 
specimen (12, 112, 212). 

18. The method of claim 12, wherein mounting (228) a neutron source (14, 54, 1 14, 214, 
254) adjacent the material specimen (12, 1 12, 212) comprises mounting (228) an isotopic 
neutron source (54, 254) adjacent the material specimen (12, 1 12, 212). 

19. The method of claim 12, wherein continuing to collect and store positron annihilation 
data (22, 122, 222) is performed while the material specimen (12, 1 12, 212) is in service. 

2 1 . The method of claim 1 2, further comprising: 

retrieving stored positron annihilation data (22, 122, 222); and 

processing (84, 236; Figs. 6, 9; ffl|26, 42, 47, 49, 53-54, 67, 70, 75, 77) the positron 

annihilation data (22, 122, 222) in accordance with a Doppler-broadening algorithm (40; Figs. 
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2, 6; 26-27, 38, 42, 47-49, 51, 70, 75) to produce output data indicative of an accumulation 
of lattice damage over time. 

23. The method of claim 12, further comprising removing the neutron source (14, 54, 1 14, 
214, 254) before collecting positron annihilation data (22, 122, 222). 
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GROUNDS OF REJECTION TO BE REVIEWED ON APPEAL 

1. Whether claims 1-10, 12-19 and 21-23 are patentable under 35 U.S.C. §112, first 
paragraph, as complying with the enablement requirement. 

2. Whether claims 12-19 and 21-23 are unpatentable under 35 U.S.C. §112, first 
paragraph, as not enabled for failing to include critical elements. 

3. Whether claims 1-10, 12-19 and 21-23 are patentable under 35 U.S.C. §112, 
second paragraph, as distinctly claiming the subject matter of the invention. 

4. Whether claims 12-19 and 21-23 are unpatentable under 35 U.S.C. §112, second 
paragraph, for failing to include critical elements. 
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ARGUMENT 

Opening Statement 

The Board will undoubtedly experience a profound sense of deja-vu in this appeal. It 

once again raises many of the same issues that the Board has already decided, yet must decide 

again here (as well as, in the appeal of application serial no. 10/383,096). Once again, the 

examiner has rejected all the claims in this application for indefiniteness and lack of 

enablement. Yet, in his quest to brand all the claims indefinite and not enabled, he has based 

his rejections primarily on limitations that are not part of the claims at issue in this 

application. However, even where the examiner's rejections were actually tied to claim 

limitations requiring use of a Doppler-broadening algorithm, for example, those claims are 

enabled and sufficiently definite. The Board reached the same conclusion in a related case, 

reversing this examiner's rejections of lack of enablement and indefiniteness. The Board's 

words apply equally well here: 

[W]hile the appellant's disclosure does not convey much detail as 
to the various algorithms described therein, these algorithms 
admittedly are generally known in the art. The examiner has failed 
to advance any cogent reasoning as to why the disclosure would 
not have enabled a person of ordinary skill in the art to employ 
these algorithms without undue experimentation. 

Ex parte Akers, Appeal No. 2005-0855, at 8 (Bd. Pat. App. & Interf. 2005) (slip op.). The 

Board also reversed the examiner's rejections based on indefiniteness of the claims as follows: 

The examiner views the appealed claims to be indefinite for 
reasons essentially similar to those listed above in connection with 
the enablement rejection. In this regard, the examiner points to the 
claim limitations relation to . . . the positron lifetime, Doppler- 
broadening, three-dimensional imaging and activation/analysis 

process algorithms [T]he examiner's criticism of the claim 

limitations pertaining to the algorithms . . . focus[es]on the breadth 
of these limitations. It is well settled, however, that mere breadth 
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does not equate to indefiniteness. In re Miller , 441 F.2d 689, 169 
USPQ 597, 600 (CCPA 1971). 

Id. at 9. 

The examiner's new grounds of rejection do not demonstrate that the claims are 
unpatentable either. The examiner argues that the claims are not enabled and incomplete 
because they fail to include a "critical" processing step. No part of the specification describes 
a processing step as "key" or "critical" or "essential." The fact that the specification discloses 
a step of processing does not render that step essential. Moreover, in making these rejections, 
the examiner has sidestepped portions of the specification that describe embodiments for 
collecting and storing data without processing, thus demonstrating that processing is not 
essential. 

Although steadfast in his arguments, the examiner's rejections lack what is required to 
sustain them — evidence from the specification, including the abstract and the claims 
originally filed, to show that any of these claims are indefinite or not enabled. Thus, the 
examiner's rejections must be reversed. 

I. Whether claims 1-10, 12-19 and 21-23 are patentable under 35 U.S.C. §112, first 
paragraph, as complying with the enablement requirement. 

A. Legal Standard For Rejecting Claims Under 
35 U.S.C. §112. First Paragraph 

The legal standard for determining whether the disclosure provides a sufficient 

description of the invention is whether a person reasonably skilled in the art could make or 

use the invention without undue experimentation based on the disclosure and on information 

known in the art. United States v. Telectronics, Inc., 857 F.2d 778, 8 USPQ2d 1217 (Fed. Cir. 

1988). The fact that experimentation may be complex does not necessarily make it undue if 
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the art typically engages in such experimentation. In re Wands, 858 F.2d 731, 8 USPQ2d 
1400 (Fed. Cir. 1988). See Ex parte Lemelson, 2002 WL 32334419, *3 (Bd. Pat. App. & 
Interf. 2002) ("The fact that an invention is sophisticated and complex does not, by itself, lead 
to the conclusion that undue experimentation would be required to make and use the 
invention.") That is, the test of enablement is not whether any experimentation is required, 
but whether, if experimentation is necessary, it is undue. In re Angstadt, 537 F.2d 498, 190 
USPQ214(CCPA 1976). 

The factors to be considered when determining whether a claim limitation is enabled 
and whether any necessary experimentation is "undue" include, but are not limited to: (1) the 
breadth of the claims; (2) the nature of the invention; (3) the state of the prior art; (4) the level 
of one of ordinary skill; (5) the level of predictability in the art; (6) the amount of direction 
provided by the inventor; (7) the existence of working examples; and (8) the quantity of 
experimentation needed to make or use the invention based on the content of the disclosure. 
In re Wands, 858 F.2d at 737; MPEP 2164.01(a). It is improper to conclude that a disclosure 
is not enabling based on an analysis of only one of the above factors while ignoring one or 
more of the others. See, e.g., Ex parte Noguchi, 2002 WL 1801470 (Bd. Pat. App. & Interf. 
2002) (supplying information directed to one factor is insufficient); Ex parte Kopetzki, 2001 
WL 11197754 (Bd. Pat. App. & Interf. 2001) (merely reciting facts without reference to the 
factors is insufficient); MPEP 2164.01(a). Moreover, whether the specification is enabling 
must be determined with reference only to the limitations claimed. CFMT, Inc. v. Yieldup 
Int'l Corp., 349 F.3d 1333, 1338 (Fed. Cir. 2003). 

With regard to the burden of proof required to support a rejection under §112, the 
Patent Office is required to assume that the specification complies with the enablement 
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provision of §112 unless it has acceptable evidence or reasoning to suggest otherwise. See, 
e.g., In re Marzocchi, 439 F.2d 220, 169 USPQ 367 (CCPA 1979). The Patent Office thus 
must provide reasons, supported by the record as a whole, why the specification is not 
enabling. Then and only then does the burden shift to the applicant to show that one of 
ordinary skill in the art could have practiced the claimed invention without undue 
experimentation. Gould v. Missinghojf, 229 USPQ 1 (D.D.C. 1985), aff'd in part, vacated in 
part, and remanded sub. nom., Gould v. Quigg, 822 F.2d 1074, 3 USPQ2d 1302 (Fed. Cir. 
1987). Mere conclusionary statements as to the level of ordinary skill in the art are not a 
sufficient basis for a rejection under 35 U.S.C. §112. In re Brebner, 455 F.2d 1402, 173 
USPQ 169 (CCPA 1972). In addition, the law does not require, and indeed prefers, that a 
patent specification omit that which is well-known. In re Buchner, 929 F.2d 660, 18 USPQ2d 
1331 (Fed. Cir. 1991). 

B. Response to Examiner's Rejections 

The examiner failed to carry his prima facie burden of demonstrating that claims 1-10, 
12-19 and 21-23 are not enabled as required by §112, first paragraph. First, the examiner's 
rejections are based on limitations that are not part of the claims. Second, where the examiner 
did base his rejections on at least one limitation actually contained in the claims, the examiner 
failed to analyze the enablement issue with reference to the eight required factors. Finally, 
given that Doppler-broadening (and other) algorithms are well-known in the art, the level of 
detail contained in the specification sufficiently enables the claims. 

1. Rejections are improperly based on limitations not in the claims . 

As an initial matter, the examiner argues that "the claims require the use of algorithms 
(i.e., Doppler-broadening and positron lifetime), detector(s), a neutron source and a data 
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processor. Thus to enable an artisan to make and use an operative embodiment, one must first 
know which particular algorithm(s) to use, and how and in what manner the data from the 
detectors are to be analyzed in the 'data processor'." Office Action (1/16/07), pp. 6-7; Final 
Office Action (3/1/06), p. 2. A fundamental deficiency in the examiner's argument is that 
none of the claims requires use of a data processor. And, none of the claims requires use of a 
positron lifetime algorithm. Since the examiner's rejections are based on limitations not 
found in the claims, they should be reversed. 

a. Independent claim 12 and dependent claims 13-19 and 23 

Independent claim 12 and dependent claims 13-19 and 23 require neither processing in 
accordance with an algorithm nor use of a data processor. Therefore, the examiner's §112 
rejections of these claims for lack of enablement of those limitations are improper on their 
face and should be reversed on that ground alone. 

b. Independent claim 1 and dependent claims 2 and 21 

Now the appellant turns to independent claim 1 and dependent claim 21, both of which 
contain the limitation, "processing ... in accordance with a Doppler-broadening algorithm." 
Claim 2 claims the method of claim 1 further comprising "calculating positron lifetime data;" 
however, claim 2 does not require use of a positron lifetime algorithm. The examiner argues 
that the claims are not enabled because undue experimentation is necessary to select 
algorithms (including Doppler-broadening, positron lifetime, 3D imaging and selective 
activation algorithms) that will (1) obtain accurate detection of lattice defects (i.e., not too few 
and not too many); (2) obtain reasonable statistical analyses; (3) determine which algorithm 
(in conjunction with selection of "a suitable data processor") "actually provides a reasonable 
determination of the presence of enough lattice defects to be of concern;" and (4) determine 
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"whether there are present enough lattice defects to cause metal fatigue, etc " Office 

Action (1/16/07), p. 10; Final Office Action (3/1/06), p. 6. The examiner further argues that 
the appellant has "refused to disclose how these algorithms have to be combined," referring to 
Doppler-broadening and positron lifetime algorithms. Office Action (1/16/07), p. 11; Final 
Office Action (3/1/06), p. 7 (underlining omitted). However, these limitations are not in the 
claims and, therefore cannot be considered when determining whether the specification 
enables the claims. 

None of these claims requires a data processor; positron lifetime, 3-D imaging or 
selective activation algorithms; accuracy; reasonable statistical analysis; "provid[ing] a 
reasonable determination of the presence of enough lattice defects to be of concern;" or 
determining "whether there are present enough lattice defects to cause metal fatigue." Office 
Action (1/16/07), p. 10; Final Office Action (3/1/06), p. 6. Similarly, none of the claims 
requires combining Doppler-broadening and positron lifetime algorithms. As the Federal 
Circuit has held, "Title 35 does not require that a patent disclosure enable ... a perfected, 
commercially viable embodiment absent a claim limitation to that effect." CFMT f Inc. v. 
Yieldup Int'l Corp., 349 F.3d 1333, 1338 (Fed. Cir. 2003). Therefore, the examiner cannot 
evaluate enablement based on limitations not contained in the claims. Id., 349 F.3d at 1338 
(since claims contained no standard for cleanliness, claims for cleaning methods were enabled 
for any level of contaminant removal); Ex parte Lemelson, 2002 WL 32334419, *4. 

Ignoring the plain language of the claims, the examiner nonetheless insists that the 
claims require the use of algorithms, a detector, a neutron source and a data processor. Office 
Action (1/16/07), p. 16. The examiner makes this statement based on his belief that 
"Appellant does not disclose any alternative method of performing the material evaluation 
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that does not require these four elements." Office Action (1/16/07). However, that is not the 

case. Paragraph 71 of the specification describes "the third embodiment" wherein "the data 

processing system serves merely as a data collection device": 

[0071] Alternatively, the third embodiment 210 to be used in a 
long-term arrangement, such as, for example to measure the 
build-up or accumulation of lattice defects during the service 
life (or some portion of the service life) of the material 
specimen 212, the data processing system 224 could be 
simplified considerably to where the data processing system 
224 serves merely as a data collection device, collecting the 
data (e.g., the prompt gamma ray data 220 and/or the positron 
annihilation data 222) produced by the detector assembly 216, 
and storing the data in a data storage system 225 for later 
retrieval and processing in accordance with the teachings 
provided herein. . . . Thereafter, the collected and stored data 
may be retrieved (i.e., downloaded) and processed in 
accordance with the teaching herein in order to provide 
information indicative of the build-up or accumulation of 
lattice defects during the collection period. 

See infra at pp. 29-30 (quoting the abstract which describes a similar embodiment). Contrary 

to the examiner's arguments, no processing using an algorithm is disclosed. 

With respect to his rejections of claim 2, the examiner argues that a data processor and 

a positron lifetime algorithm must be included in claim 2 as essential elements; otherwise, 

claim 2 will not be enabled. Office Action (1/16/07). However, this argument cannot be 

sustained in light of passages from the specification stating that "the data processing system 

224 may process the prompt gamma ray data 220 and positron annihilation data 222 in 

accordance with a positron lifetime algorithm . . ." See Specification, ^ 70 (emphasis added). 

The examiner argues that because no other method for calculating positron lifetime is 

disclosed, "there must be other undisclosed means" for calculating positron lifetime, 

insinuating that the appellant has failed to disclose the best mode (which is denied). First, the 

disclosure of only one embodiment for making and using the invention satisfies the 
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enablement requirement. MPEP 2163.01(b). Second, the examiner is not entitled to assume 
that these allegedly "undisclosed means" - to the extent they exist at all (which is denied) - 
are better than that disclosed in the specification, for they just as well could be worse. And, 
finally, the examiner has provided no evidence that the best mode has not been disclosed. 

c. Figure 2 

Similarly, but without reference to any claims, the examiner objects to the sufficiency 
of the disclosure in Figure 2, regarding a feedback arrangement between the data processing 
system 24 and the algorithms 38, 40. Office Action (1/16/07), p. 12; Final Office Action 
(3/1/06), p. 8. However, the examiner has not tied these arguments to the elements of any 
claim to show how such elements are not enabled. CFMT, Inc., 349 F.3d at 1338 (level of 
disclosure necessary for enablement varies with the scope of the claims). Therefore, the 
examiner's arguments in this regard cannot be used to demonstrate a prima facie case of non- 
enablement. 

d. Ex parte Lemelson 

The case of Ex parte Lemelson, 2002 WL 32334419, presents similar facts to the ones 

here and is persuasive in its analysis of enablement. In that case, the invention was directed to 

a method and apparatus for controlling the travel of a powered vehicle, comprising calculating 

whether the vehicle and an object are on a collision course by measuring the distance to and 

relative velocity between the two, and using fuzzy logic interference rules to determine a 

combination of steering and acceleration to avoid a collision. The examiner rejected the 

claims for lack of enablement on grounds almost identical to the ones here: 

The rejection also states that the disclosure is insufficient for 
teaching how to distinguish one object from another object or how 
to determine distance and relative velocity of each object. . . . The 
examiner observes that the object detection and recognition and the 
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fuzzy logic rules require complex algorithms and extensive 
processing that are not sufficiently described in the disclosure 
[citation omitted]. In the final rejection the examiner added that 
the complexity of putting all the features into one real-time system 
required undue experimentation. 

Id., at *2. The appellants argued that fuzzy logic and image analysis were techniques well- 
known in the art. The Board reversed the examiner's rejections. First, the Board held that the 
examiner improperly based his rejections on a real-time operation requirement not included in 
the claims. It also held, "Second, the examiner's findings of undue experimentation are mere 
conclusions based on the examiner's own speculations. The examiner essentially finds that 
there would be undue experimentation because the claimed invention requires sophisticated 
and complex operations. The fact that an invention is sophisticated and complex does not, by 
itself, lead to the conclusion that undue experimentation would be required to make and use 
the invention." Id. at *3. The examiner here argues that the instant case differs from the 
Lemelson case because he did show undue experimentation and he did show that the elements 
associated with the rejection are required (both of which are denied). However, arguing that 
the result of the case is different does not distinguish Ex parte Lemelson from this case on the 
facts. The Board's analysis in the Lemelson case applies equally well here, demonstrating that 
the examiner here failed to carry his burden. 

e. Dependent claims 2-10 and 22 

In addition, claims 2-10 and 22 depend from claim 1 and are also allowable at least for 
the reasons that claim 1 is allowable. 

2. The eight required factors were not addressed . 

While the examiner sets forth a laundry list of variables that he argues one of ordinary 
skill in the art must consider when choosing an appropriate Doppler-broadening algorithm, 
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positron lifetime algorithm, or other algorithm, the examiner has nonetheless failed to meet 
his prima facie burden. As an initial matter, again, (with the exception of Doppler- 
broadening) none of the claims at issue here requires processing in accordance with positron 
lifetime or other algorithms, so those "limitations" may not be considered in determining 
whether the claims are enabled. In addition, the examiner failed to undertake the required 
analyses to demonstrate lack of enablement. 

Only claims 1 and 21 contain a limitation that the examiner argues was not enabled. 
That limitation is "processing ... in accordance with a Doppler-broadening algorithm. As 
stated above, while claim 2 claims "calculating positron lifetime data," it does not require use 
of a positron lifetime algorithm. However, in rejecting claims 1, 2 and 21 (as well as all the 
other claims) for non-enablement, the examiner was required to analyze each of the eight 
required factors set forth in In re Wands, 858 F.2d 731, 737 (Fed. Cir. 1988), with reference to 
the claims, which he did not. Nor did the examiner supply facts to support each and every 
factor required by In re Wands. The Board has repeatedly reversed similar rejections, holding 
that the examiner failed to meet his prima facie burden. See, e.g., Ex parte Noguchi, 2002 
WL 1801470 (Bd. Pat. App. & Interf. 2002) (supplying information directed to one factor is 
insufficient); Ex parte Kopetzki, 2001 WL 11197754 (Bd. Pat. App. & Interf. 2001) (merely 
reciting facts without reference to the factors is insufficient). Even assuming that the 
examiner put forth any "facts" (which is denied), the "mere citing of a series of facts, without 
relating them to the criteria for determining undue experimentation set forth in In re Wands, is 
not sufficient to establish a prima facie case of nonenablement." Ex parte Kopetzki, 2001 WL 
1 1 197754, at *3. In addition, by reviewing the variables that one of ordinary skill in the art 
would know to consider, the examiner makes the appellant's point that one of ordinary skill in 
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the art would well know how to pick a suitable Doppler-broadening algorithm (or positron 
lifetime algorithm, assuming one were even required) and select appropriate constants, 
making any experimentation routine and not undue. 
3. The specification is enabling . 

Finally, given that Doppler-broadening algorithms are well-known, the level of detail 
provided in the present specification is sufficiently enabling, as the Board determined with 
reference to the appeal in the divisional case issued as U.S. Patent No. 7,058,153. See supra at 
p. 12. Moreover, even if a positron lifetime algorithm were required by claim 2 (which is 
denied), positron lifetime algorithms are also well-known, rendering the level of detail 
provided in the specification enabling as to that limitation, too. 

The appellant specifically states in the specification that the Doppler-broadening 
algorithms and positron lifetime algorithms are well-known in the art and could be readily 
implemented by persons having ordinary skill in the art after having become familiar with the 
teachings of the present invention. Specification, fflf [0043]-[0049], Because these algorithms 
are well-known to those of ordinary skill in the art who are familiar with their applications, 
the specification need not and does not describe them further. It is undisputed that these 
algorithms are well-known. Office Action (1/16/07), p. 18. Since that fact is undisputed, the 
examiner was bound to establish that the amount of experimentation in using these algorithms 
familiar to those of ordinary skill in the art would be undue. As explained above, he failed to 
make such a showing because he did not analyze the undue experimentation factors from In re 
Wands. 

During prosecution, the appellant provided examples of issued patents to show the 
general level of detail required to enable similar inventions for those of ordinary skill in the 
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art. See Amendment (1/13/06), p. 12. In his Summary of Argument presented in the 
appellant's Amendment dated January 13, 2006, the appellant stated that "the best measure of 
level of disclosure required in this particular field of endeavor can be obtained by reviewing 
issued U.S. patents in the same field." Amendment (1/13/06), p. 7. The examiner argued that 
U.S. Patent No. 6,178,218 to Akers ("Akers") and U.S. Patent No. 4,064,438 to Alex, et al. 
("Alex") cannot be used by analogy to show the appropriate level of detail required for 
enablement. Office Action (1/16/07), p. 13; Final Office Action (3/1/06), p. 9. However, the 
examiner's arguments are not well-founded. 

Specifically, the examiner argued that Akers cannot be used because his Doppler- 
broadening algorithm "is not modified by feed back from the data processing system. . . ." 
Office Action (1/16/07), p. 13; Final Office Action (3/1/06), p. 9. However, the claims at 
issue in this appeal do not contain as limitations either use of a data processing system or 
modification of the Doppler-broadening algorithm through feedback. Basing non-enablement 
rejections on limitations not found in the claims is improper. 

Further, the examiner argues that Alex cannot be used to demonstrate enablement by 
analogy because Alex uses a positron lifetime algorithm for a different purpose than do the 
claims of the present invention. As an initial matter, none of the claims at issue here requires 
use of a positron lifetime algorithm. While claim 2 claims "calculating positron lifetime 
data," it does not require use of a positron lifetime algorithm either. In any event, the 
examiner misses the point. The point is that the description in Alex enables his claim 
limitation of a method "wherein the variations in annihilation characteristics within the tested 
material are detected by positron lifetime measurements." The description of positron 
lifetime measurements in Alex contains even less detail than the present specification, merely 
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referring to one page of a previously published book and one issued patent. Compare Alex, 
Col. 6, Ins. 60-66 with Specification fflf [0043]-[0046]. If such a brief description in Alex 
were considered enabling, then the more detailed description of the present specification must 
also be considered enabling to one of ordinary skill in the art. Finally, the other references 
previously cited by the examiner (i.e., Derlet, Bandzuch, Zhu, and Shaffer) illustrate that 
persons having ordinary skill in the art are well-aware of the appropriate use and description 
of Doppler-broadening and positron lifetime algorithms, as well as the various issues and 
complexities involved in selecting, modifying and combining any of these algorithms, 
rendering any experimentation required routine and not undue. 

The examiner argues that "the specification indicates that there are different known 
Doppler-broadening algorithms but, fails to disclose which of these known algorithms would 
be suitable for use in the present invention." Office Action (1/16/07), p. 7; Final Office 
Action (3/1/06), p. 3. That is not the case. The specification specifically discloses a suitable 
algorithm when it states that "in one preferred embodiment of the invention, the Doppler- 
broadening algorithm 40 may comprise the Doppler-broadening algorithm described in U.S. 
Patent No. 6,178,218 Bl . . . Specification, f[0048]. Contrary to the examiner's suggestion, 
there is no experimentation involved at all in selecting the precise algorithm the appellant 
identified in the specification. 

The examiner reads the specification as requiring that " both the positron lifetime 
algorithm and Doppler-broadening algorithm, each requires two detectors ." Office Action 
(1/16/07), p. 12; Final Office Action (3/1/06), p. 8. He then states the specification also 
discloses that one detector may be used. Without reference to any limitations in the claims, 
the examiner argues that appellant's "disclosure is thus clearly insufficient and non-enabling 
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as to which specific known positron lifetime algorithm can be used with only a single detector 
and, which specific Doppler-broadening algorithm can be used with only a single detector. " 
Office Action (1/16/07), p. 12; Final Office Action (3/1/06), p. 8. The appellant now 
understands that "the examiner is not arguing the number of detectors that are enabled by the 
specification. The issue is the lack of enablement for the specific algorithm to apply when 
using only a single detector, and the undue experimentation resulting therefrom." Office 
Action (1/16/07), p. 19. 

Independent claims 1 and 12 claim a "method for evaluating a material specimen, 
comprising ... a detector . . . However, neither one of these claims contains as a limitation 
use of a specific algorithm. Thus, once again, the examiner has premised his non- 
enablement rejections on limitations that are not part of the claims. That is improper and the 
rejections must be reversed for that reason alone. CFMT f Inc. v. Yieldup Int'l Corp., 349 F.3d 
1333, 1338 (Fed. Cir. 2003). (since claims contained no standard for cleanliness, claims for 
cleaning methods were enabled for any level of contaminant removal); Ex parte Lemelson, 
2002 WL 32334419, *4. See supra at pp. 16-21. 

In addition, the examiner was required to analyze each of the eight required factors set 
forth in In re Wands, 858 F.2d 731, 737 (Fed. Cir. 1988), with reference to the claims, which 
he did not. Nor did the examiner supply facts to support each and every factor required by In 
re Wands. Even assuming that the examiner had recited any facts (which is denied), merely 
reciting facts without reference to each of the eight factors is insufficient to sustain a rejection 
for non-enablement. Ex parte Kopetzki, 2001 WL 1 1 197754 (Bd. Pat. App. & Interf. 2001). 
See the discussion supra at pp. 21-23. 
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Claims 2-10 and 22 ultimately depend from claim 1 and are also allowable at least for 
the reasons that claim 1 is allowable. Claims 13-19, 21 and 23 ultimately depend from claim 
12 and are also allowable at least for the reasons that claim 12 is allowable. 

Having failed (1) to base his rejections on limitations in the claims and (2) to analyze 
the In re Wands factors, the examiner failed to meet his prima facie burden to demonstrate 
that the claims are not enabled. In any event, as the Board previously determined with 
reference to a related case, descriptions of the various algorithms used in detecting positron 
annihilation events are sufficiently enabling. 

II. Whether claims 12-19 and 21-23 are unpatentable under 35 U.S.C. §112, first 
paragraph, as not enabled for failing to include critical elements. 

A. Legal Standard For Rejection Claims 
Under 35 U.S.C. §112. First Paragraph 

The legal standard set forth above at pages 14-16 is incorporated herein by reference. 

In addition, a claim that omits matter disclosed to be essential in the specification may 
be rejected under §112, first paragraph, as not enabled. In re Mayhew, 527 F.2d 1229, 1233 
(C.C.P.A. 1976); MPEP 2164.08. However, the specification as a whole must disclose the 
omitted element or step to be critical to the claimed invention. MPEP 2164.08. As stated in 
the MPEP, "Features which are merely preferred are not considered to be critical. In re Goffe, 
542 F.2d 564, 567, 191 USPQ 429, 431 (CCPA 1976)." MPEP 2164.08. 

As the MPEP emphasizes, "[A]n enablement rejection based on the grounds that a 
disclosed critical limitation is missing from a claim should be made only when the language 
of the specification makes it clear that the limitation is critical for the invention to function as 
intended. Broad language in the disclosure, including the abstract, omitting an allegedly 
critical feature, tends to rebut the argument of criticality." MPEP 2164.08. See also Ex parte 
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Kato, et al, 2002 WL 465390 (Bd. Pat. App. & Interf. 2002) (specification must show the 
omitted limitation to be critical); Ex parte Nakano, et al. 9 1998 WL 1735647 (Bd. Pat. App. 
& Interf. 1998) (same). 

B. Response to Examiner's Rejections 

The examiner rejected claims 12-19 and 21-23 as not enabled based on his belief that 
the claims improperly omit the "processing of positron annihilation data by Doppler- 
broadening algorithm and/or positron lifetime algorithm, which are critical or essential to the 
practice of the invention." Office Action (1/16/07), p. 3. However, in making his rejection, 
the examiner failed to analyze the lack of enablement factors as required by In re Wands and 
he failed to consider the specification as a whole, including disclosures omitting the allegedly 
critical feature. 

a. Claim 21 

As an initial matter, the appellant notes that claim 21 contains the limitation 
"processing the positron annihilation data in accordance with a Doppler-broadening algorithm 
to produce output data indicative of an accumulation of lattice damage over time." Since 
claim 21 already includes the precise limitation that the examiner is arguing needs to be 
included in the claims as critical, the examiner's rejection of claim 21 must be reversed. 

b. Claims 12-19,22-23 

Turning to the remainder of the claims here rejected, the examiner failed to address the 
In re Wands factors to demonstrate lack of enablement, as required. See supra at pp. 14-16. 
In addition, there is no language in the specification which describes as "critical" or 
"essential" a step of processing data by Doppler broadening or positron lifetime algorithm. In 
addition, as described more fully below, the specification (in portions not cited by the 
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examiner), including the abstract, describe embodiments of the invention that do not include a 
step for processing data by means of a Doppler broadening or positron lifetime algorithm: 

1 . The abstract discloses an embodiment of the invention that omits a step 

of using an algorithm to process data, as follows: 

A method for evaluating a material specimen comprises: 
Mounting a neutron source and a detector adjacent the material 
specimen; bombarding the material specimen with neutrons 
from the neutron source to create prompt gamma rays within 
the material specimen, some of the prompt gamma rays being 
emitted from the material specimen, some of the prompt 
gamma rays resulting in the formation of positrons within the 
material specimen by pair production; collecting positron 
annihilation data by detecting with the detector at least one 
emitted annihilation gamma ray resulting from the annihilation 
of a positron; storing the positron annihilation data on a data 
storage system for later retrieval and processing; and 
continuing to collect and store positron annihilation data, the 
continued collected and stored positron annihilation data being 
indicative of an accumulation of lattice damage over time. 

2. Although the cited example from the abstract demonstrates that the 

claims at issue are enabled, other embodiments in the specification demonstrate that a step of 

processing data using an algorithm is not critical either. For example, paragraph 71 discloses 

an embodiment of the invention for collecting and storing the data, without using an 

algorithm to process the data, as follows: 

[0071] Alternatively, the third embodiment 210 to be used in a 
long-term arrangement, such as, for example to measure the 
build-up or accumulation of lattice defects during the service 
life (or some portion of the service life) of the material 
specimen 212, the data processing system 224 could be 
simplified considerably to where the data processing system 
224 serves merely as a data collection device, collecting the 
data (e.g., the prompt gamma ray data 220 and/or the positron 
annihilation data 222) produced by the detector assembly 216, 
and storing the data in a data story system 225 for later 
retrieval and processing in accordance with the teachings 
provided herein. . . . Thereafter, the collected and stored data 
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may be retrieved (i.e., downloaded) and processed in 
accordance with the teaching herein in order to provide 
information indicative of the build-up or accumulation of 
lattice defects during the collection period. 

(Emphasis added). 

In addition, paragraph 70 of the specification notes that such 

algorithms may (or may not) be used: 

[0070] The data processing system 224 is operatively 
associated with the detector system 216 and receives the 
prompt gamma ray data 220 and positron annihilation data 222 
produced by the detector system 216. As was briefly described 
above, the data processing system 224 may process the prompt 
gamma ray data 220 and positron annihilation data 222 in 
accordance with a positron lifetime algorithm (e.g., positron 
lifetime algorithm 38) . . . [s]o processing the prompt gamma 
ray data 220 and the positron annihilation data 222 results in 
positron lifetime data. In addition, the data processing system 
224 may also process the positron annihilation data 222 in 
accordance with a Doppler-broadening algorithm (e.g., 
Doppler-broadening algorithm 40) in the manner already 
described for the first embodiment 10. 

(Emphasis added). 

3. The appellant further notes that original claims 12-19 and 22-23 did not 
contain the allegedly omitted steps. See the discussion infra at pp. 39-40. 

4. As is the case here, where the examiner has failed to provide 
unequivocal evidence from the specification that the allegedly omitted step was described as 
essential, the Board has reversed rejections under §112, first paragraph. See Ex parte Kato, 
2002 WL 465390 at *5 ("examiner failed to cite any passage of the specification or other 
statements of record that would establish that any essential element as been omitted from the 
claim under appeal"). While the examiner here selectively cited from the specification 
passages to show that algorithms were disclosed in certain embodiments, that disclosure by 
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itself is not sufficient to establish that critical elements were omitted from the claims, 
especially where other passages from the specification, including the abstract, describe 
embodiments without such allegedly omitted elements. See Ex parte Nakano, 1998 WL 
1735647 at *4 ("The mere fact that other steps have been disclosed in the preferred 
embodiment does not render each and every step thereof an essential step.") 

The examiner's rejections of claims 12-19 and 21-23 under §112, first paragraph, 
must be reversed. 

111. Whether claims 1-10, 12-19 and 21-23 are patentable under 35 U.S.C. §112, 
second paragraph, as distinctly claiming the subject matter of the invention. 

A. Legal Standard For Rejecting Claims 

Under 35 U.S.C. §112, Second Paragraph 

Section 112, second paragraph, requires that the claims particularly point out and 

distinctly claim what the applicant regards as his invention. 35 U.S.C. § 112. The standard is 

one of reasonably particularity and completeness. MPEP 2173.02. "The indefiniteness 

inquiry focuses on whether those skilled in the art would understand the scope of the claim 

when the claim is read in light of the rest of the specification." Union Pacific Res. Co. v. 

Chesapeake Energy Corp., 236 F.3d 684, 691 (Fed. Cir. 2001); Credle v. Bond, 25 F.3d 1566, 

1576 (Fed. Cir. 1994) (all that is required under § 1 12, Tf2 is that one of ordinary skill in the 

art would understand the claim language based on the specification and the drawings). The 

degree of precision necessary for adequate claims depends on the nature of the subject matter. 

Miles Laboratories, Inc., v. Shandon, Inc., 27 USPQ2d 1123 (Fed. Cir. 1993). As such, the 

Patent Office does not permit per se rejections of certain language as being indefinite under § 

112, If 2. MPEP 2173.02 ("Office policy is not to employ per se rules to make technical 
rejections.") 
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The examiner bears the initial burden of presenting evidence or reasoning to explain 
why persons skilled in the art would not understand the claim language based on the 
specification. See In re Oetiker, 977 F.2d 1443, 1445 (Fed. Cir. 1992). To meet that burden 
of proof, the MPEP requires that "[d]efiniteness of claim language must be analyzed, not in a 
vacuum, but in light of: 

(A) The content of the particular application disclosure; 

(B) The teachings of the prior art; and 

(C) The claim interpretation that would be given by one possessing the ordinary level 
of skill in the pertinent art at the time the invention was made." MPEP 2173.02 (emphasis 
added). See In re Johnson, 558 F.2d 1008, 1015, 194 USPQ 187 (CCPA 1977). In addition, 
the Federal Circuit has repeatedly held that the enablement requirement and the written 
description requirement of 112 are not to be conflated. "We interpret 35 U.S.C. §112, fl to 
require a written description requirement separate and apart from the enablement 
requirement." In re Curtis, 354 F.3d 1347, 1357 (Fed. Cir. 2004). See In re Hyatt, 708 F.2d 
712, 715 (Fed. Cir. 1983) (the fact that a claim is not enabled does not render it indefinite). 
"After evidence or argument is submitted by the applicant in response, patentability is 
determined on the totality of the record with due consideration to persuasiveness of 
argument." In re Oetiker, 977 F.2d at 1445. 

B. Response to Examiner's Rejections 

The examiner never established a prima facie case of indefiniteness. Indeed, the 
examiner's rejections here are indefinite, not the appellant's claim language. In short, by 
confusing breadth with indefiniteness, the examiner has issued improper rejections under 
§112, second paragraph. Moreover, even assuming that the examiner did carry his initial 
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burden (which is denied), the appellant demonstrated that the claim language is definite, 
which evidence the examiner ignored. The indefiniteness rejections should be reversed. 
1. The only specific rejection was traversed . 

The only specific claim language that the examiner rejected during prosecution was 
u at least one" emitted gamma ray, which he rejected as indefinite in his office action of May 
10, 2005. Claims previously containing that limitation have all been amended to "a plurality 
of emitted gamma rays, thereby traversing the examiner's indefiniteness rejection. The 
examiner acknowledged that this rejection had been traversed in his October 17, 2005 office 
action. Despite this acknowledgement, in the same October 17, 2005 office action, the 
examiner newly rejected all of the claims as "vague and indefinite" without even specifying 
which claim language would not be understood by one of ordinary skill in the art. In the final 
office action, the examiner summarily dismissed the appellant's argument and evidence 
regarding definiteness as "not convincing" and merely incorporated by reference his previous 
action of October 17, 2005 which incorporated his previous action of May 10, 2005 which 
incorporated his arguments directed to the enablement issue. Similarly, in his non-final office 
action of January 16, 2007, the examiner rejected the claims as "vague and indefinite" again, 
with reference to the same October 17, 2005 office action. The examiner's only response to 
the appellant's intial appeal brief was to assert that "appellant's arguments are: a) simply 
repetitions of items that have been adequately addressed above." Office Action (1/16/07), p. 
19. 

Thus, when reduced to its essence, prosecution in this case has produced only one 
specific rejection for indefiniteness. The only specific rejection having been traversed, the 
examiner's rejection for indefiniteness should be reversed on those grounds alone. 
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2. The examiner failed to meet his prima facie burden . 

Once his initial indefiniteness rejection was traversed, the examiner never specifically 
and directly identified the claim language that he considered "vague and indefinite." Rather, 
as stated above, he just bootstrapped his indefiniteness rejections onto his lack of enablement 
rejections. The examiner's indefiniteness rejections are therefore improper. 

Not only did he fail to identify any deficient claim language, but also the examiner 
failed to perform the indefiniteness analyses required by MPEP 2173.02. For example, since 
the examiner never identified the claim language at issue, he certainly did not interpret the 
claim language, as required by the MPEP. Further the examiner has not provided any 
persuasive reasoning why any specific claim language is vague and indefinite. 

Even assuming that such language were specified, by relying only on arguments he 
made to support his non-enablement rejections to support his indefiniteness rejections, the 
examiner failed to carry his burden. Indeed, the Federal Circuit has specifically stated, "[I]f 
the 'enabling' disclosure of a specification is not commensurate in scope with the subject 
matter encompassed by a claim, that fact does not render the claim imprecise or indefinite or 
otherwise not in compliance with the second paragraph of §1 12; rather the claim based on an 
insufficient disclosure (§112, first paragraph) and should be rejected on that ground." In re 
Hyatt.lOS F.2dat715. 

Having been based entirely on his non-enablement rejections, appellant therefore 
believes that the examiner's indefiniteness rejections extend to the scope of the Doppler- 
broadening algorithm claim limitation of claims 1 and 21, as well as claims 2-10 and 22 
which depend from claim 1. (The appellant has already explained the error in the examiner's 
non-enablement rejections as improperly based on limitations not found in the claims. See 
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supra at pp. 16-21). However, the examiner's rejections of all the claims as indefinite cannot 
be sustained on that basis. 

a. Independent claim 12 and dependent claims 13-19 and 23 

First, claims 12-19 and 23 do not require the use of any algorithm and the examiner 
has put forth no evidence or reasoning as to why any specific claim limitation of those claims 
is imprecise. Having made out no case whatsoever that any language in claims 12-19 and 23 
is indefinite, the examiner's rejections in that regard must be reversed. 

b. Independent claim 1 and dependent claim 21 

Second, with respect to the indefiniteness rejections of claim 1 (and its dependent 
claims 2-10 and 22) and claim 21 presumably based on the use of a Doppler-broadening 
algorithm, the words of the Board in the related appeal apply just as well here: "[T]he 
examiner's criticisms of the claim limitations pertaining to the algorithms . . . focus on the 
breadth of those limitations. It is well settled, however, that mere breadth does not equate to 
indefiniteness. In re Miller , 441 F.2d 689, 169 USPQ 597, 600 (CCPA 1971). Accordingly, 
the examiner's contention that the appealed claims are indefinite is unpersuasive." Ex parte 
Akers, supra, at 9-10. 

3. The examiner failed to consider evidence of definiteness . 
a. Independent claim 1 and dependent claims 2 and 21 

Although the examiner never identified the specific claim language he considered to 
be indefinite, the applicant surmised that his rejections must have related to the claims 
language "processing ... in accordance with a Doppler-broadening algorithm" (claims 1 and 
21) and "calculating positron lifetime data" (claim 2). (Neither of these terms is present in 
independent claim 12 or claims 13-19 and 23). Based on that understanding, the appellant 
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responded with evidence and argument, explaining why the "various algorithms are not vague 
[or] indefinite . . . i.e., a the person having ordinary skill in the art would understand the 
algorithms and possess the knowledge required to evaluate the many algorithms available and 
select those that would be appropriate for the desired application." Amendment (1/13/06), p. 
13. Such evidence and argument demonstrated that the examiner's rejections did not meet his 
prima facie burden. 

Specifically, the appellant argued in his Amendment (1/13/06), pp. 12-13 that the 
claim language was not indefinite because: 

(i) Persons having ordinary skill in the art would understand the algorithms and 
possess the knowledge required to evaluate the many algorithms available and select those 
that would be appropriate for the desired application. 

(ii) With regard to the Doppler-broadening algorithm, the written description clearly 
states that several different types of Doppler-broadening techniques have been developed and 
could be used. The written description also states that the Doppler-broadening algorithm may 
comprise the algorithm disclosed in U.S. Patent No. 6,178,218. 

(iii) The written description specifically states that systems for detecting positron 
lifetimes, as well as the algorithms utilized thereby, are well-known in the art and could be 
easily provided by persons having ordinary skill in the art after having become familiar with 
the teachings of the present invention. 

Instead of addressing these arguments and presenting evidence why persons of 
ordinary skill in the art would not understand the meaning the claim terms set forth above, the 
examiner pronounced the arguments "not convincing." Final Office Action (3/1/06) at p. 2. 
Without any further explanation, he merely incorporated his previous rejections by reference 
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into his final office action. Final Office Action (3/1/06) at p. 12. No clarification was 
forthcoming in the most recent office action, mailed January 16, 2007, either. In a similar 
situation where, in the examiner's answer, the examiner responded that the appellants' 
arguments were "not persuasive," the Board held that the examiner failed to satisfy his 
burden. See Ex parte Voordum, 2002 WL 1801 166, *4 (Bd. Pat. App. & Interf. 2002). 

Even assuming that the examiner had specified the claim language set forth above as 
indefinite (which is denied), one of ordinary skill in the art would well understand the scope 
of "processing ... in accordance with a Doppler-broadening algorithm" (claims 1 and 21) and 
"calculating positron lifetime data" (claim 2) in light of the specification. The examiner has 
failed to carry his burden to demonstrate otherwise. 

Claims 3-10 and 22 depend from claim 1 and are allowable at least for the reasons that 
claim 1 is allowable. 

b. Independent claim 12 and dependent claims 13-19 and 23 

As mentioned above, claims 12-19 and 23 contain no limitations directed to any 
algorithmic claim limitation that the examiner evidently considers indefinite. Therefore, 
having failed to raise any specter that these claims are indefinite, the examiner's rejections as 
to these claims should be summarily reversed. 

IV. Whether claims 12-19 and 21-23 are unpatentable under 35 U.S.C. §112, second 
paragraph, for failing to include critical elements. 

A. Legal Standard for Rejecting Claims Under 
35 U.S.C. §112, Second Paragraph 

The legal standard for rejecting claims as indefinite under §112, second paragraph is 

set forth above at pages 31-32 and is incorporated herein. 

In addition, a claim which fails to interrelate essential elements of the invention may 
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be rejected as indefinite under §112, second paragraph. MPEP 2172.01. However, the Board 
has cautioned that an "'incompleteness' rejection should be extremely rare . . . Ex parte 
Trummer, 1997 WL 33106615, *4 (Bd. Pat. App. & Interf. 1997). The Board explained that 
this "rejection could be applied in almost every case since it could be said that every part of a 
disclosed combination is somehow essential to the overall purpose of the invention. This 
would infringe an applicant's right to claim what he regards as his invention." Id The 
examiner bears the burden of demonstrating by reference to the specification as a whole that 
the allegedly omitted element is, in fact, critical. Ex parte Nakano, 1998 WL 1735647 at *3. 
"The mere fact that other steps have been disclosed in the preferred embodiment does not 
render each and every step thereof an essential step." Id. The fact that an originally filed 
claim does not include the allegedly omitted elements "indicates the original intent to claim 
the system without either element." Ex parte Trummer, 1997 WL 33106615, at *4. 

B. Response to Examiner's Rejections 

The examiner rejected claims 12-19 and 21-23 "as being incomplete" for failing to 
include what he believes to be the essential step of "processing the positron annihilation data 

(e.g., see claim 12) and the prompt gamma ray data (e.g., see claim 13) " Office Action 

(1/16/07), p. 5. However, the examiner failed to demonstrate that the appellant described the 
allegedly omitted step as essential, 
a. Claim 21 

As an initial matter, the appellant notes that claim 21 contains the limitation 
"processing the positron annihilation data in accordance with a Doppler-broadening algorithm 
to produce output data indicative of an accumulation of lattice damage over time." Since 
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claim 21 already includes the precise limitation that the examiner is arguing needs to be 
included in the claims as a critical step, the examiner's rejection of claim 21 must be reversed, 
b. Claims 12-19, 22-23 
Turning to the remainder of the claims here rejected, there is no language in the 
specification which describes as "critical" or "essential" a step of processing data by Doppler 
broadening or positron lifetime algorithm. In addition, as described more fully below, the 
specification (in portions not cited by the examiner), including the abstract, describe 
embodiments of the invention that do not include a step for processing data by means of a 
Doppler broadening or positron lifetime algorithm. The appellant has set forth above 
examples from the specification, including the abstract, describing embodiments of the 
invention that do not include processing positron annihilation data or prompt gamma ray data 
See supra at pp. 29-30. In addition, the appellant notes that originally filed claims 12-19 and 
22-23 did not include a processing step. For example, independent claim 12, as originally 
filed, claimed: 

"A method for evaluating a material specimen, comprising: 

mounting a neutron source adjacent the material specimen; 

mounting a detector adjacent the material specimen; 

bombarding the material specimen with neutrons from the neutron 
source to create prompt gamma rays within the material specimen, some of 
the prompt gamma rays being emitted from the material specimen, some of 
the prompt gamma rays resulting in the formation of positrons within the 
material specimen by pair production; 

collecting positron annihilation data by detecting with the detector at 
least one emitted annihilation gamma ray resulting from the annihilation of a 
positron, the detector producing the positron annihilation data; 

storing the positron annihilation data on a data storage system for later 
retrieval and processing; and 
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continuing to collect and store positron annihilation data, the continued 
collected and stored positron annihilation data being indicative of an 
accumulation of lattice damage over time." 

No processing step is included. 

The Board in Ex parte Trummer, 1997 WL 33106615, at *4, found dispositive the fact 
that an originally filed claim did not include the allegedly omitted element. In that case, the 
examiner rejected one claim under §112, second paragraph, because that claim did not include 
controllers, interface module, blockable coupling element, control lines and other structure the 
examiner regarded as essential. The Board reversed. The Board's analysis applies equally 
well in this case: "The Examiner does not point to any evidence that Appellants regard their 
invention to be something other than what is claimed. The fact that the originally filed 
independent claim 1 did not include a controller or an interface module indicates an intent to 
claim the system without either element." Id. Similarly, here, the fact that the abstract, an 
alternative embodiment disclosed in the specification, and claims 12-19 and 22-23 do not 
include the allegedly omitted processing step demonstrates that the appellant intended to 
claim his invention without that step. The examiner therefore has failed to carry his burden of 
demonstrating that these claims are incomplete. The rejections of these claims under §112, 
second paragraph, must be reversed. 



CONCLUSION 

While the examiner argues that the claims here are deficient for being indefinite, non- 
enabled and lacking critical elements, to the contrary, it is the rejections here that are 
deficient. Thus, the examiner's non-enablement rejections were largely based on limitations 
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not found in the claims. Where the rejections were based on limitations actually found in the 
claims, the examiner failed to perform the enablement analyses required by In re Wands, 
supra. With respect to indefiniteness, the examiner failed to identify any specific claim 
language that he considered to be indefinite. Merely piggybacking indefiniteness rejections on 
non-enablement rejections is improper, as the Federal Circuit has held. None of these 
deficiencies were cured by the examiner's new omitted element rejections either. The 
examiner's arguments that certain claim limitations were "critical" and had to be included in 
the claims were belied by the specification, which disclosed embodiments that did not contain 
the allegedly omitted step and which did not describe the allegedly omitted limitation as 
"critical" or "essential." This case has been pending for over three years and has gone in and 
out of appeal. Being neither indefinite nor non-enabled, these claims are patentable and the 
appellant is entitled to have them issued. 



Respectfully submitted, 
FENNEMORE CRAIG, P.C. 




Susan E. Chetlin, PTO Reg. No. 59,722 
1700 Lincoln Street 
Suite 2625 
Denver, CO 80203 
Telephone: (303)291-3200 



Date: March 22, 2007 
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CLAIMS APPENDIX 

1 . A method for evaluating a material specimen, comprising: 
mounting a neutron source adjacent the material specimen; 
mounting a detector adjacent the material specimen; 

bombarding the material specimen with neutrons from the neutron source to create prompt 
gamma rays within the material specimen, some of the prompt gamma rays being emitted from 
the material specimen, some of the prompt gamma rays resulting in the formation of positrons 
within the material specimen by pair production; 

collecting positron annihilation data by detecting with the detector a plurality of emitted 
annihilation gamma rays resulting from the annihilation of positrons, the detector producing the 
positron annihilation data; 

processing collected positron annihilation data in accordance with a Doppler-broadening 
algorithm; and 

continuing to collect and process positron annihilation data to measure an accumulation of 
lattice damage over time. 

2. The method of claim 1, further comprising: 

collecting prompt gamma ray data by detecting with the detector a plurality of emitted 
prompt gamma rays, the detector producing the prompt gamma ray data; 

calculating positron lifetime data from the positron annihilation data and the prompt 
gamma ray data; and 
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continuing to collect positron annihilation data and prompt gamma ray data and calculate 
positron lifetime data to measure an accumulation of lattice damage over time. 

3. The method of claim 1, wherein said mounting a neutron source adjacent the material 
specimen comprises mounting the neutron source to the material specimen. 

4. The method of claim 3, wherein said mounting a detector adjacent the material specimen 
comprises mounting the detector to the material specimen. 

5. The method of claim 4, further comprising positioning a shield adjacent the neutron 
source to absorb stray neutrons. 

6. The method of claim 5, further comprising positioning a moderator between the neutron 
source and the material specimen. 

7. The method of claim 6, further comprising positioning a reflector adjacent the neutron 
source to reflect neutrons toward the material specimen. 

8. The method of claim 1 , wherein mounting a neutron source adjacent the material 
specimen comprises mounting an isotopic neutron source adjacent the material specimen. 
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9. The method of claim 8, wherein mounting an isotopic neutron source adjacent the 
material specimen comprises mounting a neutron source of 252 Cf. 

10. The method of claim 1, wherein continuing to collect and process positron annihilation 
data to measure an accumulation of lattice damage over time is performed while the material 
specimen is in service. 

11. (Canceled). 

12. A method for evaluating a material specimen, comprising: 
mounting a neutron source adjacent the material specimen; 
mounting a detector adjacent the material specimen; 

bombarding the material specimen with neutrons from the neutron source to create prompt 
gamma rays within the material specimen, some of the prompt gamma rays being emitted from 
the material specimen, some of the prompt gamma rays resulting in the formation of positrons 
within the material specimen by pair production; 

collecting positron annihilation data by detecting with the detector a plurality of emitted 
annihilation gamma rays resulting from the annihilation of positrons, the detector producing the 
positron annihilation data; 

storing the positron annihilation data on a data storage system for later retrieval and 
processing; and 

continuing to collect and store positron annihilation data, the continued collected and 
stored positron annihilation data being indicative of an accumulation of lattice damage over time. 
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13. The method of claim 12, further comprising: 

collecting prompt gamma ray data by detecting with the detector a plurality of emitted 
prompt gamma rays, the detector producing the prompt gamma ray data; 
storing prompt gamma ray data on the data storage system for later retrieval and processing; and 

continuing to collect and store prompt gamma ray data, the continued collected and stored 
prompt gamma ray data being indicative of an accumulation of lattice damage over time. 

14. The method of claim 12, wherein said mounting a neutron source adjacent the material 
specimen comprises mounting the neutron source to the material specimen. 

15. The method of claim 14, wherein said mounting a detector adjacent the material specimen 
comprises mounting the detector to the material specimen. 

1 6. The method of claim 1 5, further comprising positioning a shield adjacent the neutron 
source to absorb stray neutrons. 

1 7. The method of claim 1 6, further comprising positioning a moderator between the neutron 
source and the material specimen. 

1 8. The method of claim 12, wherein mounting a neutron source adjacent the material 
specimen comprises mounting an isotopic neutron source adjacent the material specimen. 
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19. The method of claim 12, wherein continuing to collect and store positron annihilation data 
is performed while the material specimen is in service. 



20. (Canceled). 

2 1 . The method of claim 1 2, further comprising: 
retrieving stored positron annihilation data; and 

processing the positron annihilation data in accordance with a Doppler-broadening 
algorithm to produce output data indicative of an accumulation of lattice damage over time. 

22. The method of claim 1, further comprising removing the neutron source before collecting 
positron annihilation data. 

23. The method of claim 12, further comprising removing the neutron source before 
collecting positron annihilation data. 

24. (Canceled). 



-46- 



EVIDENCE APPENDIX 
None. 



-47- 



RELATED PROCEEDINGS APPENDIX 
1 . Ex parte Akers, Appeal No. 2005-0855 (Bd. Pat. App. & Interf. 2005) (slip op.). 
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The opinion in support of the decision being entered today was not written 
for publication and is not binding precedent of the Board. 



UNITED STATES PATENT AND TRADEMARK OFFICE 



BEFORE THE BOARD OF PATENT APPEALS 
AND INTERFERENCES 



Ex parte DOUGLAS W. AKERS 



Appeal No. 2005-0855 
Application No. 10/269,807 



HEARD: May 18, 2005 



Before PATE, MCQUADE, and BAHR, Administrative Patent Judges . 
MCQUADE, Administrative Patent Judge . 

DECISION ON APPEAL 
Douglas W. Akers originally took this appeal from the final 
rejection (mailed June 16, 2003) of claims 1 and 8 through 39, 
all of the claims pending in the application. After the 
appellant's main and reply briefs (filed October 3, 2003 and 
January 13, 2004) and an examiner's answer (mailed December 1, 
2003) had been made of record, the examiner issued an Office 
action (mailed February 25, 2004) reopening prosecution and 
entering superseding rejections of all of the claims. Pursuant 
to 37 CFR § 1.193(b) (2) (ii), the appellant then submitted a 
request that the appeal be reinstated (filed April, 14, 2004) and 
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a supplemental brief (filed August 2, 2004). The examiner, 
implicitly granting the request, entered a second answer (mailed 
September 14, 2004) and forwarded the application to this Board 
for review of the new rejections of claims 1 and 8 through 39. 

THE INVENTION 

The invention relates to "methods . . . for performing non- 
destructive testing of materials using positron annihilation* 
'(specification, page 1). Representative claim 1 reads as 
follows: 

1. A method, comprising: 

determining whether a specimen to be tested includes at 
least one positron emitter therein that will be activated in 
response to photon bombardment; 

selecting a positron emitter to be activated; 

determining a threshold photon energy required to activate 
the selected positron emitter; 

determining a half-life of the selected positron emitter; 

and 

when the half-life of the selected positron emitter is less 
than a selected half-life, then performing a rapid 
activation/analysis process, said rapid activation/analysis 
process comprising: 

activating for an activation time the selected positron 
emitter by bombarding the specimen with photons having energies 
at least as great as the threshold photon energy; 

detecting for a detection time gamma rays produced by 
annihilation of positrons with electrons in the specimen; and 

repeating said steps of activating for an activation time 
and detecting for a detection time until detecting a sufficient 
number of gamma rays to determine at least one material 
characteristic of said specimen; 
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when the half-life of the selected positron emitter is 
greater than or equal to the selected half -life, then performing 
a normal activation/analysis process, said normal 
activation/analysis process comprising: 

activating the selected positron emitter by bombarding the 
specimen with photons having energies at least as great as the 
threshold photon energy; and 

detecting gamma rays produced by annihilation of positrons 
with electrons in the specimen. 

THE REJECTIONS 

Claims 1 and 8 through 39 stand rejected under 35 U.S.C. 
§ 112, first paragraph, as being based on a specification which 
fails to comply with the enablement requirement. 

Claims 1 and 8 through 39 also stand rejected under 35 
U.S.C. § 112, second paragraph, as failing to particularly point 
out and distinctly claim the subject matter the appellant regards 
as the invention. 

Attention is directed to the main, reply and supplemental 
briefs and the second answer for the respective positions of the 
appellant and the examiner regarding the merits of these 
rejections. 1 



1 In the Office action which reopened prosecution, claims 1 
and 8 through 10 additionally stood rejected under 35 U.S.C. 
§ 112, first paragraph, as being based on a specification which 
fails to comply with the written description requirement. Upon 
reconsideration, the examiner has withdrawn this rejection (see 
page 2 in the second answer) . Presumably, the examiner also has 
withdrawn any reasoning relating to the two remaining rejections 
which was set forth in the Office action but not restated in the 

(continued...) 
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DISCUSSION 

I. Preliminary matter 

In the supplemental brief (see pages 3, 4 and 17) and second 
answer (see pages 11 and 12), the appellant and the examiner 
debate the propriety of the Office action which reopened 
prosecution subsequent to appeal. As this matter is not directly 
connected with the merits of any issue involving a rejection of 
claims, it is reviewable by petition to the Director rather than 
by appeal to this Board (see In re Henqehold , 440 F.2d 1395, 
1403-04, 169 USPQ 473, 479 (CCPA 1971), and hence will not be 
further addressed in this decision. 

II. The 35 U.S. C. § 112, first paragraph, re-iection 

The dispositive issue with respect to the enablement 
requirement of 35 U.S.C. § 112, first paragraph, is whether the 
appellant's disclosure, considering the level of ordinary skill 
in the art as of the date of the application, would have enabled 
a person of such skill to make and use the claimed invention 
without undue experimentation. In re Strahilevitz . 668 F.2d 
1229, 1232, 212 USPQ 561, 563-64 (CCPA 1982). In calling into 
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second answer ( see Ex parte Emm , 118 USPQ 180, 181 (Bd. App. 
1957)). 
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question the enablement of the disclosure, the examiner has the 
initial burden of advancing acceptable reasoning inconsistent 
with enablement. Id. 

The examiner considers the appellant's disclosure to be non 
enabling in four respects: 

. . . There is neither an adequate description nor 
enabling disclosure as to what is encompassed by the 
term, "activating a positron emitter" . A "positron 
emitter" is inherently already activated, i.e., it is 
radioactive. 

At best, the use of the term "activating a 
positron emitter" is superfluous. 

At worst, the term would imply that a positron 
emitter is either being transformed to another positron 
emitter or its energy level is further raised by the 
energy of the activating photon. There is no support 
in the specification for either one of these two 
alternatives [second answer, page 3J; 

There is neither an adequate description nor 
enabling disclosure as to how and in what manner 
potential interferences in the data are accounted for 
in the analysis. For example, Claim 1 recites the step 
of bombarding the specimen with photons at least as 
great as the threshold photon energy required to 
activate the selected positron emitter and detecting 
gamma rays produced by annihilation of positrons with 
elections in the specimen. 

. . . [T]here is no support as to how one would 
differentiate between the signals from the selected 
positron emitter and from the non-selected ones [second 
answer , pages 6 and 7 ] ; 

There is neither an adequate description nor 
enabling disclosure as to how and in what manner one 
should select an algorithm from a plurality of 
available algorithms, modify/manipulate the selected 
algorithm and evaluate the constants in the selected 
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algorithm to fit Appellant's situation. Appellant 
himself admits in the specification that there is not 
only one known algorithm but rather a plurality of 
known positron lifetime algorithms (e.g. see paragraph 
0058), Doppler broadening algorithms (paragraph 0057) 
and three-dimensional imaging algorithms (e.g., see 
paragraph 0060) . 

. . . There is no support as to how and in what 
manner one selects the specific algorithm to use for 
the two activation/analysis processes, to evaluate 
requisite constants and to modify the selected 
algorithm to Appellant's situation [second answer, page 
7]; 

and 

There is neither an adequate description nor 
enabling disclosure as to how and in what manner one 
can determine when the half-life of the selected 
positron emitter is less than a selected half-life 
(sic), where the composition of the specimen is unknown 
and it is this composition that is desired to be 
determined. The same lack of support exists for the 
determination of when the half-life of the selected 
positron emitter is more than a selected half-life 
(sic) [second answer, page 8]. 

A fair assessment of the appellant's disclosure shows that 
the references therein to Activating 0 a positron emitter through 
photon bombardment of the specimen would be understood by the 
artisan as synonymous with "producing" or * forming" the positron 
emitter from a stable precursor. For example, the specification 
states that "photons 16 from the photon source 12 produce such 
neutron-deficient isotopes . , . referred to herein in the 



alternative as ^positron emitters'" (page 7, emphasis added), and 
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that "one way for generating positrons is through the formation 
within the specimen 18 of neutron-deficient isotopes, i.e., 
positron emitters" {page 20 , emphasis added) . In the same vein, 
the specification describes Tables I and II (on pages 21 and 22) 
as identifying n those isotopes that may be converted into 
positron emitters by photon bombardment" (page 20, emphasis 
added). Thus, even if the appellant's description of 
"activating" a positron emitter is somewhat inaccurate as urged 
by the examiner, it is so only in the most hyper- technical sense 
Read in context, the appellant's terminology would not pose any 
enablement problem from the perspective of a person having 
ordinary skill in the art. 

The examiner's concerns as to possible interference between 
respective signals produced by the selected positron emitter and 
non-selected positron emitters and as to the half-life of the 
selected positron emitter also are unfounded, primarily because 
they are not commensurate with the scope of the invention 
disclosed and claimed by the appellant. In this regard, the 
claims do not exclude the activation of multiple positron 
emitters, and indeed the appellant's specification (see pages 8, 
13 and 14) actually allows for such activation. The claims also 
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do not require the composition of the specimen to be unknown, and 
are not directed to a method for determining the composition of 
the specimen. 

Finally, while the appellant's disclosure does not convey 
much detail as to the various algorithms described therein, these 
algorithms admittedly are generally known in the art. The 
examiner has failed to advance any cogent reasoning as to why the 
disclosure would not have enabled a person having ordinary skill 
in the art to employ these algorithms, without undue 
experimentation, to the extent required by the appealed claims. 

Hence, the examiner's position that the appellant's 
disclosure, considering the level of ordinary skill in the art as 
of the date of the application, would not have enabled a person 
of such skill to make and use the claimed invention without undue 
experimentation is not well taken. Therefore, we shall not 
sustain the standing 35 U.S.C. § 112, first paragraph, rejection 
of 1 and 8 through 39. 

III. The 35 U.S.C. 5 112. second paragra ph, rejection 

The second paragraph of 35 U.S.C. § 112 requires claims to 
set out and circumscribe a particular area with a reasonable 
degree of precision and particularity. In re Johnson , 558 F.2d 
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1008, 1015, 194 USPQ 187, 193 (CCPA 1977). In determining 
whether this standard is met, the definiteness of the language 
employed in the claims must be analyzed, not in a vacuum, but 
always in light of the teachings of the prior art and of the 
particular application disclosure as it would be interpreted by 
one possessing the ordinary level of skill in the pertinent art. 
Id. 

The examiner views the appealed claims to be indefinite for 
reasons essentially similar to those listed above in connection 
with the enablement rejection. In this regard, the examiner 
points to the claim limitations relating to (1) the activation of 
a positron emitter, (2) the positron lifetime, Doppler 
broadening, three-dimensional imaging and activation/analysis 
process algorithms and (3) the determination of the half-life of 
the selected positron emitter (see pages 9 through 11 in the 
second answer) . For the reasons explained previously, the 
limitations pertaining to the activation of a positron emitter 
would have been readily understood by the artisan when read in 
light of the underlying specification. Furthermore, the 
examiner's criticisms of the claim limitations pertaining to the 
algorithms and the half-life of the selected positron emitter 



( ( 

Appeal No. 2005-0855 10 
Application No. 10/269,807 



focus on the breadth of these limitations. It is well settled, 
however, that mere breadth does not equate to indef initeness . 
In re Miller , 441 F.2d 689, 169 USPQ 597, 600 (CCPA 1971). 
Accordingly, the examiner's contention that the appealed claims 
are indefinite is unpersuasive. 

Thus, we shall not sustain the standing 35 U.S. C. § 112, 
second paragraph, rejection of claims 1 and 8 through 39. 

SUMMARY 

The decision of the examiner to reject claims 1 and 8 
through 39 is reversed. 
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The reliability of a least squares analysis of positron annihilation lifetime data represented as a sum of convoluted exponentials is 
investigated with computer-stimulated test data. The method of analysis, which is available as a computer program called 
" SPLMOD" [R.H. Vogel, Technical Report EMBL-DA08, European Molecular Biology Laboratory, Heidelberg, FRG, (1986)], 
avoids direct determination of the instrument resolution function by employing the decay curve of a reference material with a 
well-known single lifetime. The performance of the algorithm for extracting positron annihilation Lifetime information was evaluated 
by using several measures of the information content of the decay curve described recently by Schrader and Usmar [in: Positron 
Annihilation Studies of Fluids, ed. S. Sharma (World Scientific, Singapore, 1988)). The sensitivity of the algorithm to systematic 
errors including errors in the resolution function, shifts in the zero time channel and contamination of the reference decay curves by 
additional lifedme components was evaluated. The method provides excellent estimates of lifetime parameters when the reference 
decay Curve accurately reflects the resolution function of ihc sample data. However, the algorithm is extremely sensitive to errors in 
the instrument resolution function. 
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1. Introduction 

The analysis of positron annihilation lifetime data is 
an extremely difficult undertaking. The decay curve is 
usually represented by: 

C(r)= L«,A,exp(-X,0, (1) 

where a, is the number of positrons annihilating with 
natural lifetime, Analysts of sums of exponentials 
is notoriously difficult, particularly when the number of 
components, n, is unknown* In general, even for data 
with arbitrarily small noise levels, there exist a large 
number of solutions, [ot /s X /s nl, that fit the data to 
within the noise level. These difficulties are com- 
pounded by the fact that in a positron annihilation 
lifetime experiment one cannot observe the function 
C(f) directly, only its convolution with the instrument 
resolution function, R(i): 

y(t)~R(t) + C(t). (2) 

The difficulties associated with the analysis of eq. (2) 
are therefore two-fold: the choice of an appropriate 
representation of C(f) and the determination of Ihe 
function, R(t). 



The problems associated with analysis of data with 
eq. (1) are very well known. In general, we have no prior 
knowledge of the number of components, /i, in eq. (1). 
If n is too large the solutions become very unstable to 
noise in the data. Because of this difficulty, it is com- 
mon to assume only two or three components in eq. (1). 
This has the effect of stabilizing the solution to noise 
and in some simple systems may provide a physically 
meaningful representation of the decay. However, in 
more complex systems setting n too small will result in 
an inadequate representation of the system, resulting in 
a loss of valuable information and a set of a, and X, 
pairs that may have little or no physical meaning. 
Because of the non-linearity of C(r) in the X„ it is not 
possible to extend linear hypothesis tests to define an 
upper value for n '. However, Provcncher [1) has de- 
scribed a hypothesis test which incorporates a correc- 
tion for non-linearity that has been found to be very 
effective. It is employed in the programs DISCRETE, 
which solves sums of exponentials and SPLMOD which 
solves more general sums of one parameter functions, 
including exponentials and convoluted exponentials. The 
use of eq. (1) assumes that a discrete sum of exponen- 
tials is a physically meaningful representation of the 
data. For many simple, homogeneous systems this is 
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indeed the case, however, for complex, heterogeneous 
sysiems and systems in which distinct structural do- 
mains persist for limes that arc long relative to positron 
lifetimes, the lifetimes spectrum is expected to be more 
complex and may consist of a continuous distribution 
of lifetimes over some interval. For such sysiems an 
alternative approach is to replace the sum in cq. (1) by 
an integral [2-4]: 

C(t)-j~\a(\) exp(-Xf) dX = Aa(\)}, (3) 

\C(t) is simply the Laplace transform of \a{\), where 
a(X) is the annihilation rate probability density func- 
tion (pdf). The fraction of positrons annihilating with 
rates between X and X + dX is then given by a(X) dX 
(assuming <x(X) is appropriately normalized). In princi- 
ple, a(X) can be recovered by inverse Laplace transfor- 
mation of C(t). The problems of solution uniqueness 
and stability to noise associated with analysis of eq. (1) 
also plague the inversion of eq. (3). However, recent 
advances in the solution of Fredholm integral equations 
provide a number of reliable inversion methods [4-9]. 

Schrader [2,3] has pioneered the application of 
Laplace transform methods in the analysis of positron 
annihilation lifetime data and the application of infor- 
mation theory approaches to the determination of the 
information content of the experimental data. A par- 
ticularly important result of this work is the expression 
of the information content (Shannon Number) of the 
data in terms of spectrometer parameters (i.e. the num- 
ber of channels, the channel width, FWHM and the 
total number of annihilation events observed). Not only 
does this knowledge define the complexity of the func- 
tion, C(f), that is required to properly represent the 
data, but it also serves as a guide in the design Of 
experiments to maximize the information con Lent of the 
data. 

The accurate determination of the instrument resolu- 
tion function, R((), is a major problem in analysis of 
positron lifetime data. Probably the most widely used 
method is to determine R(t) from the prompt spectrum 
of a 60 Co source. However, the resolution function is 
very dependent on the distribution of energy deposited 
in the scintillators. As a result, the function determined 
with the 1.17 MeV and 1.33 MeV coincident 7- rays 
from 60 Co will not correctly represent the resolution 
function when data are collected with a 22 Na source. 
One approach, that attempts to overcome these difficul- 
ties, is to determine the resolution function from the 
lifetime data themselves. An analytical function, such as 
a sum of Gaussi&ns or a modified Gaussian, is assumed 
for R(t) and its parameters are treated as variables to 
be determined along with the lifetime parameters of 
interest [10,11]. However, ihere are several fundamental 
objections to this approach. The addition of further 
parameters to be determined by the analysis can severely 
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compound the problem of solution uniqueness. In ad- 
dition, the approach requires the choice of an analytical 
function to represent the resolution properties of the 
instrument. An incorrect choice will, of course, transmit 
errors to the lifetime estimates. 

The resolution function may also be determined in a 
separate experiment from the decay curve of a standard 
sample with a known single lifetime, such as a well-an- 
nealed metal [12]. In principle, this method can yield 
very good estimates of the instrument resolution func- 
tion. However, in practice, armihilation in the source 
and positron trapping at the surface can give rise to 
additional components. An alternative approach is to 
employ a 207 Bi source. 207 Bi emits two y-rays separated 
by a time interval of 186 ps, the energies of which are 
almost identical to those employed in a positron life- 
time experiment with a 22 Na source. Although the 
method does not introduce unwanted source compo- 
nents, it is not yet clear if the energies of the Y^rays 
emitted by 207 Bj are sufficiently similar to those occur- 
ring with **Na that the types of error associated with 
the use of 60 Co are entirely avoided. 

At present, the use of standard samples with 22 Na 
sources appears to offer the best hope for reliable 
determination of the instrument resolution function. We 
have been interested in the use of a deconvolution 
method due to Gauduchon and Wahl [13) and 
Wijnaendts van Resandt et al. [14] for this purpose. A 
preUminary report of its use in the analysis of positron 
armihilation lifetime data with numerical Laplace inver- 
sion methods has appeared [4] but no tests of its perfor- 
mance have been presented. Our purpose in this paper 
is to explore the reliability of this deconvolution method 
in the analysis of positron annihilation lifetime data. 
The efficiency of the algorithm for extracting lifetime 
information is evaluated by using several measures of 
the information content of the decay curve described by 
Schrader and Usmar [3]. In this work, C(t) is repre- 
sented by a sum of exponential terms as in eq. (1) and 
the analysis is performed with the FORTRAN program, 
SPLMOD, developed by Vogel [15]. Although eq, (1) is 
more restrictive than eq. (3), it can be a valid represen- 
tation of annihilation lifetimes in simple, homogeneous 
systems and is considerably simpler for our purpose 
because the measure of information content of the data, 
the Shannon Number, is directly related to the number 
of components in the lifetime spectrum. 



2. Deconvolution of the instrument resolution function 

Deconvolution of the instrument resolution function 
follows the method of Gauduchon and Wahl [13] and 
Wijnaendts van Resandt et aL [14], which avoids the 
direct determination of R(i) by measuring, y r (t) % the s 
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decoy curve of a reference material with a well-known 
single lifetime. \,: r 

y r (t) = *(/)• <;(/), (4) 

where 

Ct(')*=« r \ exp(-A r r). ( 5 ) 
Convolving eq. (2) with C r (r) we obtain: 
C"r<Ov<0-*<0*C(r). ( 6 ) 
Here wc assume that MO and >/(,) have the same 
resolution function, R{t). This would be true if the 
reference and sample decay curves could be measured 
simultaneously, but this is not possible in a positron 
annihilauon lifetime experiment. Any instrument drift 
occurring between measurement of ihe sample and ref. 
erence decay curves will therefore introduce error into 
the lifetime estimates. This, of course, is true of all 
methods thai involve a separate measurement of the 
resolution function of the instrument and can only be 
avoided by improving the stability of the instrument. 
Eq. (6) can be rewritten as: 

a f \ r {ex P (-\ r r))* JK ( f )= >f ( / )^(r) 1 (7) 
which on taking Laplace transforms gives: 

)^HP)=MP)C(P), (8) 

where the tilde denotes Laplace transformation and p 
is the Laplace variable. Eq. (8) is solved for y(p): 

?(p)~«7 1 K\KMp)C(p)+pMp)c(p)}. 

(*) 

which on inverting Laplace transforms gives: 

*^(0+-Kr(0-[^(0 + C(0)5(r)]}. 

(10) 

If C(f) is represented by eq. (1), eq. (10) becomes: 

HO 

<") 

If X f is known and a r is absorbed into the a it eq. (11) 
must be solved for the 2n parameters and' A/. If X 
is not known, there arc In + 1 parameters. . 
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The current version of SPLrylOD (version 3) distributed 
by Vogel [15] can perform the following types of analy- 
sis of relevance to Positron Armihilation Lifetime (PAL) 
spectroscopy: 

(a) Solution of eq. (12) when the resolution function is 
known; i e. with /(X„ t k ) *(/,). exp(- A,.,,) and 
S\K*k) 88 1 (to account for the background). 

(b) Solution of eq. (12) using a reference material with a 
known single arinihilation lifetime, X r . i. e . with 
f(K ^)-<*,-MM'*)*«p(-*,r and 
8x('k)r 1. 

'(c) Solution of eq. (12) using a reference material with a 
single annihilation lifetime which need not be 
known; i.e. with /(X,, t k ) = H P (r*)*exp(- A,r*), 

ftO*)-l and gi(t k )**y r (t k ). 

We have employed option (b) in our simulations be- 
cause it provides more stable solutions than option (c). 

SPLMOD performs a weighted least squares analysis 
in which /(X„ t k ) is approximated by a sum of normal- 
ized cubic B-splines. The program employs Gauss- 
Newton muumization, with an extensive search of 
parameter space to determine the initial parameter val- 
ues. The convolution, ^ P (^)*exp(-X / r) t is calculated 
[16] USin8 a reCmiVC formula of Grinvald and Steinberg 

/(A,. '* + i)»{/(X„ t k )+Q5Aty T (t k )) 

Xcxp(^X^r) +0.5A/y r (f^ J ), (13) 

where At is the channel width of the multi-channel 
analyzer in which the data are collected. SPLMOD 
performs several analyses in which n in eq. (1) takes the 

cdt U wU1; 2 ' 3% \7 "-w where n mAx is set by the user. 
SPLMOD supplies a number of statistical measures 
including Fisher probabilities and random runs prob- 
abilities, to aid in the choice of an appropriate solution. 
Measures of the information content of the data may 
also be employed to define an appropriate solution [3 J. 

Note that it Is not necessary to independently define 
the zero time point of the spectrometer. This mforroa- 
tion is already contained in the reference decay curve. 
However, there is an option in SPLMOD to correct for 
small changes in the zero time point due to instrument 
drift occurring between measurement of sample and 
reference decay curves. 



»™t° f ^ ° f convoluted ^PonentiaLs with 4. Simulation of decay curves 



SPLMOD' 

SPLMOD analyzes data that can be represented as a 
sum of one parameter functions: 



(12) 



Positron annihilation decay curves have been simu- 
lated by convoluting eq. (1) with a Gaussian resolution 
function with FWHM equal to 2<n/inT- 



A(0--4=- exp -(;/<,)' 



(14) 
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This convolution is given by: 
y(t) = 0.5 Z exp{ -X,/ + (\,o/2) 2 } 



xerfc(A,cr/2-r/o). (15) 
where erfc is the compliment of the error function: 
e r fc(x) - l - cri(x) = 1-1/' exp(-r 2 ) dr. (16) 

Conversion of the continuous /unction. y(t) % to a table 
ol channel counts, y k , was performed by numerical 
integration, A background and appropriate Poisson 
noise were added to each channel. The reference data 
were simulated in the same way except that no back- 
ground was added. The reference and sample spectra 
were calculated over the same range of channels in each 



5. Information content of the simulated data 

Schradcr and Usmar [3] have given a detailed discus- 
sion of the determination of the information content 
(Shannon Number) of decay curves and have provided 
an approximate expression for the Shannon Number in 
terms of instrument setting of the spectrometer. The 
Shannon number, in the context of PAL spectroscopy, 
is a measure of the number of resolvable lifetimes that 
can be geometrically fitted between the largest and 
smallest resolvable lifetimes. This measure is an ap- 
proximation because it depends on ones definition of 
the relative scatter in the decay curve and on the values 
thai are chosen for the largest and smallest resolvable 
lifetimes. Schrader and Usmar [3] define the Shannon 
Number, F t as; 

f - \n{a^At/a 2 FWHM), (17) 

where a 5 and a 2 are constants, FWHM is the full 
width at half the peak height of the resolution function, 
N ch is the number of channels to the right of the 
aero-lime channel, jV 0 , employed in the analysis (a 
maximum value for is defined by the number of 
channels for which the count is greater than the square- 
root of the background count) and w m4A is defined by: 

tW^TT" 1 ln(v^27/n), (18) 

tj is a measure of the effective relative noise in the decay 
curve and is defined as: 



k-N u 



(19) 



w mftX has a special significance in the eigenfuncUon 
expansion of Laplace transforms which we need not be 
concerned with here [3,5). Eq. (18) is derived by assum- 
ing that the shortest resolvable lifetime is equal to 



a 2 F^HM and the longest resolvable lifetime is equal to 
a^^At. A resolution factor, p, is defined as the ratio 
of resolvable lifeiimcs and it follows that F is the 
number of resolvable lifetimes that can be geometrically 
fitted between a 2 FWHM and a^N^Au TTie constants 
£i and a 2 need to be determined. In practice, the 
shortest resolvable lifetime is probably less than FWHM 
and the longest resolvable lifetime is longer than N^At. 
so a 2 < 1 and a, > 1. Schrader and Usmar [3] suggest 
a l «• 3 and a 2 =» j as reasonable values when a t FWHM 
> 34 r, but they suggest replacing a 2 FWHM by 34/ if 
a i FWHM<34f. 

The resolution factor is given by: 

/> = exp(V«nu*) (20) 

and is therefore dependent on the definition of the 
effective relative noise of the decay curve (eq. (19)). 
Estimates of a lt a 2 and the correctness of cq. (19) are 
discussed below. 



6. Analysis of simulated data 

Table 1 lists the results of an analysis of simulated 
PAL data together with the information parameters for 
the decay curve and will serve as a point of reference for 
our subsequent discussions. The simulated data contain 
three components. The reference has a single lifetime of 
125 ps. A total of 800 channels were analyzed, starting 
30 channels (0.75 ns) to the left of the peak. A ploi of 



Tahiti 1 

Calculated information content of simulated data and results 
of its analysis with SPLMOD 



FWHM 290 ps 


P 


2.6 


At 25 ps 


V 


0.015 


Reference lifetime 125 ps 




3.25 


Total events I0 6 


maximum 
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°\ 02 f 


Resolvable lifetimes (ns) 




min. 


max. 


1 1 3.8 


0.29 


11 


2 0.5 5.2 


0.15 


22 


3 0.33 6.0 


0.1 


33 




Lifetime* (nS> 




Simulated 0.123 


0.400 


2.000 


Fitted 0.124 


0,397 


1.997 




Alpha 




Simulated 0130 


0.470 


0.400 


Fitted 0.125 


0.474 


0.4O1 




Background 




Simulated 


100 




Fitted 


100.5 
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Fig. 1. Example of the fil to the simulated data obtained wilh SPLMOD. (a) Simulated decay curve calculated with the lifetime 
parameters given in table 1. (b) Weighted residuals about ihe zero lime channel for ihe fil given in table I, 



the simulated data and the residuals is shown in fig. 1. 
The fit to the data is clearly excellent with no evidence 
of any oscillations of the residuals about the peak. The 
smallest ratio of simulated lifetimes is 3.2, while the 
calculated resolution factor is 2.6 SPLMOD has little 
difficulty in resolving the three components, and pro- 
vides good estimates of the parameters for the short-lived 
component (125 ps). The information content of the 
data, measured by ihe Shannon number, is clearly suffi- 
cient to define all the components even with a x and a 2 
assumed equal to one. 

These results already provide some clues as to the 
appropriate value for a 2 . It must be £ 0.33 because the 
short-lived component is well defined. Table 2 lists the 
results obtained when FWHM is increased with 4f — 
12.5 ps. The range of channels employed for analysis 
was chosen such that the first and last channels were 



Table 2 

The effect of FWHM on the resolution of lifetime components. 
Simulation parameters are the same as given in tabic 1 except 
that At » 12.5 ps 







Simulated 
data 


FWHM 
500 ps 


600 ps 


700 ps 


Lifetime 


1 


0.125 


0.121 - 


0.116 


0.122 




2 


0.400 


0.400 


0.398 


0.326 




3 


2.000 


2.000 


1.998 


1-633 


Alpha 


1 


0.130 


0.129 


0.125 


0.069 




2 


0.470 


0.471 


0.475 


0.526 




3 


0.400 


0.399 


0.400 


0.406 



sufficient to sample ihe baseline on both sides of the 

peak. Good estimates of the parameters are obtained 
even with FWHM values as large as 500 to 600 ps 
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(cr 2 - 0.25 to 0.20). Beyond this value, SPLMOD is still 
able lo identify three components but the parameter 
values become very poor. Increasing the total count to 
10 1 does not improve these estimates. At a FWHM of 
290 ps, a value that is typical of modem PAL spec- 
trometers employing plastic scintillators, the shortest 
lifetime that cart be resolved by SPLMOD is approxi- 
mately 50 ps (i.e. a 2 = 0.17). An estimate of a, was 



determined by varying the number of channels em- 
ployed in the analysis. In all cases the start channel was 
sufficiently far to the left of the peak (-0.74 ns) to 
sample the baseline. With ^ = 1940, the longest life- 
time that can be resolved is approximately 100 ns 

We have examined the ability of the program to 
resolve lifetime components and were particularly inter- 



Table 3 

The resolution of lifetime components with SPLMOD- Variations in the effective relative noise and resolution factor were generated 
by changing the total number of events. Simulation parameters were; FWHM «=» 290 p$, At — 12.5, reference lifetime =125 ps and 
background -100 





Simulated 


Fitted 


Fitted 


Fitted 




Total events 




10 6 


10* 


io- 




Relative noise 




0.024 


0.054 


0.074 




w m*< 




3.0 


2.4 


2.4 




Resolution factor 




2.9 


3.6 


4.1 




F(a, = 4, o 2 -0.25) 




6.1 


4.5 


2.7 




Maximum N ch 




959 


451 


64 




Lifetime 1 


0.125 


0.122 


0.113 


0.124 




(OS) 2 


0.400 


0.386 


0.313 


0.609 




3 


1.280 


1.213 


0.895 


2.786 




4 


4.096 


3.947 


3.338 






Alpha 1 


0.250 


0.245 


0.223 


0.341 




2 


0.250 


0.247 


0.210 


0.345 




3 


0.250 


0.248 


0.258 


0.314 


j 


4 


0.250 


0.260 


0.310 







Table 4 

The effect of systematic errors on lifetime parameters estimated by SPLMOD. Simulation parameters were FWHM = 290 ps, 
Ai «= 12.5 ps» background « 100, total events =» 10 7 and reference lifetime -125 ps 





Simulated 


Incorrect reference FWHM 


Zero channel 




data 


285 ps 


295 ps 


shift 


Lifetime 1 


0.125 


0.127 


0.138 


0.134 


(ns) 2 


0.400 


0.409 


0.410 


0.409 


3 


2.000 


2.004 


2.004 


2.005 


4 




0.003 


100.0 


0.002 


Alpha 1 


0.130 


0.663 


0.148 


0.247 


2 


0.470 


0.171 


0.449 


0.402 


3 


0.400 


0.150 


0.395 


0.351 


4 




0.016 


0.008 


0.001 


Background 


100 


100.4 


91.1 


99.9 




Simulated 


Contamination of reference by additional 






data 


component with Hf etlrac - 450 ps 




Contaminant % 




5 


10 


15 


Lifetime 1 


0.125 


0.122 


0.120 


0.118 


(ns) 2 


0.400 


0.382 


0.365 


0.347 


3 


2.000 


2.001 


2.003 


2.005 


Alpha 1 


0.130 


0.131 


0.135 


0.137 


2 


0.470 


0.476 


0.474 


0.A77 


3 


0.400 


0.393 


0.391 


0.386 


Background 


100 


100.3 


100.2 


100.1 
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csted in determining if the value of the resolution factor 
calculated according to eq. (20) correctly predicted the 
point at which SPLMOD failed to resolve lifetime com- 
ponents. An example of the results for such a test is 
shown in table 3. The simulated data consisted of four 
components with successive lifetimes increasing by a 
factor of 3.2. When the total number of events is 10*, 
the calculated resolution factor is 2.9 and the Shannon 
number, F % is 6.1 SPLMOD is able to resolve the four 
components and provides reasonably good estimates of 
the lifetime parameters. When the total number of 
events is reduced to 10 V the effective relative noise is 
doubled, the resolution factor increases to 3.6 and the 
Shannon number decreases to 4.5. SPLMOD is still able 
to resolve four components but the lifetime parameter 
values become poor. When the total number of events is 
reduced to 10 4 , the relative noise increases to 0.074, the 



resolution factor increases to 4.1 and the Shannon 
Number decreases to 2.7. SPLMOD is unable to resolve 
the four components. Eq. (19) thus appears to provide a 
good estimate of the effective relative noise in the data 
as measured by the predictive power of the calculated 
resolution factor. 



7. Systematic errors 

We have examined the sensitivity of the SPLMOD 
algorithm to various systematic errors in the data. The 
most important sources of systematic error arise from: 
(i) variations in the resolution function for the reference 
and sample data, (H) differences in the zero time point 
for the reference and sample, and (iii) contamination of 
the reference with additional lifetime components. 
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The fitted lifetime parameters are extremely sensitive 
to variations in the FWHM of the resolution function. 
The lifetime parameters obtained from an analysis of a 
three component sample (FWHM = 290 ps) with refer- 
ence decays having incorrect FWHM values (285 and 
295 ps) arc given in table 4. The lifetime parameters 
determined by SPLMOD when the FWHM for the 
reference is less ihan that employed for the sample are 
very poor. The intensities of all components are grossly 
in error. In addition, a spurious shori-Hved component 
appears. When the FWHM of the reference is made 
larger than thai employed for the sample, a spurious 
long-lived component appears (100 ns). Estimates of the 
other lifetime parameters are in error, but arc not nearly 
as poor as those determined with the smaller reference 
FWHM. Plots of the residuals are shown in fig. 2. The 
fits to the data are poor about the peak. 

A shift in the zero channel of the reference with 
respect to the sample data also generates poor estimates 
of the lifetime parameters and an additional spurious 
short-lived component. The example given in table 4 
and fig. 3 was obtained by shifting the reference decay 
curve by -2.5 ps with respect to the sample decay 
curve. This shift is a fifth .of the channel width em- 
ployed. SPLMOD can correct for small shifts in the 
zero time channel by assuming that v(r + c) has been 
used instead of y(t) and then expanding eq. (11) in a 
Taylor series about y(t) to the first order term. This 
gives an additional linear term in eq. (12). However, in 
this case the correction did little to improve the fit to 
the data. 

The analysis assumes the reference to have a single 
lifetime component, an ideal which is unlikely to be 
exact for real systems. Annihilation in the source and 
trapping at the surface of the reference material un- 
doubtedly introduce additional lifetime component*. 



The effect of introducing an additional component with 
a lifetime of 450 ps into the reference decay curve is 
shown in table 4. The estimated lifetime parameters are 
relatively insensitive to the presence of the contaminant, 
even when it represents 15% of the reference decay. The 
largest error appears in the sample component which 
has a lifetime closest to the contaminant lifetime. 



8. Conclusions 

The results presented in the previous sections suggest 
that the deconvolution and least squares analysis meth- 
ods employed by SPLMOD can provide excellent esti- 
mates Of positron annihilation lifetime parameters for 
systems that can be represented by a discrete sum of 
exponentials. As we have mentioned, this representation 
is too restrictive for many systems of interest in which 
distributions of positron lifetimes are expected. For 
such systems numerical Laplace inversion methods are 
obviously to be preferred. The results of our tests of a 
modified version of the program, CONT1N [3,5,6], which 
solves integral equations with convoluted kernels will be 
published shortly. It is nevertheless clear from the pre- 
sent work, that estimates of the lifetime parameters are 
critically dependent on accurate knowledge of the in- 
strument resolution function. We doubt that the decon- 
volution method employed here is any more sensitive to 
errors in the resolution function than other algorithms 
that have been described in the literature or that arc 
currently in use. For example, lifetimes reported using 
60 Co prompt spectra to measure the instrument resolu- 
tion function are undoubtedly in error since it is known 
that this method provides a poor estimate of the true 
resolution function occurring with 22 Na sources [3]. 
There is little doubt that these problems will become 
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even more severe as we attempt to determine continu- 
ous lifetime distributions with numerical Laplace inver- 
sion techniques. 
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Figure 14. A Comparison of the Peak Region of the Raw Data and the 
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Figure 16. A Plot of Rj Values for the Raw Data vs. quench 
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ABSTRACT 



Positron Annihilation Spectroscopy (PAS) has become a useful tool for 
the nondestructive investigation of the defect structure of materials. Of the 
existing PAS measurement techniques Ooppler broadening appears to offer 
several practical advantages. These advantages include low radioactive source 
strength requirements, comparatively short experimental times, and minimal 
geometrical restrictions. Unfortunately, the technique suffers from some 
Inherent disadvantages. These Include comparatively poor instrumental energy 
resolution and a limited ability to distinguish between multiple defect types. 

In an attempt to overcome the difficulties associated with the use of 
Ooppler broadened PAS the Fast Fourier Transform/Power Spectrum (FFT/PS) de- 
convolution algorithm has been developed in this work. This algorithm is com- 
putationally efficient, noMteratl ve, and requires minimal a priori, informa- 
tion. The utility of the FFT/PS deconvolution algorithm is demonstrated for 
high purity aluminum samples containing known concentrations of vacancies or 
prismatic dislocation loops. It is shown that the FFT/PS procedure 1s capable 
of automatically correcting for accidental and/or Intentional changes in the 
instrument resolution function. In addition, based on the limited data anal- 
yzed in this work, It appears that both the sensitivity and the accuracy of 
the energy distribution Hneshape parameters are improved by deconvolution. 
Finally, it is demonstrated that if the filtering error introduced by the 
FFT/PS procedure is controlled, then 1t is possible to distinguish between 
samples containing predominately vacancy type defects and those with prismatic 
dislocation loops as the dominant defect type. 
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CHA1TER I 
INTRODUCTION 



Over 50 years ago. Dirac postulated the existence of an atomic particle 
with a mass equal to that of an electron, but which is positively charged. In 
1932 these positively charged electrons, referred to as positrons, were first 
observed by C. D. Anderson as a by-product of cosmic ray showers. Since that 
time positron research has been pursued In such varied disciplines as high en- 
ergy nuclear physics, astrophysics, and materials science. It is within the 
broad context of materials science, and more specifically in the area of non- 
destructive testing (NDT) of aluminum, that the interaction of positrons with 
matter will be examined in this work. Before going on to describe the use of 
positrons as an NDT tool in the investigation of the defect structure of 
metals, a brief review of positron physics will be presented. A comprehensive 
treatment of this subject has been written by West (1974). 

As postulated by Dirac 1n 1935. positrons and electrons are complemen- 
tary states of matter such that when the two particles encounter one another 
they will annihilate each other. Their kinetic and rest mass energies will be 
dissipated as the energy of one, two. or three gamma rays In such a way as to 
conserve momentum. It can be shown that the ratio of the two-photon to three- 
photon annihilation event is 372:1 (Stevens 1980). while the ratio of the two- 
photon to one-photon events is ~2.6xlo6 : i, (West 1974). With these ratios as 
justification the remainder of this work will consider only the annihilation 
event yielding two photons. If both the electron and positron were at rest, a 
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positrons in vacancy type lattice defects. This effect was first observed by 
MacKenzie et al_. and Dekhtyar et aU 1n 1964. Schultz et al_. (1981) describe 
the changes in the annihilation characteristics due to the localization of a 
positron In a monovacancy type defect as follows: 1) Because the localized 
positron has a smaller overlap of its wave function with that of the high ener- 
gy core electrons than 1t does when diffusing through the nearly perfect lat- 
tice, the annihilation rate with high energy core electrons decreases. 2) The 
density of conduction electrons in the vacancy is also decreased resulting in a 
decreased annihilation rate with conduction electrons. However, the percentage 
reduction is less than in the case of core electrons so that the net result is 
an inc^ase in the fraction of annihilation events involving low energy conduc- 
tion electrons, thus resulting in a decrease In the width of the Doppler broad- 
ened annihilation radiation distribution. 3) The potential for electrons in 
the vacancy is higher than normal resulting in the conduction electrons travel- 
ling more slowly. This in turn also decreases the observed Doppler broadening. 
4) A second order effect is that the overlap of the positron wave function with 
weakly bound outer-core electrons is reduced less than is the overlap with 
inner-core electrons. This too leads to less broadening of the observed 511 
KeV line. S) The uncertainty principle requires that localization of the posi- 
tron be accompanied by an increased spread in its momentum distribution (and a 
corresponding decrease in its de Broglie wavelength). As a result, the contri- 
bution of the positron to the center-of-mass velocity increases corresponding 
to a broadening of the annihilation energy spectrum (Schultz et a].. 1981). The 
first four effects lead to a narrowing of the annihilation energy spectrum 
while the fifth effect tends to broaden the spectrum. The combination of these 
effects results in more sharply peaked (narrower) Doppler spectra for metals 



As discussed by Siegel (1980) It is important to note that the posi- 
tron is not equally sensitive to all lattice defects. It is sensitive to 
defect regions which can be characterized as having a lower-than-average elec- 
tron density, for example, vacancies, small vacancy clusters, voids, and dis- 
location lines or jogs. On the other hand, positrons are relatively insensi- 
tive to defect regions with a hlgher-than-average electron density such as 
interstitials. Indeed this defect specificity 1s one of the advantages PAS. 
measurements offer as compared with other NDT techniques such as residual 
electron resistivity (Siegel 1980). 

In addition, positrons are capable of probing the entire depth of the 
valence band (Mackenzie et a, . 1964, which is in contrast with techniques such 
as electrical and thermal conductivity, thermoelectric power. Hall coef- 
ficient, magnetic susceptibility, and Knight shift, all of which yield infor- 
mation restricted to the electrons near the Fermi surface. 

PAS measurements are an attractive NDT tool since the positron causes 
negligible perturbation of the electronic structure as it passes through a 
crystal (Johnson et a], 1978) and the positron-electron Interaction does not . 
appreciably alter the momentum distribution of the electrons (Szotek and 
Chmielewskl 1981). A further advantage of the technique is that it yields 
information regarding both the electron momentum distribution as well as the 
electron density within the sample (Siegel 1980). 

The three most common PAS techniques are lifetime measurements, angu- 
lar correlation, and Ooppler broadening. The effective electron density in 
the region from which the positron annihilates can most easily be determined 
using lifetime measurement techniques. These measurements usually use "Na as 
a source material since it simultaneously « lo psec) emits a y-ray along 
: "A th thepOSUron - This MeV Y -ray acts to start the, timing apparatus ... 



which is stopped by the arrival of the ~ 511 KeV y-ray signifying the occur- 
rence of the annihilation event. The positron lifetime depends on the con- 
centration and distribution of the electrons in the region of the material 
through which the thermallzed positron drifts before its annihilation. This 
"mean" lifetime will of course be affected by any process, such as the 
trapping of positrons at defects, which alters the local concentration of 
electrons seen by the positron (Coleman 1977). These mean lifetimes, which 
are on the order of hundreds of picoseconds, can be measured to an accuracy of 
a few picoseconds with modern equipment. 

As mentioned previously, when the annihilation event does occur the 
electron-positron pair will usually have a small net momentum in the labora- 
tory frame. This net momentum affects both the energy of the resulting y-ray 
as well as Its direction of emission as shown In Figure 1. The relationships 
between the momentum of the annihilating y-ray, the measured Doppler shift in 
energy, and the measured angle 5+ are described in Appendix B. 

1 

\ Center of Mass 
\ Velocity V 




Figure 1. A Schematic Diagram Showing the Effects of the Electron-Positron 
Pair Momentum on the Energy and Olrection of the Two Resulting 
Annihilation y-Rays in the Laboratory Reference Frame. 
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Angular correlation measurements, which measure the deviation from co- 
linearity of the two photons emitted from the annihilation process, were first 
reported by DeBenedetti and his coworkers in I960. The so called long-slit 
geometry (DeBenedetti 1950) is commonly chosen since it provides both good 
angular resolution and adequate counting rates. However, other geometries can 
be employed for specific applications. The angular distribution of a sample 
is determined by measuring the coincidence counting rate as a function of 
angular deviation from colinearlty. In well annealed metals this data gener- 
ally takes the form of a central Inverted parabola with Gaussian wings. The 
instrument resolution is good, typically ~ 0.5 mrad, but the counting rate is 
comparatively slow and the required source strengths are relatively high. 

Alternatively, Doppler broadening measurements could be used to deter- 
mine the net momentum distribution of the annihilating electron -positron pairs 
by measuring the energies of the resulting Y-rays. The center of mass momen- 
tum of the annihilating pair along the line of emission of the f-ray gives 
rise to Doppler shifts of the Y-ray energies. The distribution of these 
momentum components then results in a Doppler broadening to give the 511 KeV 
line a width of approximately 511 (V F /2c) KeV where V F (~ 108 cra/s ) is tne Fermi 
velocity of the conduction electrons and c 1s the speed of light (Hughes 
1982). Details of the experimental procedure for collecting Doppler-broadened 
spectra are given in Chapter II. 

Although lifetime measurements, angular correlation, and Doppler- 
broadening measurements are by far the most common PAS techniques there are 
two others which should be mentioned. One of the drawbacks of using conven- 
tional sources such as 22 Ua or 68 Ge is that only certain depths of the mater- 
ial under investigation are sampled extensively by these positrons. Hughes et 
ii«..0»82). give th e equati on: 
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* ■ 4^ («•) ■* 255.5 («e)ev ()-2) 
where «E is the Ooppler Energy resolution 1n eV 
and SB is the Angular Correlator. Resolution in mrad. 

The best Ooppler systems ha»e a resolution of about 1.2 KeV which 
corresponds to an ang„,ar correlation resolution of 5 „ rao . w „ } 
correlation system h„e a resolution approximately an order of magnitude 
better than this. For Ooppler broadening technics to effectively compete 
with angular correlation techniques in bign precision applications it is 
necessary to improve the effective resolution of Dopp.er PAS. Scbultt et ah 
Omi point out that theoretically a ,0-f.ld increase in effective resolving 
power of Ooppler measurements can be achieved by deconvolution. 

Oeconvolution is the mathematical technique used to separate the 
instrumental broadening effects from the true Ooppler broadened positron 
annihilation ,-ray spectrum. Instrumenta, broadening is a completed problem 
caused by many different factors. For example, physical properties of the 
detector crystal such as frontal dead layers, gentry of the active region 
and regions of poor charge collection, can all >ead to line broadening. The 
primary contribution to Instrumenta, broadening ,„ . s „,id state semiconductor 
detector is due to statistica, Muctuations in the division of absorbed energy 
between ionization and heating of the crysta, .attlce. i.e. the Fann. effect. 
(McNelles and Campbell 1975). 

The shape of the instrument resolution function, as well as the mag- 
nitude of the line broadening, is most often determined from the 5,4 KeV ,-ray 
emitted from the isotope »5 Sr . The d1ffjajUy ^ ^ ^ ^ ^ ^ 

uoppier-broadening of the 5„ KeV way Hne in the absence of defects in a 
metallic crystal is ofthe same order- (-KeV) as the detector resolution, and 
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changes Induced by defects are significantly smaller. Fortunately, the defect 
induced changes in the deconvolved way energy distribution are much more 
dramatic than in the instrumental ly smeared raw data (McNelles and Campbell 
197b). In addition, removing instrumental effects offers another advantage. 
Such a procedure allows direct comparison of experimental results obtained in 
different laboratories. If results are reported without deconvolutlon they 
must be nondimensionallzed in some way or else the Instrument resolution func- 
tion must also be reported. In either case the interpretation of results is 
complicated in the absence of a deconvolution procedure. 

In light of the above discussion it was decided that the results pre- 
sented in this work would be for deconvoluted Ooppler spectra. Unfortunately, 
the mathematics involved are somewhat complicated. Because of this, most of 
the existing deconvolution schemes are iterative in nature and hence computa- 
tionally inefficient. One of the major goals of this work was to develop a 
deconvolution scheme for Doppler broadened positron annihilation measurements 
which is model Independent, does not require iteration, and which yields 
reproducible results. 

In addition to the need for deconvolution, DBPAS also currently suf- 
fers from an Inability to distinguish various defect types. For example, 
vacancies and dislocations tend to alter the annihilation lineshape in quali- 
tatively the same manner. That is, in both cases the spectrum becomes more 
narrow as the defect concentration increases. Most of the lineshape param- 
eters currently used to characterize Doppler broadened spectra are unable to 
distinguish the effects of different defect types because they describe only 
gross changes in lineshape. For example, Mackenzie's S-parameter, which is 
the most commonly used lineshape parameter, essentially measures the width of 
ri?L^ r ^^i^^li ^IjZOl^Hence, for any Qlven spactr- ^ 
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an analysis yields a single value. This is useful if there 1s only one degree 
of freedom in the series of experiments. If two variables are allowed to 
change, however, it is impossible to separate their effects. Clearly, 
additional parameters are required. The other parameters selected must be 
uncorrected with the S-parameter. That is, 1f the second parameter chosen Is 
a function of the first parameter only, then no new information Is gained. 
Campbell (1977) and Schultz et aU (1981) have discussed the use of other 
possible parameters to describe llneshapes. Although these parameters may be 
more sensitive to variations in specific defect concentrations they are 
correlated with S, and therefore of little help in differentiating defect 
types. 

In light of the above discussion it can be seen that 1t would be de- 
sirable to find a family of lineshape parameters which is capable of disting- 
uishing various types of metallurgical defects. For example, one of the 
groups of parameters to be examined in this work is the central moments of the 
T-ray energy distribution. It 1s well known that the values of these moments 
give Information regarding the shape of the distribution. In fact, it can be 
shown that if all of the moments are known one can construct the distribution. 
The even moments give information about the width of the distribution while 
the odd moments vanish for symmetrical distributions such as the ones to be 
examined in this work (the symmetry question Is addressed in Chapter III). 
The second central moment Is called the variance and its square root is the 
familiar standard deviation. The fourth central moment is called the Kurtosis 
or flatness. The advantage of moment analysis is that the information ob- 
tained is not redundant. For example, it is very possible to have two dis- 
tributions with equal standard deviations (second moments) and yet widely 
y£5 C£!!L/^l ~~?~ --=: ~}L of 1 inesha P e parameters suggests the 
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possibility of separating the effects of multiple metallurgical defect types 
using DBPAS. 

The above brief review of positron annihilation spectroscopy as an NOT 
tool for the study of metallurgical defects indicates that the technique with 
the greatest potential for practical use outside the laboratory is Uoppler- 
broadenlng. This technique offers the advantages of low source strength re- 
quirements, the ability to use 6«Ge as a source material 1n order to reduce 
surface effects, and few geometrical experimental constraints. Doppler broad- 
ening currently suffers from relatively poor energy resolution and an Inabil- 
ity to distinguish different metallurgical defect types. Based on these ob- 
servations, the purpose of this work is twofold: 1) to develop an efficient, 
model independent, noniterative deconvolution method, and 2) to develop a 
family of Hneshape parameters capable of distinguishing various defect 
types. 



CHAPTER II 
EXPERIMENTAL PROCEDURES 



This chapter 1s divided into two sections. The first section des- 
cribes the Ooppler broadened PAS instrumentation and the methodology employed 
to collect Doppler spectra. The second section outlines the thermomechanical 
treatments used to create the desired defect concentrations in the metallic 
samples which were studied in this work. 

The Popper PAS System 

The heart of the Doppler broadened PAS instrumentation system is the 
single crystal high purity Germanium (HPGe) detector. This detector is mount- 
ed at the end of a cold finger in a cryostat and, to reduce thermal noise, is 
operated at liquid nitrogen temperatures. The cylindrical detector crystal 
used in this work is 46.7 m in diameter and 46.2 mm In length. When the 
annihilation Y-rays impinge upon the crystal their energy is absorbed by two 
mechanisms. The first is the creation of electron-hole pairs within the 
crystal. These electron-hole pairs are subsequently swept apart by the 
applied detector bias (4000 V) and flow through an external circuit creating a 
current pulse whose integral with respect to time has an amplitude propor- 
tional to the incident y-ray energy. The second absorption mechanism is heat- 
ing of the crystal. The ratio of these two absorption mechanisms is known as 



the Fanno factor and the statistical fluctuations in the division of absorbed 
energy between ionization and heating is the primary contribution to Instru- 
mental line broadening in solid state serai conductor detectors. For further 
details of the absorption process see Moore (1983). 

The resulting current pulse from the detector Is rather weak and un- 
suitable for cable transmission. Because of this a preamplifier (Ortec Model 
120-5) is mounted directly adjacent to the detector output. The detector and 
preamp are O.C. coupled, as is the entire system, and powered from the same 
A.C. source to prevent ground loops. The preamp is a charge sensitive device 
capable of accepting signals corresponding to the Y-ray energy range of 0-10 
MeV. Its output 1s a voltage pulse with an amplitude proportional to the 
integrated charge supplied by the detector and hence proportional to the inci- 
dent Y-ray energy. The magnitude of the voltage pulse amplitude 1s related to 
the incident Y-ray energy by the approximate conversion factor of 280 mV per 
MeV. Therefore for a typical 511 KeV Y-ray the output of the preamp would be 
approximately 143 mV. In addition to converting the detector charge pulse to 
a voltage pulse with an intensity suitable for cable transmission the preamp 
also shapes the resulting voltage pulse so that it is acceptable to the linear 
amplifier. 

The linear amplifier used in this work is an Ortec Model 571. It 
accepts the shaped pulses from the preamp and amplifies them to suitable volt- 
ages for transmission to the multichannel analyzer. The output pulses are 
active filter shaped (skewed Gaussian approximations) and have amplitudes In 
the 0-10 V range. Unipolar outputs are used to optimize the signal-to-no1se 
ratio (resolution). Since the main amplifier Is operated within its linear 
range the amplitude of the output voltage pulse is proportional to the inci- 
dent Y-ray en ergy as measured b y the HPGe detector. 
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The output from the linear amplifier is routed to an Ortec Model 408A 
biased amplifier. This device accepts the output pulses from the main ampli- 
fier and allows one to selectively view a:.y region of the initial 0-10 MeV en- 
ergy range. In this work the biased amplifier was used to select the energy 
window of 483 to 548 KeV. The output pulses from the biased amplifier serve 
as the Input to the analog-to-d1g1tal (ADC) converter located within the Nor- 
land IT-53UO multi-channel analyzer (MCA). In this device the amplitude of 
the voltage pulse 1s converted to a digital representation of the energy of 
the detected T-ray and subsequently dynamically displayed on the MCA's CRT in 
the form of a histogram of energy versus count intensity. The MCA can be pro- 
grammed to accept input signals until a predetermined number of counts 1n the 
peak channel is reached. At this time the experiment is terminated and a com- 
plete Uoppler broadened spectrum 1s displayed on the CRT. The energy distri- 
bution is then ready for transmission to the computer for subsequent data man- 
ipulation. 

The two most Important parameters which characterize the quality of a 
Doppler PAS system are the width, in units of energy, of the Instrument reso- 
lution function and the energy resolution, in units of eV per channel, of the 
MCA. The width of the instrument resolution function (IRF) is, in essence, a 
measure of the quality of the detector crystal, although 1t Is also affected 
by the associated amplification electronics as well as the ADC. The IRF is 
most often determined by examining the response of the system to the monoener- 
getic (E = 614 KeV) photons emitted from the Isotope «5 Sr# If the instrumen- 
tation was perfect the energy distribution, as recorded by the MCA, would con- 
sist of a single bar in the histogram in the channel corresponding to 514 KeV. 
However, in any real instrumentation this line will be broadened for reasons 
^ 1 ^:^^'~ J - amo ^-^ ^WgtaTjHw^n^J s, mo st., of ten-re^.., 
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ported in terms of the energy width of the 85 Sr line as miilire<1 acr£)$s t||f 
two channels which contain half the number of counts in the peak channel. 
This energy measure is commonly named the full width half maximum (FWHM). 
Figure 2 shows typical IRF for the Doppler PAS system used In this work. As 
can be seen in this figure the FWHM of the IRF Is about 1.28 KeV. 

The second parameter which describes the quality of the Doppler PAS 
instrumentation is the energy resolution of the MCA which is determined by the 
following procedure. First, a Doppler broadened spectrum is collected and the 
channel corresponding to its centroid is determined. The centroid channel is 
assigned the energy value 511 KeV. Next, an energy spectrum Is obtained for 
the Isotope a5 Sr. The channel corresponding to the centroid of this distribu- 
tion is calculated and assigned the energy value 514 KeV. The MCA energy 
resolution is then easily calculated by dividing the energy difference (3 KeV) 
by the number of channels between the centroids. In this work the 511 KeV 
channel was number 439 while channel number 486 corresponded to 514 KeV. 
Hence, the MCA energy resolution Is 3000 eV/(486 -439) or " 64 eV per channel. 
Clearly this value can be adjusted using the biased amplifier. The reason for 
selecting the 64 eV/channel Is to obtain the best possible resolution while 
maintaining the ability to display the entire region of interest (483 to 548 
KeV) on the 1024 channel MCA. The reasons for requiring the above energy 
window will be discussed in Chapter III. 

While the MCA energy resolution is a parameter which, to a certain ex- 
tent, can be varied, the IRF is essentially predetermined by the instrumenta- 
tion. The one factor contributing to the width of the IRF which can be con- 
trolled is the data acquisition, or count rate. It is well known that in- 
creasing the count rate decreases the instrument resolution. Therefore, one 



21 



40 



m 



30 



X 

to 



o 
o 

o 

a> 

jQ 
3 



20 



■1.28 keV 



X 



10- 



• 



• 



• 



• 



512 



~?i7 

keV 



p mmmm 



516 



Fiqure 2. A Typical Instrument Resolution Function for the Dopoler Pas 
^stemjjsed_in this Work, Nate that the Full-Width- Half - 



22 

keeping the time required to collect a spectrum reasonably short. Moore 
(1983) demonstrated that with the Doppler instrumentation used in this work 
very reproducible results can be obtained with count rates of < 2 kHz. The 
count rate 1s easily adjusted by changing the sample to detector distance. 
All of the results reported in this work were obtained using count rates of 
< 2 kHz with an associated experiment time of ~ 30 minutes resulting in a peak 
height of 40K counts and ' 3 x 10* total integrated counts. 

The positron source isotope in all experiments was ^Ge which was 
chosen to avoid surface effects. The mean depth of penetration in aluminum 
for such positrons is ~ 0.50 nr.. The activity of the source was initially 50 
uC. The half-life of 68 6 e is ~ 275 days and hence the activity of the source 
did not change significantly for any single set of experiments since a typical 
set of measurements requires less than a week to accumulate. It should be 
pointed out that the activity of the source does not measurably affect the 
observed annihilation spectra, provided the count rate remains constant. 

Although the count rate can be easily adjusted by changing the source 
to sample distance such a change introduces an additional error source. This 
error source 1s related to the source/sample geometry. It has been demon- 
strated that the observed spectra may depend on the angle with which the gamma 
rays impinge upon the detector crystal (Sroedskjaer 1983). It has also been 
shown that in high precision experiments among samples of the same material, 
one should use samples of nearly identical thickness to reduce the error asso- 
ciated with small angle scattering (Smedskjaer 1983). 

The "Ge source was deposited on a thin (~ 0.05 m) nickel foil for 
handling convenience. This results In what Is known as a source contribution 
to the annihilation lineshape. That is, in each experiment some of the posl- 
■ ..-■■.■■. ^ 
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events will contribute to the overall annihilation spectrum. However, the 
ratio of source annihilation events to sample annihilation events 1s small (< 
1%), and the source contribution is neglected in the remainder of this work. 

In all measurements the source and sample were arranged in what 1s 
commonly called the sandwich configuration. In this configuration the posi- 
tron source is placed between two identically treated metallic samples. The 
samples are thick enough, typically a few millimeters, to guarantee that 
virtually all of the positrons annihilate within the material being examined. 



Specimen Preparation 

The material examined in this work is high purity, large grain size 
aluminum purchased from the Materials Research Corporation. The aluminum was 
zone refined to a nominal purity of 99.9995% before being vacuum melted and 
cold fabricated to size. The chemical analysis, in PPM of metallic inter- 
stitial, is listed in Appendix C. The aluminum was supplied as a 13 mm 
diameter rod from which ~ 1.5 mm thick disks were machined. Each of these 
disks had an ~ 1.5 mm hole drilled through it near its edge. The purpose of 
this hole was to allow the disks to be suspended from a titanium wire while 
within the vertical tube furnace. After the disks were machined to size, they 
were given a full anneal heat treatment consisting of 2 hours at 400°C fol- 
lowed by furnace cooling overnight. The as annealed aluminum disks were the 
starting material for all subsequent thermal treatments. 

The thermal treatments used in this work are very similar to those 
outlined by Cotterill etaK (1972). In this work the authors describe a two 
step time/temperature cycle which allows one to produce, In a controlled 
f - a A h l 0 /!jL./ir. s ^j m ^.g s yaca ncies as the domi na|g -dgfgilL typjg4_.and_t.hej>, 



samples with dislocation loops as the major defect type. The actual thermal 
treatment employed in this work is as follows: 

1) Heat the samples to the desired quenching temperature and hold 
them at temperature for one hour to establish thermal equilibrium. 

2) Quench the samples from the furnace into - 293°K water. 

3) Dry the samples within 5 to li seconds. 

4) Quench from - 293°K to a liquid nitrogen bath at 77°K. This 
procedure has been shown to produce primarily vacancies, the concentration of 
which is a direct function of the quenching temperature. 

&) After DBPAS characterization, the aluminum disks are then annealed 
at 3S3°K for thirty minutes, which results in the production of prismatic 
dislocation loops by vacancy coalescence. 

6) The samples are cooled to 77°K and the PAS measurements corres- 
ponding to the dislocation loops are made. 

7) Finally the disks are annealed at S23°K for two hours to remove 
the dislocations and return the material to the annealed state. 

Several of the steps in this procedure require further explanation. 
The quenching procedure is taken from Cotterill et a]_. (1972) and their 
results together with the following discussion forms the basis for the 
assumption used in this work that the PAS measurements made after step 4 are 
for a fairly uniform distribution of vacancies of concentration approximately 
equal to the equilibrium concentration at the quenching temperature. 

Quenching experiments are usually performed on thin metallic foils so 
that the quench rate remains constant across the sample. However, as dis- 
cussed previously, constraints imposed by the positron annihilation experi- 
ments put a lower limit on the sample thickness of ~ 1.5 m . For a complete 
review of the complications which can evolve durino <m»nrhinn . 
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sample thickness see, for example, Cotterill et al_. (1965). The main problem: 
to be addressed here are the formation of vacancy clusters and the loss of 
vacancies at the sample surfaces. The second problem was investigated in a 
paper by Wilkes and Harrison (1972) in which they observed the quenching of 3 
mm thick disks of aluminum and reported a very low vacancy loss. Since the 
samples used in this work are only ~ 1.5 mm thick and since the quench pro- 
cedure is similar to that employed by Wilkes and Harrison, 1t will be assumed 
that the vacancy loss during the quench can be ignored. As for the formation 
of vacancy clusters, the positron lifetime measurements of Cotterill et al. 
(1972) demonstrate that for this quenching procedure the probability of form- 
ing bivacancies, or larger clusters, is small. 

It is very important to dry the samples thoroughly between the water 
quench and the liquid nitrogen quench. This is necessary to avoid the forma- 
tion of a thin layer of ice on the specimen surface which would contribute an 
undesirable component to the measured PAS spectrum. As pointed out by 
Cotterill et aj_. (1972) the annealing of the vacancies during this short dry- 
ing procedure may be neglected since the mean lifetime of vacancies (at room 
temperature) in a specimen quenched from 873°K is - 100s and increases for 
lower quenching temperatures. It 1s, of course, desirable to keep the drying 
time as short as possible. 

Step 5 of the procedure, the formation of prismatic dislocation loops 
from the collapse of disk shaped vacancy clusters during a low temperature 
anneal of aluminum, was reported by Hirsch et an in 1958. Cotterill et al. 
(1972) demonstrate that the appropriate thermal treatment consisted of a 30 
minute anneal at 353°K. The theoretical relationship between the vacancy 
concentration after the quench, C v , and the concentration of dislocation 
cores, C d (gi ven as an atomic fraction), is presented in Appendix 6. 
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The apparatus used to perform the thermal treatments is shown in 
Figure 3. A quart; tube passes through a vertical tube furnace and is 
attached to the room temperature water reservoir by a -quick disconnect" 
coupling. Just al)ove the water )eve , .„ ^ reservofr h arson ^ 

port. Ail thenna, treatments were performed in an argon atmosphere to avoid 
oxidation. The other e« of the tube is connected to , flow meter which is 
used to monitor the argon flow rate. The other side of the flow meter is 
vented to atmosphere. The temperature is «asured by a chromel-alumel ther»- 
couple positioned inside the tube furnace and adjacent to the aluminum disks. 
These disks are supported by , titanium wire, the other end of which is con- 
nected to an iron slug. The slug, and hence the specimens, are held in posi- 
tion in the furnace by a magnet located outside of the tube. The actual 
quench is performed simply by pu.Mng the magnet away from the tube and allow- 
ing the entire specimen-wire-slug arrangement to fa., by gravity into the 
quench oath. It is. of course, important for the systea, to be designed so 
that the specimens fall freely through the tube and quick.y reach the quench 
bath. If the specimen-wire-slug arrangement does not fall freely the te^era- 
ture of tne samples when they enter the water will be significantly lower than 
the temperature measured by the thermocouple. 
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CHAPTER III 
DECONVOLUTION OF OOPPLER PAS MEASUREMENTS 



As discussed in the Introduction, compared to angular correlation (AC) 
measurements, the major disadvantage of the Doppler PAS technique is Its in- 
herently poor energy resolution. In order for Doppler methods to be competi- 
tive with AC methods 1n high resolution applications the Instrumental line 
broadening must be removed from the experimentally obtained Ooppler spectra. 
Mathematically, the technique which accomplishes this separation is deconvolu- 
tion. In this chapter the deconvolution of Ooppler PAS measurements will be 
discussed. 

The chapter 1s divided into four sections. In the first section the 
mathematics of deconvolution are presented. This section justifies the use of 
deconvolution techniques to separate Instrumental effects from the desired 
true Ooppler-broadened energy spectrum. In addition, the inherent advantages 
and mathematical difficulties associated with the use of Fourier analysis to 
solve the problem are addressed. 

In the second section of this chapter the methods introduced by other 
Investigators to solve the incorrectly posed deconvolution problem, along with 
their inherent advantages and disadvantages, are discussed. The third section 
presents the FFT/power spectrum deconvolution method (Schaf fer et a].. 1984); 
and describes its use in detail. In the final section of this chapter the 
FFT/power spectrum method is demonstrated using a neutron irradiated high 
purity single crystal sapphire sample. 
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The Mathematics of Deconvolution 

The first objective of this section is to demonstrate the validity of 
using the convolution equation to describe the relationship between h(x), the 
experimentally observed energy distribution, f(x). the true Doppler broadened 
eneryy distribution, and r(x), the instrument resolution function (IRF). In 
Figure 4(a) the instrument resolution function r(x) is shown. Recall that this 
function represents the instrumental broadening of what should be a single 
energy line at the centroid of r(x). The function r(x) is assumed to be norm- 
alized to unit area. In order to understand the relationship between these 
three functions, f(x) should be viewed as a histogram comprised of segments of 
width ax. If it can be assumed that the instrumental response is linear (the 
detector and all associated electronics are operated within their linear range) 
then each segment of the f(x) histogram will be broadened by the unit area IRF 
r(x). Hence, the value of the observed function h(x) at x* will have contribu- 
tions from many of the f(x) histogram segments in the neighborhood of x'. Fig- 
ure 4(b) is a schematic representation of the contribution to h(x') from the 
histogram section of f(x) at a distance xj from the origin. This contribution 
Ah^x') will be the product of the area f(x])Ax and r(x'- Xl ), that is: 

Ahi(x') » f(x])r(x' - Xl )Ax (3.1) 
The total value of h(x') will be the sum of the contributions from all pos- 
sible histogram segments in f(x). Thus, 

Mx') = I Ah^x') = I f( Xl )r(x' - Xi )Ax (3.2) 
i =- 00 i =-«• 

or as the width of the segments approaches zero the summation in Eq. 3.2 is 
replaced by an integral to yield, 

GD 

h(x') = / f(x)r(x' - x)dx (3.3) 
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Clearly the same procedure holds for any value of x In h(x) and the final form 
of the relationship, with a change of variables, Is: 

00 

Mx) = f(y)r(x - y)dy (3#4) 

Eq. 3.4 is the convolution equation which serves as the starting point for the 
remainder of the mathematics used in this work. 

The advantage of using Fourier analysis when dealing with convolution 
problems can be seen clearly when Eq. 3.4 is rewritten as: 

H(w) = F(w)-R(w) (3>5) 
where the functions in Eq. 3.5 are the Fourier transforms of the respective 
functions in Eq. 3.4 according to the general relationship: 

G(W) Vi? L e " iWX9(X,d * (3.6) 
Note that in Fourier space convolution is simply a complex multiplication. 
The equivalence of Eqs. 3.4 and 3.5 is demonstrated in Appendix D. 
Rearranging Eq. 3.5 yields 

F(w) = H(w)/R(w) {3>7) 
from which one can find f(x) using the Inverse Fourier transform 

f(x)^/>xH^ dw ^ 

Eq. 3.8 is an expression for the true Doppler broadened energy spectrum in 
terms of the Fourier transforms of the experimentally observed function h(x) 
and the instrument resolution function r(x). However, the uniqueness of the 
solution given in Eq. 3.8 requires that h(x) be known exactly. This condition 
1s never satisfied in a real experiment and this fact significantly compli- 
cates the mathematics of the problem. The observed function h*(x) will be 
related to the "true" function h(x), as defined in En. a. a f„ 
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function f(x) by: 

max |h*(x)-h(x)|<c (3.9) 
-°°<y<+» 

where e is a measure of the noise associated with the stochastic nature of the 
experiment. 

It can then be shown (see Appendix E) that there are two solutions 
fj(x) and f 2 (x) = fj(x) + C s1n(2*nx/(b-a) J with corresponding function hj and 
h 2 (as defined by Eq. 3.4) both of which satisfy the inequality (3.9). Hence, 
both f] and f 2 must be considered as solutions to Eq. 3.4 with a left hand 
side h*(x) taken from an experiment to within an error e (Turchin et al. 
1971). In fact, there is an infinite set of such solutions a subset of which 
contain indefinitely large rapidly oscillating components. 

In order to understand the mathematical source of these difficulties 
Eq. 3.8 was investigated further. If in Eq. 3.8 H(w) 1s replaced by H*(w) 

with - . 

H*( W ) =-L- / e lw V(x)dx (3.10) 

and h*(x) = h(x) + 6( x ) (3.11) 

where 6 is a stationary random process with zero mean, then it has been shown 
(Kautian 19S8) that the dispersion of f(x) is 

0(f) = <f 2 > - <f> 2 = _L f G (»? ? dw (3.12) 

fli - |R(w)r 

where 6(w) is the power spectrum of 6(x) 

and |R(w)| ♦ 0 for |w| ♦ « 
Therefore the dispersion of the solution f(x) can remain finite only if the 
power spectrum of the noise approaches zero rapidly as w approaches infinity. 
In practice this will not be the case since the noise will inevitably contain a 
"white noise" component such that G(w) remains finite as |w| ♦ « and the dis- 



33 

pension of the solution f(x) therefore becomes infinite (Turchln et aK 1971). 

Alternatively one can understand the source of the mathematical diffi- 
culties as follows. The observed function h*(x) contains high frequency noise 
components which are not physically realizable in the true function h(x). In 
Eq. 3.8 these high frequency components are divided by small eigenvalues. R(w), 
which results in larye oscillations in the solution f(x) (Turchin et al_. 1971). 

Recognizing these difficulties, previous investigators have proposed 
several deconvolution techniques which appear to successfully solve the prob- 
lem. Since the accuracy of the deconvoluted spectrum can never really be de- 
termined, except in certain carefully constructed numerical experiments, the 
criticisms of a specific method tend to relate to the amount of computer time 
required, the complexity of the programming effort, the degree to which the 
analysis may be automated, or to the need for some a priori information con- 
cerning the function f(x). As a preface to the introduction of the FFT/power 
spectrum deconvolution technique, the following section is a review of the 
prior techniques that have been developed to deconvolute Doppler PAS data. 

Previous Deconvolution Techniques 

One of the early Doppler PAS deconvolution schemes Introduced by Hotz 
et al_. (1968) is basically a matrix approach. The convolution in Its discrete 
form is 

h(x) - r(y - x)f(y)Ay (3J3) 

where Ay is the incremental width (resolution) 
and the other functions have their usual meanings 
This can be rewritten in matrix form as 
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b i " A j1 c j (3.14) 
where b is the discrete experimental vector (spectrum) 

c is the unknown discrete vector before instrument broadening 
and A is a rectangular matrix in which each row is the detector re- 
sponse function appropriate to the corresponding component of c. 
in the method developed by Hotz et al_. (1968) the matrix A is determined using 
« 5 Sr and the vector b is the output of the Doppler system MCA for a given 
experiment. In this technique an initial estimate of the unknown vector c is 
made, say cO), and a corresponding vector bO) is calculated from Eq. 3.14. 
The quality of C 0) is determined by finding X 2 for the fit between b and b<D 
as defined by 

X2 » ftb.d) - b 1 )/b i (3 . lb) 

Next cO) is modified in a manner described in the original work to create a 
second estimate, c<2). of the desired function c. Then bW is found from Eq. 
3.14 and a new X* value determined. The process Is repeated until x2 is mini- 
mized. This process is extremely involved and inefficient. In addition the 
initial estimate C 0> as well as the modification criteria used to generate 
subsequent cO) estimates can influence the final result. 

In more recent years the matrix approach has been examined in greater 
detail (Aulick and GalUe 1981). Although this work i S not related to Doppler 
PAS measurements, the mathematical problem is the same. It Is demonstrated 
that 3.14 should be rewritten as: 

b = * (3.16) 

where b = b + e 

and e is the measurement error 
- The least squares solut ion of Eg. 3.16 is then 
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Atb " c (3.17) 
where A* is the pseudo-inverse (Penrose 1956) of A. 

However, due to the fact that b contains noise Eq. 3.16 is incorrectly posed 
and A will be ill-conditioned. Hence, the elements of A t are large and the 
error in b will be magnified, obscuring the true solution which is 

C = Atb (3.18) 
Aulick and Gallie (1981) point out that to avoid this problem one needs to 
calculate solutions to Eq. 3.16 in a way which will damp out the error in b. 
One way of accomplishing this goal is to supply some additional a priori in- 
formation about the shape of the desired vector c. The additional Information 
recasts the problem as a constrained optimization problem. This technique is 
called regularization. 

An example of the type of a priori information commonly supplied 
appears in the deconvolution technique used by Reddy and Carrigan (1970). In 
their approach the resolution function is assumed to be the sum of two 
Gaussians with the same centroid but differing in their amplitudes and widths. 
The shape of f(x) (returning to the more common notation) is assumed to be a 
truncated parabola based on the free electron model (Reddy and Carrigan 1970). 

Introducing these constraints on the shapes of r(x) and f(x) into the 
least square solution problem yields a physically plausible result. That is. 
the function f(x) which is obtained by this method does not contain indefl- ' 
nitely large and rapidly oscillating components typical of the results obtain- 
ed by evaluating Eq. 3.8 directly. The difficulties associated with this 
method are related to the validity of the a ^rion assumptions. For example, 
the resolution function r(x) is constrained to be symmetric. In addition as- 
suming the shape of the true Ooppler energy distribution. f(x). to be a simple 

Other, inves^igatoxs.^. : _„, 
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(Mackenzie etal_. 1976) suggest the modeling of f(x) using the superposition 
of Gaussian and inverted parabolic functions as representing contributions 
from core and valence electrons respectively. While it Is not the purpose of 
this work to discuss the validity of the various electron distribution models, 
the very fact that there exists more than one analytic prediction of the shape 
of f(x) is strong reason to develop a deconvolution technique which places 
minimal constraints on the shape of f(x). 

One of the most widely used deconvolution techniques in Doppler PAS 
studies is an Iteratively based method Introduced by Jackman et al.. (1974). 
In this method an initial estimate of the desired result f(x), say f,(x). Is 
selected. The estimate f,(x) is then convoluted with r(x) to obtain h^x), 
the first estimate of the experimentally obtained energy distribution h(x). 
The function h^x) is then compared point -by-point with h(x) and if satisfac- 
tory agreement is obtained f,(x) is assumed to be a reasonable approximation 
to the true solution. If unsatisfactory agreement is obtained then fj(x) is 
modified to produce f 2 (x) and the process is repeated. 

The initial estimate of f(x), the modification procedure for obtaining 
f 2 (x) from f,(x), and the criteria for stopping the iteration procedure are 
all current topics of investigation In the literature. In the original method 
the experimentally obtained function h(x) served as the initial estimate 
f](x). The modification algorithm was 

f 2(*i) " fl(xi) x ( I h(x,)/ I hWx,)) (3.19) 

j=i-2 J j=1-2 J 

Subsequent work demonstrated that the convergence rate of the iteration pro- 
cess could be improved if the summations in Eq. 3.19 were taken over nine 
channels instead of the original five and if the summation ratio was squared 
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cedure is stopped when the comparison of h(x) and h N (x) yields a x2 value less 
than some predetermined minimum (Jackman et aU 1974) or when the sum of the 
squares of the deviations falls below some predetermined level (Dannefaer and 
Kerr 197b). 

The main advantage of these iteratively based techniques is that a so- 
lution is obtained without any restriction on the shape of f(x). While this 
is indeed an attractive feature of the method it also suffers from some inher- 
ent disadvantages. These disadvantages include the amount of computer time 
required, the possible dependence of the result on the initial estimate of 
f(x). and the possibility of introducing large oscillations into the solution 
if the stopping criterion 1s selected with too small an error bound. This 
last difficulty Is of course related to the mathematical nature of the incor- 
rectly posed problem. 

The Iterative deconvolution technique essentially Ignores the statis- 
tical nature of the problem. An alternate approach is to treat Eq. 3.4 as a 
problem of mathematical statistics and introduce a regularizer (a priori in- 
formation) in the form of some probability distribution for f(x) (Turchin et 
aU 1971). This is the method of statistical regularizatlon which offers the 
advantage of error analysis. The result of this process is the most probable 
f(x) for a given h(x) measured with mean square errors S(x) which are indepen- 
dent of x and normally distributed with zero mean (Goldanskii et a]_. 1978). 
The regularizer 1s a form of a priori information concerning the shape of f(x) 
and hence somewhat restrictive. However, such an approach appears to yield 
reasonable results if a solution can be obtained using minimal a priori infor- 
mation concerning the smoothness of f(x). 

In 1981 the Fourier cut-integral method was introduced by Szotek and 
e l^ k : W I«- i^i^^onl^^-^^&^cXttenr-as^ -. 



38 

^'-P 0 e ' 'm** (3.20) 



f(x) .-LA e (1wx) 



The difference in the equations is the introduction of finite limits of inte- 
gration. If p 0 is then chosen small enough the solutions obtained are physi- 
cally plausible. The mathematical explanation for this result is that the in- 
sertion of finite limits of Integration forces the dispersion of f(x) to 
remain finite. (Refer to Eq. 3.12 and the associated discussion). The prob- 
lem associated with this method is that the optimum value of p 0 is not readily 
predictable. Approximations for f(x) are obtained using different values for 
p 0 . and one must then select the "best" approximation to the true f(x). This 
is typically done by convoluting the function with r(x) and then comparing the 
result with h(x). As p 0 becomes too large the solution begins to exhibit 
large oscillations as would be expected from the mathematical theory. 

In 1982 Verxerk published an extensive mathematical treatment of the 
deconvolution problem. In this work he discusses the relative strengths and 
weaknesses of many of the currently employed deconvolution techniques. Ver- 
kerk then proposes a method which in many ways is very similar to the Fourier 
cut-integral method. The crux of this method is that the high frequency 
Fourier coefficients of h(x) do not contain significant information and conse- 
quently only a limited number (say K) of Fourier coefficients H(w) contain the 
essential information from measurement of h(x). It should be pointed out that 
if H(w) contains N Fourier coefficients then deleting the N-K high frequency 
coefficients in the reconstruction of f(x) is equivalent to inserting finite 
limits of integration in Eq. 3.4 as in the Fourier cut-integral method. 

The deletion of the high frequency coefficients is equivalent to a 
rectangular low-pass filter and it is well known that the application of such 
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a filter to F(w) will cause nonphysical oscillations in the approximation to 
f(x). 

In summary, there are two different types of errors associated with 
the estimate of f(x): 1) statistical errors, which are amplified by the de- 
convolution operation and 2) systematic errors due to filtering, which are 
largest in the peak region and may cause nonphysical oscillations 1n the tails 
of the deconvoluted signal. Error type 1 Increases as the number of Fourier 
coefficients used in reconstruction increases, while error type 2 decreases 
(Verkerk, 1982). Thus, by using such an approach the problem 1s reduced to 
selecting the optimum number of coefficients. One disadvantage of Verkerk's 
method, as with the Fourier cut-integral method, 1s that the procedure must be 
repeated several times for different numbers of coefficients until the optimum 
number is determined. It would be desirable to develop a simple one-step 
method to determine the appropriate number of coefficients to use in recon- 
structing f(x) from F(w). and indeed this 1s one of the primary advantages of 
the FFT/power spectrum deconvolution method which 1s described in the follow- 
ing section. 

The FFT/Pow er Spectrum Deconvolution Method 

The deconvolution method developed 1n this work is a fast Fourier 
transform (FFT) based method which is simple to execute, fast, easily auto- 
mated, and assumes little or no a priori information concerning the shape of 
f(x). The distinctive feature of this method is that the number of Fourier 
coefficients retained in the reconstruction of f(x) is determined by examining 
the power spectrum of the experimental data. Using the observation that the 
power spectrum falls steeply to a minimum before rising due to the white noise 
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coefficients at this first minimum. The logic behind this approach is rather 
simple. Retaining fewer coefficients than this implies a smoothing of the 
deconvolved spectrum and the loss of potentially valuable structural informa- 
tion. Adding more coefficients gains nothing because the s1gnal-to-noise ratio 
is of order one at the minimum. In addition, truncating the coefficients at 
the point where the white noise begins to dominate forces the dispersion of 
f(x) to remain finite. 

With this brief introduction it is now possible to list the six steps 
which together constitute the FFT/power spectrum deconvolution method: 

1) Both the raw data h(x) and the resolution function r(x) are 
background corrected. 

2) The functions h(x) and r(x) are Fourier transformed using an 
appropriate FFT algorithm. 

3) The Fourier coefficients of h(x) are examined and those up to and 
including the first minimum in the power spectrum of h are 
retained, the rest are replaced by zeros. 

4) The Fourier transform of f(x) is obtained by the point by point 
division of the Fourier transform of h(x) by the transform of r(x). 

5) f(x) is obtained by the inverse FFT. 

6) The solution is checked by forming the convolution of f(x) with 
r(x) and comparing the result with h(x). 

Each of these basic steps will now be described in further detail. 
Semiconductor detectors usually produce asymmetric y-ray peaks with a non- 
linear underlying background. This is caused primarily by three processes: 
1) incomplete charge collection. 2) electron escape from the active volume, 
and 3, scattering of ways (Chaglar et al, 1981). This background component 
should^ ^ 
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the deconvolution procedure 1s performed. The background correction routine 
employed in this work involves a least squares fitting of a linear function in 
both the low and high energy tail regions. The region on the low energy side 
of the peak spans from 486 to 497 KeV, while the high energy region extends 
from 634 to 542 KeV. The straight line segments fitted in these regions are 
extrapolated to SU3 and 518 KeV respectively where they are joined by a third 
line segment to produce a simple background approximation (Jackman et al « 
1974). The linear least squares analysis algorithm used in this work is a 
simplification of a higher order least squares routine by Gerald (1978). As 
in the work of Jackman et al . (1974), use is made of only the high energy half 
of the background corrected spectrum since it is known that incomplete charge 
collection of solid state detectors results in low energy tailing. In theory 
the spectrum should be symmetric about 511 KeV since statistically one should 
observe equal numbers of red and blue Doppler shifted Y-rays. 

Before proceeding to the details of the FFT algorithm one important 
point requires explanation. In order to correlate MCA channels with a cor* 
responding energy range, one must know the MCA energy resolution per channel 
as well as which channel corresponds to the center of the energy distribution 
and hence 51 1 KeV. The former value is determined as described in the experi- 
mental procedures. However, the center of the distribution must be calculated 
for each experimental spectrum. A simple centroid calculation is unacceptable 
since the higher background on the low energy side would bias the result* 
Similarly, simply selecting the channel with the maximum number of counts is 
undesirable since, due tn the statistical nature of the experiment, the 
centroid channel and the peak channel are not necessarily the same. The tech- 
nique used in this work is to determine the eleven channel region which con- 
-jAj 1ns th JL 1 a Hj est " umber of counts and j :a11 t he m iddle channel of t his region 
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the center of the distribution. In most cases, the channel corresponding to 
the center as calculated above was also the channel with the maximum number of 
counts. 

Proper arrangement of the data string Is required before proceeding to 
the FFT algorithm. Since this algorithm requires periodic data, the raw data 
are rearranged in the following way. The first datum point of the rearranged 
set corresponds to the center of the raw data string. The remaining rear- 
ranged data points numbered from 2 to ((N/2) + l) correspond to the remaining 
values of the background corrected high energy half of the raw spectrum. It 
should be pointed out that N must satisfy the relationship N = 2" 1 where m is 
an integer. Since the real data string rarely has exactly ((N/2) + l) data 
values it 1s usually necessary to insert a string of zeros at the end of the 
data set to fulfill this requirement. In order to achieve periodicity, data 
point N should correspond to data point 2. In general the points from ((N/2) 
+ 2) to N are defined in the following way: 

9((N/2) + P) = g((N/2) + 2 - P) (3<21) 
where P varies from 2 to (N/2). In the experiment performed in this work m 
was equal to eleven which corresponds to a total of 2048 data points (includ- 
ing the zerofill). 

Once the functions h(x) and r(x) have been background corrected and 
properly arranged they are passed to the FFT algorithms. The FFT algorithms 
employed in this work are those developed by Cooley, Lewis, and Welch (1979). 
There are a number of other such algorithms available which should work 
equally well. In every case it is recommended that the Fourier coefficients 
be calculated with double precision accuracy. The advantage of using FFT's in 
place of discrete Fourier transforms (DFT) can be seen from the following 
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tipUcations while the FFT requires 2N Log 2 N or ~ 4.5 x 10* complex multipli- 
cations. In this example the number of calculations required is reduced by 
about two orders of magnitude. 

The justification for using only the Fourier coefficients up to and 
including the one corresponding to the first minimum in the power spectrum of 
h(x) has been briefly discussed. Similar techniques are commonly employed in 
electronic signal processing and it can be shown that the number of coeffi- 
cients retained is a measure of the frequency resolution. In the present case 
the resolution refers to the width of the lineshape components which can be 
resolved. The procedure by which the computer calculates the number of re- 
tained Fourier coefficients from H(w) is as follows. Values associated with 
the power spectrum are obtained by multiplying each Fourier coefficient H(K) 
by its complex conjugate. However, since the function h(x) is even, the 
Fourier coefficients H(w) are all real. Truncating H(w) at the K correspond- 
ing to the first minimum in the power spectrum is equivalent to the require- 
ment that the set of retained coefficients be monotonical ly decreasing, which 
is easily programmed by testing 

H(w + 1) < H(w) (3.22) 
In this analysis only the first (N/2) + 1 Fourier coefficients are considered. 
The first coefficient corresponds to the D.C. component of zeroth frequency 
and the N/2 + 1 coefficient corresponds to the highest frequency available. 
It is well known that the remaining coefficients correspond to the negative 
frequencies and are in fact the complex conjugates of the first (N/2 + 1) 
coefficients following the relationship g(N/2 + P) = g*(N/2 + 2 -P) for P = 2 
to N/2. However since all the coefficients are real, no new information is to 
be gained by examining the negative frequencies. 
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Once the appropriate number of coefficients, say I. has been deter- 
mined the remaining N-I coefficients of H(w) are set equal to zero before 
dividing H(w) po1nt-by-point by R( w ) and proceeding to obtain f(x) by the in- 
verse Fourier transform. The consequences of this truncation, or filtering 
operation are well known. See, for example, the discussion by Verkerk (1982). 
In general the process results in some systematic errors in which the peak 
values are underestimated and nonphyslcal oscillations may appear in the de- 
convoluted function f(x). For Doppler-broadened positron annihilation spectra 
taken with our instrumentation this filtering error is negligible for Doppler 
spectra having more than 40K peak counts. 

To check the accuracy of the function f(x) requires the convolution of 
this function with r(x) and the comparison of the resulting function with 
h(x). The details of practical convolution techniques have been discussed by 
Lam et aK (19B1), and the procedure used in this work is essentially the 
same. 

A complete listing of the Fortran FFT/power spectrum deconvolutlon 
program is given in Appendix F. 

Test of the Method 

In order to explore the validity of the FFT/power spectrum method the 
following test was performed. The resolution function r(x) was obtained from 
the 614 KeV gamma ray of 85 Sr> wnile the funaion h(x) ^ ^ ^ ^ 

tron irradiated high purity single crystal sapphire sample. This sample had 
been irradiated in EMU at 740»C to a fluence of 3 x 10*1 neutrons/cm* with E 
< 0.1 MeV. Both r(x) and h(xj were then processed as described in the prev- 
ious section. For this test of the method f( x ) was estimated from F(w) using 
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tral analysis. 

Figure 5 is a plot of the power spectrum for both F(w) and H(w). The 
filled circles represent the spectrum of the experimentally obtained function 
h(x). Note the rapid decrease of about six orders of magnitude in power with 
increasing, frequency. This power spectrum is about what would be expected 
from a wide-band-width stochastic process contaminated with white noise. This 
white noise component is associated with the band of frequencies to the right 
of the first minimum in the power spectrum of h(x). Although the overall 
signal-to-noise ratio is very high, for frequencies 1n this white noise band 
the signal-to-noise ratio is less than one and therefore it is impossible to 
retrieve any structural information about the function f(x) in this high fre- 
quency range. This suggests that only those coefficients to the left of the 
white noise band, that is to the left of the first minimum, be retained. To 
use less than this number of coefficients in reconstruction is undesirable 
since it is known that the errors associated with truncation, or filtering, 
increase with a decrease in the number of coefficients used. 

The amplification of the statistical errors in h(x) by the deconvolu- 
tion procedure can be clearly seen in Figure 5. The open circles represent 
the spectrum of the deconvoluted function f(x). Note the drastic increase in 
the power associated with the high frequencies of f(x) as compared with the 
corresponding frequencies in h(x). However, the low frequency coefficients to 
the left of the first minimum are well behaved. These observations are the 
foundation on which the FFT/power spectrum method is based. 

It is not surprising to find that the number of retained Fourier re- 
construction coefficients, as determined by power spectrum analysis, is a 
function of the number of peak counts in the raw data h(x). Figure 6 is a 
semi -log plot of gUt^shlfi f^^^^j^^Q ^^Jh U figure.sug- 
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Figure 5. A Plot of the Power Spectrum of the Raw Data, h(x), (Filled 
Circles), vs. the Deconvoluted Data, f(x), (Open Circles). 
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gests that the relationship over this range of peak counts is an exponential 
of the form: 

y « m £n x (3>23) 
where y is the number of coefficients, x 1s the number of counts in the peak, 
and in is approximately equal to 7 for this data set. This curve demonstrates 
the high price one must pay in terms of data collection time to increase sig- 
nificantly the number of retained coefficients. Fortunately, our experience 
has shown that with the HUO coefficients which are usually obtainable with 
40K counts in the peak (30 minutes of data collection) the results obtained 
are very satisfactory. 

Figure 7 is a comparison of the deconvoluted function f(x) and the raw 
data h(x) after background correction. The curves have been normalized to 
equal area. As expected, the deconvoluted curve is both sharper and smoother 
than the raw data. There appears to be a slight oscillation in f(x) in the 
tail region. The quality of this estimate of f(x) can be examined in the 
usual way by convolution with r(x) to obtain an approximation h'(x) to the 
original experimental function h(x) as is shown in Figure 8. The overall 
quality of the fit is so good that on this scale it 1s difficult to see any 
differences between the raw spectrum and the convoluted spectrum. The one 
exception occurs in the peak region where several of the raw data points of 
h(x) lie well below their counterparts in h'(x). Careful examination of the 
shape of h(x) in this region reveals some rather spurious changes in the slope 
of the spectrum. Physical constraints on the annihilation process suggest 
that, assuming no posltroniura formation, the annihilation spectrum should be 
relatively smooth in this region. This suggests that these slope changes in 
h(x) are related to statistical counting errors. To test this idea two other 
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Une of the theoretical advantages of deconvolution is a significant 
improvement in the resolution of UBPAS measurements. Figure 9 is a plot of 
the high energy half of a typical 85 Sr peak before and after deconvolution. 
The closed circles represent the raw data after background correction, while 
the open circles show the shape of the resolution curve after deconvolution 
using 243 coefficients. This is the number of coefficients selected by the 
FFT/PS method corresponding to the first minimum In the power spectrum for 
this data set. Note that the full width half maximum (FWHM) value for the raw 
data is, -1280 eV while the deconvoluted curve has a FWHM of -320 eV. This 
corresponds to a factor of four improvement in resolution. However, this 
result is somewhat misleading in that for aluminum the number of retained 
coefficients is approximately 60-70. The triangles in Figure 9 show the 
resolution function after deconvolution using 70 retained coefficients. In 
this case the FWHM is -576 eV corresponding to an improvement in resolution by 
a factor of -2.2. Clearly it can be seen that the improvement in resolution 
obtained by deconvolution using the FFT/PS method is a function of the number 
of retained coefficients, which in turn 1s directly related to the statistical 
quality of the raw data set (Schaffer et al. 1984). This Implies that the 
data collection process should be carried out for sufficiently long times in 
order to obtain good counting statistics. (The counting error is proportional 
to 1/VN, where N is the number of counts in the channel). 

Unfortunately, counting for long times introduces another error 
source. This error source is related to instrumental electronic drift due, at 
least in part, to temperature fluctuations. There are two major ways in which 
other researchers have dealt with this problem. First, the accuracy of re- 
sults can be improved by operating the entire OBPAS system within a temper- 
-aAure^conj^ 
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conjunction with the first, is to digitally stabilize the electronics. This 
digital stabilization is accomplished by simultaneously monitering the 511 KeV 
annihilation line along with a suitably located monoenergetic gamma ray source 
such as the 478 KeV gamma ray from he. In this procedure a spectrum is 
collected for -100S and then the centrold of the 47b KeV peak is calculated. 
Subsequently, several other spectra are collected and their corresponding 
centroids also calculated. Next, each spectrum is linearly shifted such that 
all of the 478 KeV centroids coincide. Finally the 511 KeV peaks are summed 
to produce a DBPAS spectrum with as many total integrated counts as is 
necessary to obtain suitable counting statistics. 

While temperature and/or digital stabilization should improve the re- 
producibility of DBPAS measurements, it is still possible to obtain meaningful 
results using unstabilized equipment. In this work neither of the stabiliza- 
tion techniques were employed, and the major consequence of such a procedure 
is an observable spread in the values of the lineshape parameters calculated 
for repeated DBPAS measurements of the same metallurgical sample. This prob- 
lem will be discussed in greater detail in the sections of this chapter deal- 
ing with the separation of multiple defect types. 

In order to minimize the change in the metallurgical defect structure 
of the samples as a function of time, all of the DBPAS measurements in this 
work were made at 77°K. Performing the measurements at liquid nitrogen 
temperatures raises an important theoretical question: is the assumption that 
the positrons are essentially at rest, with respect to the electrons, at the 
time of annihilation still valid? Previous experiments by Kubica and Stewart 
(1975) demonstrate that even at temperatures as low as 4.2°K the positrons 
reach thermal velocities prior to annihilation. 
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llneshape at 77°K an experiment was performed in which six D8PAS spectra were 
obtained for a well annealed pair of aluminum samples. Three of the measure- 
ments were made at 77°K and three at room temperature. The small deviation 
(of order between the six spectra showed no correlation with measurement 

temperature. It was therefore concluded that decreasing the sample tempera- 
ture from room temperature to 77°K would not affect the positron-electron 
annihilation procedure in any experimentally observable manner. 

As mentioned in the experimental procedures the samples examined in 
this work were high purity aluminum disks heat treated to control defect type 
and concentration. It is assumed in the remainder of this work, without 
direct TEM evidence but based on the aryuments presented In the experimental 
section, that one group of samples contained predominately monovacandes while 
a second group of samples possessed a defect structure dominated by prismatic 
dislocation loops. For each pair of disks three or four annihilation energy 
spectra were obtained, from which an array of lineshape parameters were sub- 
sequently extracted numerically. The lineshape parameters examined In this 
work include the first five even central moments of the energy distributions 
defined by the relationship: 

m(I) i| U i " ' X)1 f(x 1 } < 4a) ' 

where H is the number of channels in the spectrum 

x is the channel corresponding to the spectrum centroid 
f(xj) is the number of counts in channel i 
and 1 1s a " even integer (2,4,6,8,10) corresponding to the first 

through the fifth even central moments 
A second group of parameters is defined by the relationship: 
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i=j 

i=x 



with j defined as j = x + 5(1-1) 

where f( Xi ) is the number of counts In channel i 

is the channel corresponding to the spectrum centroid 
and I is an integer from 1 through 12 

for example, ,(,, 1s the rat10 „ f t „ e nu-)er of coums ^ fjrst ^ 
channel regions on the Mgb energy sld e of the centroid (actuaUy )t 1ne)(Wes 
the centroid, to the total nu*er of counts ,„ the high energy half of the 
spectra. Sim„ar,y. x<12) is the ratio of the nunfcer of counts In the 
twe.fth f,»e channe. Inter,,! (to the right of the centroid, to the tota, 
number of counts in the high energy half of the spectrum. 

The third, and final, group of , fneshaps paraTOters cxmim<1 ape ^ 
conventional "s" parameters (Mackenzie et a,. IW0) . These ftnmtm m 
defined as the ratio of the number of counts ,„ some centra, energy region to 
the tota, number of counts in the energy spectrum. The width of the energy 
window used to define the s parameter is se,ected so that the »a,ue of 5 „ ,„ 
the range 0.4 - 0.6 (Smedskjaer and F,u S s .983). ,„ this study three 
S-parameters (S(1344). SU472). and S(1600)) were examined with respective 
central energy windows of 1.344, 1.472, and 1.600 KeV. 

These parameters were calculated for both the background corrected raw 
data and for the decon.oluted data. Hence a total of forty lineshape 
parameters were calculated for each UBPAS measurement. 

Figure 10 is a plot of the parameter x(l) for the deconvolved data. 
The open circles represent the values of this parameter for the "as guenched" 
samples (predominately monovacancies). The fined circles correspond to the 
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low temperature anneal samples (prismatic dislocation loops), and the tri- 
angles represent the samples after the high temperature anneal (the essen- 
tially defect free state). Note that the numerical value of x(l) increases 
with quenching temperature for both defect types. This trend is in agreement 
with theoretical predictions since an increasing defect concentration results 
in a higher probability of annihilation with lower kinetic energy conduction 
electrons, corresponding to a more narrow energy spectrum, and hence a larger 
value for x(l). The equations used to predict the theoretical defect concen- 
trations for both the vacancies and the prismatic dislocation loops as a func- 
tion of quenching temperature are discussed in Appendix 6. 

Figure 10 demonstrates that DBPAS measurements are sensitive to the 
presence of both vacancies and dislocations in high purity aluminum. This 
result is, of course, a necessary prerequisite for the possibility of DBPAS 
measurements being able to distinguish between the defect types. From this 
plot it is possible to perform some semiquantitative analysis and estimate 
some physical properties of aluminum. For example, one can estimate values of 
the "cross section" for the trapping of a positron by both vacancies and 
dislocations in aluminum. The accuracy of these estimates is limited by the 
reproducibility of the DBPAS measurements as well as the need for some 
independent information such as the lifetime of positrons as a function of 
defect type. 

In 1972 Cotterill et aU published a paper in which they examined the 
average positron lifetime, x avg , on a set of aluminum samples very similar to 
the ones studied in this work. Their work clearly demonstrated that positron 
lifetime experiments were sensitive to both vacancies and dislocations in 
aluminum. Before discussing the similarities between their work and the 
•-Syxient ^ experiments, it is^ import 
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the parameters T avg and x(l). I„ a review article by Sieger (1980) n is 
shown that: 

T avg s (1 - At)T 5 ♦ Att t (4#3) 
where t 5 is the lifetime in the bulk material 
T t is the lifetime in the trapped state 
and A t is the probability that the positron annihilates from the 

trapped state 
Similarly, Sieger demonstrates that: 

F " (1 - A t )F b ♦ Vt (4-4) 
where F is any characteristic shape parameter obtained from the analysis 
of OBPAS measurements 
F b ^ the value of the parameter in the bulk state 
and F t is the value of the parameter in the trapped state 

Since the parameter x(l) can be shown to be equivalent to an - S - type 
parameter, and since an «S" type parameter is a common characteristic line- 
shape parameter, Eq. (4.4) can be rewritten as: 

X(l) - (1 . At)x( i, b + AtX(1)t (4 b) 

where x(l) b Is the value of the parameter in the bulk state 
and x(l) t is the value of the parameter in the trapped state 
Clearly a comparison of Eq. (4.3) and Eq. (4.5) demonstrates the similarity 
between x avg and x(l). 

With this in mind it is not surprising to find that plots of r avg vs. 
quench temperature (T Q ) and x(l) vs. T Q share many common characteristics. 
(The t avg plot used for comparison is Figure 2 in the work of Cotterill et 

1972). In both figures the Tg at which the vacancy concentration is 
large enough to affect the measured parameter is ~20U'C. The saturation value 
° f b ° th >y ^ d X <- 1 L f0r vacancies occurs at„T u , ■ 600°C. , The curve.throuq* 
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the points corresponding to vacancies is sigmoidal in shape for both the 
lifetime and deconvolved DBPAS measurements. As for the points corresponding 
to samples with prismatic dislocation loops, the value of both parameters 
initially deviates from their bulk values for samples annealed after a quench 
from ~400<»C. Neither parameter reaches a saturation value for dislocations 
within the temperature range examined in this work. 

Figure 10, therefore, shows not only the ability of deconvoluted DBPAS 
measurements to detect the trapping of positrons by both vacancies and dislo- 
cations in aluminum, but also qualitatively verifies the theoretically pre- 
dicted equivalence of the lifetime parameter T avg and the Doppler parameter 
x(D. 

Figure 11 is a plot of the same lineshape parameter. x(l), vs. Tq for 
the background corrected raw data. i.e. the same data as shown in Fig. 10 but 
before the deconvolution procedure. Note the general similarity between these 
two figures. Careful examination of the raw data shows that the points 
corresponding to the 400°C quench are systematically displaced. That is, the 
entire group of five data points seems to be "too high" as compared with the 
rest of the data. To understand the cause of this shift it is necessary to 
examine the shape of the instrument resolution function R(x). 

During the course of these experiments the instrument resolution was 
periodically checked by collecting an energy spectrum for 85 Sr# Tne mnoener . 
getic line is broadened due to the physical characteristics of the DBPAS sys- 
tem as previously described. This broadened energy line serves as the resolu- 
tion function. K(x), in the FFT/PS deconvolution procedure. It also describes 
the amount by which the undistorted energy spectrum is broadened by the in- 
strumentation. 
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Figure 12 shows two such resolution spectra. The open circles repre- 
sent the instrument resolution corresponding to the DBPAS measurements for the 
samples quenched from 400»C. The filled circles represent a resolution func- 
tion which is characteristic (within statistical error) of the resolution 
functions for all of the remaining data points. This comparatively large 
shape change in the function R(x) was apparently caused by a power failure in 
the laboratory. Note that the Function R(x) corresponding to the 400«C 
samples (henceforth termed SR8S0LD) is more sharply peaked, and therefore more 
narrow, than the other resolution function (SR85NEW). 

The physical significance of the change in resolution function is as 
follows. An observed energy spectrum will be broadened in two ways. The 
first way is by a decrease 1n the concentration of vacancy type defects, in 
the sample being examined. The second broadening mechanism Is, of course, the 
broadening due to the instrumentation. Now. if the second source of broaden- 
ing remains essentially constant, then changes in the width of the observed 
spectra may be attributed to changes in defect concentration. On the other 
hand, if the broadening due to the instrumentation decreases, and this error 
source is not properly removed, the observed spectrum will be artificially 
narrow and its apparent defect concentration will be too great. Thus, if the 
data for the 400°C samples are collected with an artificially narrow resolu- 
tion function, i.e. SR8501D, their apparent defect concentration will be 
greater than they should be. 

Of particular significance is, that while a change in experimental 
conditions can alter the shape of the resolution function and thereby distort 
the values of lineshape parameters calculated for the raw data, the deconvo- 
lution routine, when supplied with the appropriate R(x) function, auto- 
? a M<^!y thaLfprj-eerodjica-,. 
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bility the deconvolution of DBPAS measurements should be considered. In 
addition, as mentioned previously, since deconvolution removes the Instru- 
mental effects it should yield lineshape parameters which are directly 
comparable between different laboratories. 

Despite the advantages of deconvolution it 1s important to remember 
that this technique does, of course, have its limitations. In UBPAS 
measurements there are several sources of experimental error, each of which 
has a corresponding time scale. The magnitude of this time scale is the 
overriding factor in determining the applicability of deconvolution techniques 
in minimizing the specific error in question. For example, any change in 
experimental conditions which can be treated as a step function, that is, by 
changing from one constant value to another constant value in a very short 
period of time, can be easily treated by deconvolution. Such "error" sources 
include the accidental (e.g. power failure) or intentional changes in the 
electronic amplification of the energy signal. To correct for. such a change 
one simply needs to collect a resolution function before and after the event 
and use the appropriate R(x) in the deconvolution algorithm. 

On the other hand, error sources having a time scale approximately 
equal to the time required to collect a data set (tens of minutes) can not be 
corrected by deconvolution alone. Examples of such possible error sources 
include line voltage fluctuations, temperature variations, and electronic 
drift. These error sources may cause the centroid of the resolution function 
to vary during the collection period, resulting in a composite R(x) which is 
artificially broad. In addition, even if the raw data set to be deconvoluted 
and the function R(x) are collected in rapid succession the amount of centroid 
drift for the two energy spectra may well be uncorrected. This type of error 
^r^^^^^^ent-^ns of correction. The DBPAS measurements 
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should be performed in a temperature controlled laboratory. Also, digital 
stabilization, as described previously, should be used to control centrold 
drift. 

The final error type is the small time scale or high frequency class. 
An error source such as sixty-cycle electrical noise, if it were present, 
should average to zero over the time required to collect a data set. In 
addition, a high frequency error source should have a similar effect on both 
the resolution function and the DBPAS eneryy distribution. 

In order to judge the effectiveness of deconvolution in eliminating 
experimental error it is necessary to define some quantitative measure of this 
error. The predominant manifestation of experimental error appears as a 
spread in lineshape parameter values calculated for multiple DBPAS measure- 
ments on the same samples. Figure 13 is a plot of the parameter x(l) vs. Tg 
for the deconvoluted data, note the spread in values of x(l) for a given 
defect configuration. 

In an attempt to evaluate the performance of the deconvolution algor- 
ithm the following parameters are introduced (see Figure 13 for a physical 
interpretation). 

1) AP V /P 0 is a measure of the sensitivity of this parameter to the 
existence of vacancies. 

2) 6p y /AP v is a measure of the error associated with the spread in 
values of this parameter for vacancies. 

3) aPj/Pq is a measure of the sensitivity of this parameter to the 
existence of dislocations. 

4) SPjyaPj^ is a measure of the error associated with the spread in 
values of this parameter for dislocations. 
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Clearly, it is desirable for the sensitivity (AP/P 0 ) to be as large as 
possible and the error (6P/AP) as small as possible for both defect types. A 
comparison of these quantities calculated from the deconvoluted and the raw 
data is shown in Table 1. Note that for all four of the lineshape parameters 
listed the sensitivity ratio is higher and the error ratio is lower for the 
deconvoluted data than for the raw data. This analysis encourages the further 
use of the deconvolution procedure to improve the quality of the results ob- 
tained from OBPAS measurements. 

A tentative explanation for the apparent improvement in sensitivity 
and accuracy of the deconvoluted results is based upon trends in the shape of 
the energy spectrum before and after the application of the FFT/PS analysis. 
Figure 14 shows the high energy half of the peak region for an A* sample 
quenched from ~60U«>C. The deconvoluted spectrum is, of course, more narrow 
than the raw data but note that it is also smoother (due to the filtering of 
the high frequency components). One explanation for the comparatively large 
error ratios of the raw data is that they are related to the statistical 
counting errors 1n the peak region of the spectrum. If this is indeed the 
case it is not surprising to observe a decrease In the error ratio after the 
filtering of the high frequency components. However, the most probable 
explanation for the increased sensitivity resulting from deconvolution is 
related to the removal of the "broadening" action of the experimental instru- 
mentation. Recall that the major effect of the instrumentation is to spread 
counts from one channel into several adjacent channels. Therefore, the true 
lineshape changes in the annihilation energy distribution, caused by the 
trapping of positrons at defects, are spread out over a much larger number of 
channels in the raw data. The deconvolution procedure corrects for this 
effect and in so. doing restores., some of the details of the- M true" annihilation 
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Table 1 



-2SiiJ x (2) S(1344) $(1600) 



^ EC V JT" .0360 .0263 .0325 

o 

DEC, Hffl .315 .320 .331 



.215 
.0279 

.581 



UEC i - 0124 -0110 .0119 .0106 
6P 

ND v TT -235 .235 . 2 31 

ND V ^ -0348 .0165 .0256 .0207 

6P 

ND 1 -ST 21 - 471 .938 .555 .585 

ND l JT .0155 .0077 .0100 .0078 
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energy distribution. Th e net result is an increase In defect sensitivity. It 
should be mentioned that the magnitudes of the enhancement in sensitivity and 
accuracy due to deconvolution noted in this work are consistent, but small. 
However, from this work alone it is not possible to estimate the effect of 
temperature and electronic drift upon the results. It would be extremely 
enlightening to repeat this analysis employing a temperature controlled and 
digitally stabilized system. 

Before moving on from the topic of parameter sensitivity and accuracy, 
one more pertinent observation should be discussed. Figure 15 is a plot of 
the error ratio. 6P/AP. calculated separately at each quenching temperature 
(the parameter being examined is x(l)). Several trends are apparent. First, 
it is clear that the error associated with the dislocations is larger than the 
error associated with vacancies. This trend may be a. direct consequence of 
the size of the denominator. .P. which is smaller, at a given quench tempera- 
ture. for samples with dislocations than for the corresponding samples with 
vacancy type defects. This observation is however somewhat deceptive in that 
the horizontal axis of Figure lb is T Q and not defect concentration. A mean- 
inyful comparison of parameter sensitivity between defect types should be made 
"sing equivalent defect densities. Unfortunately, the dislocation densities 
which can be introduced by this thermal treatment are not large enough for 
direct comparison. In addition, it should be mentioned that the trapping of 
positrons by dislocations may be somewhat more complicated than trapping by 
vacancies. This idea is discussed in detail by Martin and Paetsch (1972). 
Very briefly, it may be that the depth of the potential well associated with a 
prismatic dislocation loop is such that detrapping of the positron, due to 
local thermal excitations, becomes significant. If this is indeed the case, 
the interpretation of DBPAS .data for low prismatic dislocation loop, densities 
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becomes more complicated. Clearly, additional experimental data Is required 
to determine the exact reason for the observation that the error parameters 
are laryer for the samples containing primarily dislocation loops than for the 
samples dominated by vacancy type defects. 

Another interesting feature of Figure 15 is the relatively constant 
value of 6P/AP for vacancies. This implies that the spread In the data de- 
creases with the quenching temperature, since clearly AP decreases with T Q . 
One possible explanation for this observation is related to the shift in the 
centrold of the experimentally obtained spectrum due to temperature and/or 
electronic drift. The explanation is that a centrold shift changes the peak 
shape of a broad and flat spectrum (low defect concentration) to a lesser ex- 
tent than that of a narrow and sharp spectrum (high defect concentration). 
This is, of course, an extremely qualitative argument but it does predict that 
in a temperature/digitally controlled series of experiments the relative size 
of the error bars associated with the high defect concentration samples should 
decrease. 

The preceding discussion has demonstrated the utility of deconvolution 
in DBPAS measurements. In summary, the use of the FFT/PS method improves the 
sensitivity and the relative error of the calculated lineshape parameters. In 
addition, deconvolution automatically corrects for a wide range of error 
sources which result in changes in the instrument resolution function. Decon- 
voluted lineshape parameters should be more directly comparable between dif- 
ferent laboratories and. in addition, offer the potential ability to resolve 
fine structure in the energy spectrum which 1s unavailable from direct obser- 
vation of the raw data. There is strong evidence (see for example Hood and 
Schultz, 1980) suggesting that the use of both temperature stabilization and 
digital stabilization, to control centrold dri ft, .should. significant 1y,.im P rove 
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It is Interesting to note the conditions under which the undeconvo- 
luted raw data has yielded reasonable results in the past. The main require- 
ment is, of course, experimental conditions which yield a nearly constant 
resolution function, R(x). As one would expect, temperature and digital con- 
trol as well as performing the series of experiments as close together in time 
as possible, both tend to help yield reasonable results for the raw data. In 
addition, the use of "S" type parameters, which are essentially area integrals 
over large energy windows, will be affected only slightly by the "spreading" 
action of the instrumentation. The major advantages one gains by deconvolu- 
tion are: 1) the ability to compare more accurately results for a wide var- 
iety of experimental conditions and, 2) the possibility of observing a more 
detailed structure of the positron-electron annihilation energy spectrum. 

The Separation of Multiple Defect Types 

The previous discussion has demonstrated that deconvoluted D8PAS 
measurements are able to detect the presence of both vacancies and disloca- 
tions 1n high purity aluminum. The evidence for this statement comes from 
controlled laboratory experiments using samples with reasonably well charac- 
terized defect distributions. In order for DBPAS measurements to become use- 
ful as. an NOT tool in industrial applications it is important for the tech- 
nique to be able to distinguish between multiple defect types. The reason for 
this requirement is that in a practical application one may require informa- 
tion about a particular defect type while being comparatively unconcerned with 
other defect types. For example, if one is interested in hydrogen embrlttte- 
ment of high strength steels the "masking" of vacancies by nascent hydrogen 
may be of primary importance while slight variations in dislocation density 
may have relatively little ^significance. "' ~ 
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Before discussing the results of this investigation to determine the 
ability of DBPAS to separate defect types a brief review of previous work in 
this area will be presented. The basic concept is that no two defect types 
will have identical electron configurations associated with them. That is, 
the distribution of electrons (both in terms of electron density and momentum) 
in the neighborhood of a single vacancy will be different than in the neigh- 
borhood of a dislocation core. Qualitatively the atmosphere will be the same, 
i.e. a decrease In electron density and a lower average electron momentum, but 
quantitatively, the measured lifetimes and fine structure of the momentum dis- 
tribution should be different. 

This variation in electron distribution has been verified repeatedly 
in the case of lifetime measurements. For example, Fluss et a]_. (1978) report 
a lifetime of -244 ps for positrons trapped at a vacancy in aluminum while 
Hautojarvi et al_. (1970) give the value of -229 ps for a dislocation in the 
same material. It should be mentioned that the accurate determination of 
multiple lifetimes from a single experiment requires some rather involved cal- 
culations. The complexity of the procedure Increases rapidly with the number 
of unconstrained components one attempts to resolve and, for crystalline 
materials, the results are questionable if more than two defect types are 
considered. 

An alternate method of analyzing lifetime data is to calculate a 
single lifetime component, known as the "average" lifetime, for a sample con- 
taining multiple defect types.' This approach simplifies the calculations and 
finds its greatest utility in experiments in which the concentration of one 
defect type is varying over a wide range of values while the concentrations of 
the other defect types remain relatively constant. Of course, such an 
.. approach results in the loss of the ability to separate defec t types. 



76 

In the case of DBPAS measurements there have been fewer attempts to 
distinguish between different defect types. The most notable exception is the 
work of Mantl and Triftshauser (1978) which will be discussed in a following 
section. There are at least two major reasons for the lack of work In this 
area. The first reason is related to the comparatively poor energy resolution 
of the DBPAS technique. Because of the smoothing action of the instrumenta- 
tion the fine structure of the energy distribution may have been hidden in 
previous experiments. In addition, even a slight change in the instrument 
resolution function during a series of experiments could mask the subtle 
cnanyes 1n lineshape required to separate defect types. Although temperature 
and digital stabilization of the electronics would improve the situation, such 
a procedure does not address the issue of the smoothing action of the instru- 
mentation. It is the use of deconvoluted data which offers the possibility of 
examining the fine structure of the electron momentum distribution. 

As mentioned in the introduction, the most comnonly reported Doppler 
lineshape parameter is the Mackenzie "S" parameter (1970). This Uneshape 
parameter is defined as the ratio of counts In some central region of the 
DBPAS energy spectrum to the total number of counts in the spectrum. Since it 
is known that the energy spectrum becomes more sharply peaked with increasing 
vacancy type defect concentration 1t follows that the numerical value of M S" 
should be proportional to the total vacancy type defect concentration. The 
disadvantage of using this lineshape parameter is that it is unable to 
distinguish between multiple defect types. For example, if both monovacancies 
and edge dislocations are present in a sample both defect types will tend to 
increase the value of "s\ however, it will be impossible to separate the 
relative defect concentrations solely from the calculated value of "S". 
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Mantl and Trifthauser (1978) developed a method for separating the 
effects of multiple defect types from experimentally obtained OBPAS data. 
To accomplish this they introduced a new lineshape parameter, I c , defined 
as the ratio of the sum of the number of counts in two regions, (one on the 
low energy side and one on the high energy side of the peak) divided by the 
total number of counts in the spectrum. This lineshape parameter, which in 
this work will be referred to as the winy or "W" parameter, can be shown to 
decrease with increasing vacancy type defect concentration. When analyzed 
separately the "W" parameter suffers from the same disadvantages as the "S" 
parameter in that it is incapable of separating multiple defect types. How- 
ever, Mantl and Trifthauser recognized that the combined use of "S" and "W" 
could yield useful results. 

In Appendix H the derivation of their defect specific "R" parameter is 
recreated. To summarize, it can be shown that in theory the "R" parameter, as 
defined below, is concentration independent, but characteristic of the type of 
trapping site Involved (Mantl and Trifthauser 1978). The "R M parameter is 
defined as: 

R - |(S - S„)/(W - W B )| (4.6) 
where S is the "S" parameter for a given specimen, 

S B is the value of "S" corresponding to a "defect free" specimen, 

W is the "W" parameter for a given specimen, 

W B is the value of "W" corresponding to a "defect free" specimen. 
While such a parameter is very useful, because of its ability to distinguish 
between multiple defect types, it still does not solve the problem completely 
due to some inherent theoretical limitations. The main restriction is that 
the analysis is valid if only one type of trapping center is present, or if 
one defect type is essentially dominant (Mantl and Trifthauser 1978). In 
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practice this means the experimenter has the ability to calculate defect type 

and concentration provided each sample has only one predominant type of 

positron trap. The difficulties arise if a sample contains two defect types. 

say vacancies and dislocation cores, in roughly equal concentrations. 

In terms of the present experiments on high purity large grain size 

aluminum, because care has been taken to produce samples with only a single 

dominant defect type, one would expect the "R" parameter analysis to yield 

useful results. A family of R-type parameters can be calculated for the 

experiments in this work by examining ratios of the form 

R N ' l(Ux(1)/Ax(j))| (4.7, 
where R N is an R-type parameter 

Ax(i) 1s the difference between x(i) for the current sample and the 
value of x(i) in a defect free sample 
and i is less than j 

In the original work of Mantl and Trifthauser (1978) the energy windows 
selected to define "S" and "W were chosen so that the annihilation events 
contributing to the "S M parameter corresponded to conduction (valence) 
electrons while core electron events resulted in contributions to the M W" 
parameter. This relationship restricts the possible values of 1 and j as 
defined in Equation (4.7) above. 

The specific "R" parameter used in this study is defined as 

p - I Ax(l) I 

1 ' Mx(lU) + x(ll) + x(12)J| (4.8) 
and is defined such that the energy windows are fairly close to those used by 
Mantl and Trifthauser (1978). For comparison, Rj has a low energy window of 
-till - 511.32 KeV and a high energy window of -513.88 - 514.84 KeV compared to 
windows of -511 - 511.5 KeV and -513.6 - 516.1 KeV in the referenced work. 
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Figure 16 is a plot of Ki as a function of quenching temperature for 
the background corrected raw data. Note that although samples with primarily 
dislocation type defects have, in general, a higher R x value than the samples 
with primarily vacancy type defects, there is a significant amount of scatter 
in the data. The source of this scatter is not known exactly, but it is prob- 
ably related to the fact that the data was collected without temperature or 
digital stabilization. In fact, it is plausible that the spread in Rj values 
Is a numerical amplification of the spread in values of x(l), x(10), x(ll) and 
x(12). In addition to this error source, any error associated with changes in 
the instrument resolution function, (IRF), should also be present in Figure 
16. To test the utility of deconvolution in R-parameter analysis one would 
like to have a data set in which the second error source is present but the 
first error source, related to electronic drift, is negligible. Unfortunate- 
ly, such a data set can not be obtained with the instrumentation used in this 
study. 

In an attempt to isolate the instrument resolution effects the fol- 
lowing procedure was employed. Using the fact that all DBPAS measurements for 
a given metallurgical sample were performed in a comparatively short time 
interval (-8 hrs.), it seems reasonable to assume that the spread in calcu- 
lated lineshape values is a result of the lack of temperature and digital 
stabilization rather than of changes in instrument resolution. This hypoth- 
esis is supported by the observation that for a given set of consecutive 
experiments, the IRF's measured before and after the experiments showed only 
very slight changes. If this assumption is true, one may be able to arti- 
ficially select a subset of the original data which qualitatively corresponds 
to a series of experiments with reduced electronic drift error. In this work 
each group of 08?AS. measurements, ...corresponding to a sample subjected to 
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a given thermal treatment, was independently examined and an average value for 
each relevant lineshape parameter was calculated. The corresponding R x values 
calculated from the average lineshape parameter values are shown in Figure 17. 
These values are, as in Figure 16, for the background corrected raw data. 

A similar procedure was then employed for the deconvoluted data. That 
is, average x(i) values were calculated for each specimen and then used to 
calculate R x values. These Rj values for the deconvoluted data are shown in 
Figure 18. It is quite apparent that the spread in. R-parameter values has 
been increased by the FFT/PS deconvolution algorithm. The explanation for the 
increased scatter in the data is not inmediately obvious. In a previous 
section of this work it was demonstrated that deconvolution improved the 
quality of the calculated x(l) lineshape parameter values. That is, the 
spread in x(l) values for DBPAS measurements performed on a given sample was 
decreased (refer to Table 1). One could argue that decreasing the scatter of 
x(l) for a given sample should Improve, or at least not adversely affect, the 
amount of scatter in the calculated *i values. The difficulty is, of course, 
that «! depends not only on x(l) but also on the values of the parameters 
x(10), x(ll), and x(12). These parameters are defined in the tall regions of 
the energy distribution and it is in this region that one must look to under- 
stand the difficulty associated with the calculation of R-type parameters for 
FFT/PS deconvoluted data. 

Figures 19 and 20 are a comparison of the energy distribution (raw 
data) for a sample quenched from 600°C and then subsequently annealed at 
523°K. Figure 19 shows the differences in the peak region of the spectra. 
Recall that x(l) is defined over the energy window of 511 - 511.32 KeV. 
In Figure 20 the vertical scale has been expanded so that the differences 
in. the.. taiL region are visible. Several, important features., of .the.. energy 
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20. A Comparison of the Tail Region of the Raw Data for an Aluminum 
^Sl e jM??ched from 600fC (0) and the Same Samp„l.e„Mt£r. the Hiqh 



d1.trlb.t1. 1„ th, ,.„ region are apparent. First _ jt mst ^ (hjt ^ 

oeneral th, data points for th, quenched sample „. b„o. their counterparts 
for th, annealed sample a5 pred)cted by tMry _ ^ stat)stjM] ^ 

counting, error from cha„„„ to channel can b, recognized. Th, former obser- 
vation is believed ,o account for th, qualitative agreement between th,- 
oretlcal predictions and th, calculate Rj values f„ r the „„ data . The 
l«„r obsorvation. statistical nois,. is th, most probabl, emanation ,a, 0 „ 9 
with electronic nois,) for th, spread in «, values as shown in Figure lb. 

For comparison. Figures 21 and 22 show th, chan 9 , in energy d1strft>u . 
tion for th, same two sables aft,r d,c„„»„,„„ 0 „ „ s ,„ 9 tne FfT/ps 
In Figure 21. « h ich shows th, „,,,< region, som, of ,n, d,sirabl, features of 
d,c.n,olut1,„ are evident. Most notably, the statistical channe, to chann,, 
nois, has been removed (compare with Figur, 1 9 ). A careful comparison „, the 
raw and deconvoluted data shows that devolution does indeed decreas, th, 
full width half maximum as discussed previously. 

Figure 22 shows, with an expanded vertical seal,, an undesirable 
feature of the FFT/PS devolution method. It was mentioned ,„ Chapter .„ 
that truncating the Fourier coefficients is equivalent to the insertion „f a 
low pass filter and it is well known that low pass filtering introduces 
oortain errors. «ec»l.i„ 9 that the Fourier analysis of the convolution 
,qu«1on, negating the statistical counting error. yi„ds th, following 
equation for the desired function f(x): 



f(x) = u(w)F(w)e-1wx dw (4>1()) 



where u(w) = 1 for w < w c 
= 0 for w > w c 
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Figure 22. 
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and w c corresponds to the frequency which is the first minimum of the 

power spectrum of h(x) 
It can be shown (Bendat and Piersol 1971) that error introduced by such a 
filter is proportional to (l/w c )sin(2wxw c ). The significance of this result 
is that the amplitude of the error is inversely proportional to w c (the number 
of retained coefficients), and the frequency of the resulting error is 
proportional to w c . 

This effect is clearly evident in Figure 22. The open circles 
represent the annealed sample for which the power spectrum analysis yielded a 
w c of fifty six (56) coefficients. For comparison, the filled circles, 
corresponding to the quenched specimen are calculated using sixty seven (67) 
retained coefficients. It is quite obvious that the error amplitude is 
decreased for the sample with the larger number of retained coefficients. A 
careful analysis of the data shows that the frequency of the error is 
proportional to w c . 

With this information, it is now possible to explain the increased , 
scatter in R\ values displayed by the FFT/PS deconvoluted data. Because the 
energy window used to define the wing parameter (513.88 - 514.84 KeV) is in 
the range where the amplitude of the filtering error is of the same order of 
magnitude as the amplitude of the true signal, this error source is 
significant. If the wavelength of the error source was small compared to the 
size of the energy window, then the error contribution would average to nearly 
zero and this difficulty would not exist. The problem arises because the 
energy window spans less than a period of the sinusoidal error source and it 
spans a different fraction of the period depending on the number of retained 
Fourier coefficients. For example, in Figure 22 the energy window of interest 
for the annealed sample spans approximately the. ; posit tve half cycle of the . . . 
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sinusoidal filtering error. This leads to a calculated wing parameter with an 
artificially high value, which in turn distorts the value of the calculated R, 
parameter. Since multiple OBPAS measurements on even the same metallurgical 
sample may yield slightly different numbers of retained coefficients, this 
error source can lead to a scatter of Rj parameters for a single sample. 
Because the number of retained coefficients varies more widely between samples 
of differing defect density (or defect structure), the spread 1n Rj values 
becomes larger for such samples. 

Having related the scatter In R x values for the deconvoluted data to 
the phase and amplitude of the filtering error, it is now possible to try to 
correct for this error source. There are at least three ways to proceed. 
First, one could simply extend the experimental collection times to get a 
higher number of total integrated counts in the energy spectrum. This would 
result in better counting statistics, and a higher number of retained 
coefficients (Schaffer et aU 1984). Since an increased number of retained 
coefficients decreases the amplitude of the filtering error, such a procedure 
should yield better R x values. Unfortunately one must pay a very high price 
in terms of collection time to significantly increase w c . 

Another approach is to use the same number of retained coefficients 
for each data set. That is. one could calculate a w c value for each data set 
to be analyzed and then deconvolute all of the data sets using the minimum 
value of w c . Using this procedure will not eliminate the error source, but it 
will be controlled so that its contribution to the calculated values of Rj 
does not vary between experiments. Clearly such an approach has several 
inherent disadvantages. The primary disadvantage is that the absolute. values 
of Rl will be a function of the particular value of w c used in the devolu- 
tion ..rocedure. One would prefer an approach which eliminates, or at least 
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greatly reduces, the filtering, error rather than simply controlling it. A 
method for accomplishing this goal will be discussed in the final chapter of 
this work. What follows here is a discussion of controlling the error source 
in the tail region by using the same number of retained coefficients for each 
data set. 

Figure 23 shows the energy window used to calculate the wing 
parameter. The vertical scale has been greatly expanded to show the fine 
detail of the structure. This plot is for the same samples as in Figure 22, 
but in this case both data sets have been deconvoluted using fifty seven (57) 
coefficients. As expected, the filtering error is in phase. In addition, as 
theory predicts, the data points corresponding to the quenched sample are 
below the points representing the annealed sample. One can extend this 
procedure to all of the other data sets and then calculate the corresponding 
Ri parameters in the manner previously described. The results are presented 
in Figure 24. When this figure is compared with Figure 18 it can be seen that 
the filtering error does influence the scatter in Rj values. By simply 
controlling this error source, i.e. by using the same number of retained 
coefficients in the deconvolution of all data sets, the quality of the 
calculated values can be improved. 

To summarize this section, the average values of Ri, along with the 
associated standard deviations for each of the two defect types, have been 
calculated for four different types of data. The four data types are: 1) the 
entire raw data set (Fig. 16), 2) the average raw data (Fig. 17), 3) the 
averages of the deconvoluted data as obtained directly from the FFT/PS 
algorithm (Fig. 18), and 4) the averages of the deconvoluted data using the 
same number of retained coefficients for each data set (Fig. 24). These 
values are shown in Table II. - 
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Table II 



Data Type 


Kavg( v ac.) 


R a vg(° is -) 


o(Vac.) 


o(D1s.) 


1 


.261 


.458 


.0677 


.1751 


2 


.240 


.440 


.0347 


.0492 


3 


.234 


.364 


.1515 


.1876 


4 


.253 


.424 


.0216 


.0290 



Note that the largest standard deviations occur for the deconvolved 
data when the maximum number of coefficients are used in the reconstruction of 
each f(x). This result is explainable in view of the previously discussed re- 
lationship between w c and the amplitude and phase of the filtering error. 
When this error source 1s controlled by using a single value of w c for the de- 
convolution of all data sets (data type 4) the calculated standard deviations 
are improved. With this modification to the FFT/PS deconvolution algorithm it 
appears that if one defect type dominates, one can distinguish between 
vacancies and prismatic dislocation loops, in high purity large grain size 
aluminum. 

To conclude the results and discussion chapter of this work the even 
central moments of the annihilation energy distributions will be discussed. 
As mentioned in the introduction, because this technique has the potential to 
distinguish between multiple defect types, this family of lineshape parameters 
is interesting from a theoretical point of view. Briefly, this possibility 
arises because the higher order central moments are increasingly dominated by 
the tail regions of the energy distribution. Thus, the low moments give 
information about the central region (S-type information) while the higher 
order moments yield W-type Information. Unfortunately, this analysis is 
extremely sensitive, to small errors in the tail regions which, are greatly 
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amplified by moment analysis. Therefore the statistical error ruled out 
moment analysis in the raw data while filtering error precluded its use in the 
deconvoluted data. It is possible, however, that an improved filtering 
technique, discussed in the following chapter, will allow moment analysis to 
become a useful tool in future work. 



CHAPTER V 
CONCLUSIONS AND FUTURE WORK 



Over the past several decades positron annihilation spectroscopy has 
evolved into an accepted tool for the study of the defect structure of mater- 
ials. Of the three most common measurement techniques Doppler broadening of- 
fers the advantages of reasonably high data collection rates, comparatively 
low source strength requirements, and few geometrical experimental con- 
straints. It is therefore the technique with the best chance of making the 
transition from a laboratory research device to an industrial NOT tool. Un- 
fortunately, DBPAS measurements suffer from some inherent limitations. The 
two most significant limitations are comparatively poor instrumental energy 
resolution and limited ability to distinguish between multiple defect types. 
Based on these observations, the two primary goals of this work were: 1) to 
develop an efficient, model independent, noniterative deconvolution method, 
and 2) to improve the ability of OBPAS to distinguish between multiple defect 
types. 

In order to improve the effective energy resolution of the DBPAS in- 
strumentation the fast Fourier transform/power spectrum analysis deconvolution 
algorithm was developed (Schaf fer et al_. 1984). The convolution integral 
equation was shown in Chapter III to describe the relationship between the ex- 
perimentally observed energy distribution, the instrument resolution function, 
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and the "true" or desired Doppler broadened energy distribution* Fourier 
analysis was then used to obtain what appeared to be a closed form solution 
for the desired function. Unfortunately, it was demonstrated that this 
solution is not unique. The statistical nature of the experiment places the 
equation into the mathematical class of incorrectly posed problems. In an 
attempt to find a reasonable approximation to the desired function the 
dispersion of the apparent closed form solution was investigated. It was 
observed that the dispersion became unbounded if the power spectrum of the 
statistical experimental "noise" remained finite for high frequencies. That 
is, the mathematical difficulties are related to the "white noise" component 
present in the experimental data. 

Since the shape of the power spectrum of the raw data appeared to be 
of critical importance, particularly if it contained a "white" component, it 
was investigated in greater detail. From this analysis it was observed that 
the overall signal-to-noise ratio for the experiment was very high, being of 
order -10 6 . Using the observation that the power spectrum falls steeply to a 
minimum before rising due to the white noise component in the experimental 
data, the FFT/PS method developed in this work truncates the Fourier coeffic- 
ients of the raw data at this first minimum. The logic behind this approach 
is rather simple. Retaining fewer coefficients than this implies a smoothing 
of the deconvoluted spectrum and the loss of potentially valuable structural 
information. Adding more coefficients gains nothing because the signal- 
to-noise ratio is of order one at the minimum. In addition, truncating the 
coefficients at the point where the white noise begins to dominate forces the 
dispersion of the solution to remain finite. 

The truncated Fourier transform of the raw data 1s then used to cal- 
. culate a reasonable estimate of the desired function.. _ It was shown*, in Chapter 
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III that this deconvolution algorithm does indeed yield reasonable results. 
It is well known that truncating the coefficients in the manner described 
results in the introduction of a filtering error. The consequences of the 
introduction of this error source were discussed in detail in Chapter IV and 
will be reviewed below. 

In order to test the utility of the FFT/PS deconvolution algorithm a 
series of high purity, large grain size aluminum samples, containing known de- 
fect types and concentrations, were prepared. The thermal treatments employed 
to prepare the samples containing various concentrations of monovacancy or 
dislocation type defects were described in detail in Chapter II. Basically 
the procedure consisted of a two stage quench from elevated temperatures fol- 
lowed by a low temperature anneal and finally a high temperature anneal. The 
DBPAS spectrum was measured after each of these three thermal treatments. It 
was assumed, without transmission electron microscopy but based on arguments 
presented in Chapter II, that each of the three treatments introduced a single 
dominant defect type. These defect types are: 1) vacancies, after the 
quench, 2) prismatic dislocation loops, after the low temperature anneal, and 
3) the bulk or "defect free" material, after the high temperature anneal. 

In Chapter IV the lineshape parameters of the DBPAS measurements were 
calculated (for both the raw and the deconvoluted data) as a function of 
quenching temperature, A comparison of these lineshape parameters revealed 
several interesting trends. First and foremost, it appears that the FFT/PS 
deconvolution method is capable of automatically correcting for accidental 
and/or intentional changes in the instrument energy resolution function (IRF). 
It must be mentioned that there are some restrictions on the types of changes 
in the IRF which deconvolution can correct and these restrictions are. dis- 
cussed in detail in Chapter IV. .Another advantage, of the use of deconvoluted 
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data is that results obtained in different laboratories should be more direct- 
ly comparable if the instrumental broadening effects have been accurately re- 
moved. This ability to directly compare numerical values for standard line- 
shape parameters should help to clear up some of the discrepancies which 
currently exist in the literature. 

While it is fairly clear that accurate deconvolution should correct 
for changes in the IRF and improve the agreement of results obtained using 
different instrumentation, some of the other trends observed from the compar- 
ison of lineshape parameter calculated before and after deconvolution were not 
entirely predictable. For example, based on the limited results presented in 
this work, both the sensitivity and the accuracy of the lineshape parameters 
are improved by deconvolution. One explanation for the improved accuracy, 
i.e. the decrease in the error ratio, of the deconvoluted lineshape parameters 
is that the filtering operation has essentially removed the statistical count- 
ing errors present in the raw data. A comparison of multiple DBPAS measure- 
ments performed on the same sample show that before deconvolution, differences 
in the statistical counting error 1n the peak region can lead to a spread in 
calculated x(l) values. However, after deconvolution, the energy distribu- 
tions show less variation from experiment to experiments, and hence yield x(l) 
values with correspondingly less scatter. 

The increased sensitivity of the deconvoluted parameters may be 
related to the removal of the "broadening" action of the experimental instru- 
mentation. Recall that the major effect of the instrumentation 1s to spread 
counts from one channel into several adjacent channels. Therefore, the true 
lineshape changes in the annihilation energy distribution, which are caused 
by the trapping of positrons at defects, are spread out over a much larger 
number of channels in the. raw data. The deconvolution proc.edu re .GQrr.efts. for 
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this effect and, in so doing, restores some of the details of the "true" 
annihilation energy distribution. This results in an increase in observed 
defect sensitivity. 

Even with the improvements in lineshape sensitivity and accuracy due 
to the FFT/PS algorithm, the OBPAS technique must have the additional ability 
to distinguish between multiple defect types if it is to become a useful tool 
in the study of industrial materials. In order to investigate the usefulness 
of deconvolution in resolving multiple defect types, a comparison was made 
between the values of the Rx parameters calculated before and after deconvolu- 
tion. If the number of reconstruction Fourier coefficients is allowed to vary 
between samples, as would be the case if the FFT/PS method is used in its 
standard form, then the calculated Rj parameter value yields little useful 
information. 

In order to explain this somewhat negative result, the consequences of 
truncating the Fourier coefficient, also known as low pass filtering, were ex- 
amined in greater detail. It is known that the truncation of the coefficients 
introduces a filtering error with an amplitude inversely proportional to the 
number of retained coefficients, w c , and a frequency which is proportional to 
w c . Ideally, one would like to greatly increase the number of retained coef- 
ficients in order to reduce the amplitude of this filtering error. Unfortu- 
nately, as w c is increased beyond the first minimum of the power spectrum the 
calculated solution will contain large, rapidly oscillating components which 
are not physically plausible. The mathematics of this argument were presented 
in detail in Chapter III. 

Since it is not possible to reduce the amplitude of the filtering 
error by arbitrarily increasing w c an alternative method was developed to 
improve the calculated values of. Rj for deconvoluted data. The first- step in -•" • 
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the procedure Is to calculate a w c value for each of the data sets to be 
analyzed. Next the minimum value of w c is selected and subsequently employed 
in the deconvolution of all the data sets. Using this procedure does not 
eliminate the filtering error, but rather controls the source such that its 
contribution to the calculated Rj values remains approximately constant from 
experiment to experiment. Although this solution is not ideal, the data 
presented in Chapter IV clearly indicate that by controlling the error source 
the reproducibility of the Rj parameter, and Its usefulness in the separation 
of multiple defect types, is improved. 

To summarize this section it can be stated that the FFT/PS deconvolu- 
tion algorithm appears to improve the quality of the data as obtained from 
UBPAS measurements in several ways. These improvements, as presented 1n this 
work include: 1) automatic correction for accidental and/or intentional 
changes in the instrument resolution function. 2) the potential of direct 
comparison of quantitative results between laboratories. 3) an improvement in 
both lineshape sensitivity and accuracy, and 4) the ability to distinguish 
between vacancy and dislocation type defects, as demonstrated for large 
grained high purity aluminum. 

To generalize these conclusions, further studies of this kind should 
be performed using a variety of different alloys with known defect type and 
concentration. In addition, it would be extremely useful to repeat the entire 
series of experiments using a OBPAS system which is both digitally and temper- 
ature stabilized. Such a data set would allow one to separate the effects of 
instrument broadening, which are theoretically correctable by deconvolution. 
from the nonsystematic effects of instrumental drift. In addition, spectra 
collected on a stabilized system could conceivably possess better counting 
statistics (due to lonyer allowable experimental ti^i.^ch \a 
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increase the number of coefficients retained in the deconvolution procedure. 

Additional experiments should be performed on the FFT/PS algorithm. 
These experiments could be designed to investigate alternative filtering 
schemes with the hope of reducing the associated filtering error. Most 
notably, the use of cosine tapering appears to be a viable approach. Once the 
filtering error has been reduced or eliminated the utility of moment analysis 
in the separation of multiple defect types should be reinvestigated. It 1s 
not known at this time if the current FFT/PS algorithm is capable of disting- 
uishing the formation of positronium in a nonmetallic material. Positronium 
would show Itself as a relatively high frequency characteristic of the annihi- 
lation energy distribution. The key factor Is determining whether the fre- 
quencies required to observe positronium are to the left of the first minimum 
in the power spectrum. This question could be investigated using both numeri- 
cally simulated data and data obtained from materials in which positronium 1s 
known to form. 

In order for Ooppler broadened positron annihilation spectroscopy to 
become a widely used industrial NOT tool it must continue to evolve in several 
ways. Its energy resolution must Improve, the ability to separate multiple 
defect types, must be further developed, and the reproducibility and reliabil- 
ity of the experimental results must increase. Although the fast Fourier 
transform/power spectrum analysis deconvolution technique has clearly not 
solved all of these problems, the results presented in this work indicate that 
this method measurably improves the quality and reproducibility of the infor- 
mation obtainable from DBPAS measurements, and serves as a sound basis for 
future refinements. 
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APPENDIX A 



In this appendix the equations used to theoretically determine the 
stopping or thermal izati on t times of positrons in metals will be discussed. 
For a more complete treatment of this subject see Brandt and Arista (1 982) • As 
briefly mentioned in the text, positrons slow down mainly by two mechanisms: 
1) ionization and excitation of electrons which is characterized by the elec- 
tronic stopping power, S e , of the material, and 2) phonon interactions charac- 
terized by the phonon stopping power, S p H* The second mechanism dominates at 
epithermal energies below some critical value E c characteristic of the mater- 
ial. In the range from the maximum source energy, E s , to E c the electronic 
stopping power dominates. This energy range can be further divided Into two 
subranges separated by the photon energy for which the material has its maximum 
stopping power, E m . In the range ^ < E < E s an accurate estimate for the 
stopping times is given by 

t! = 7.7p-l(Ei/mc2)1.2 (AJ) 

where t] is the stopping time in psec 

p is the density of the material in g/crn 3 
mc 2 is 0.511 MeV 
and £j is the initial positron energy 
since 

Ej/mc 2 < E s /mc 2 - I (A. 2) 

we can rewrite Eq. (A.l) as 

(104) 
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'< ' 7 -"-' (A.3) 

-hie, for most „.„ o, interest r esult s ,. stopp,„ a tines of - , psec ,„ ^ 
hiyh energy region. 

For E c < £ < E ro the electron gas theory of conduction electrons with 
Fermi energy E F is used to produce the equation 

S = -JL J!L \l]V2 E 5/2 

e 105 h l m J (A.4) 
in the case of a continuous slow down approximation 

t = - fJ!llV2 r E f dE 

s . z J J E i TTTTs (A.5) 
Substituting Eg. (A.4) into Eq. (A.5) and integrating yields 

t -iPl hE F 
2 ~ 8* TT (A. 6) 

fc f 

where t 2 is in picoseconds. 
In the ,ow enera, re*™ 3/2 KT < E < e c phonon interaction, d„„,.„ a te with 

0 3 2 
S „ = 2 m , W c 
PH "77 E (A.7) 

where w is the lattice deformation potential. 
Eq. (A.4) can be combined with Eq. (A.7) to define E c yielding 

.£ 3/2 . 105b 2 m 5 ' 2 E F 3 
C W 7 7~7~~ (A.8) 

Inserting Eq. (A.7) into Eq. (A.5) and integrating from 3/2 KT to E c as defined 
^ Eq. (A.8) yields 
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usiny typical values for common metals ft can be shown that 

l 3 > l 2 > *1 (A.10) 
Thus, phonon thermal izat ion times are larger than electronic stopping times. 
The total stopping time is given by 

t s = t } + t 2 ♦ t 3 (AJ1) 

with typical theoretical values of 13 psec in AC and 45 psec in Pb (Brandt and 
Arista 1982). 



APPENDIX B 



The purpose of this appendix is to describe the relationships between 
the momentum distribution of the electrons in the sample and the measured 
Doppler shifts in energy as well as the angular deviations from colinearity of 
the Y-rays. In the center of mass rest reference frame the annihilation Y-rays 
will be emitted in directions separated by 180° and each with an energy of 61 1 
KeV which is equal to the rest mass of an electron (or positron). However, in 
the laboratory frame the momentum of the electron (the positron is assumed to 
be at rest as described in the text) produces changes in both the observed y~ 
ray energy and the angle between the Y-rays. The geometry of the situation is 
shown in the diagram below, 

^ Center of Mass 
\ Velocity V 



\ 

\ 

\ 




Straightforward application of tf?e Lorentz transformations of special 
relativity combined with some trigonometric manipulation yield the result that 

7 2 
«♦ " e£) sin * o ♦ 8 (^) (8.1) 

where 6^ is the angular deviation from 180° 

v is the electron velocity 

(107) 
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c is the speed of light 

♦ 0 is the angle between the electron velocity vector and the 

direction of Y-ray emission 

y 2 2 
and is a term of order ^ which is negligible for v « c 

c c 

(For a detailed derivation of this result see Rindler 1960). The relativistic 
definition of momentum 1s 

P ■ <*n 0 v (B.2) 

where p is the momentum 

a mi - yr 1/2 

and m 0 is the particle mass 
For v « c, a ■*■ l and p - mv 
Substituting this result into Eq. (B.l) yields 

fi * " in^c* sin (B.3) 
Defining the transverse momentum component as 

PT * P sin ♦„ (B>4) 
and substituting Eq. (8.4) into Eq. (B.3) yields the final result 

2P T 

6<l " STc (B.5) 

O ' 

demonstrates that measuring the distribution of the angular deviations. 6«, is 
in essence equivalent to measuring the distribution of the momentum of the 
electrons within the sample. 

In the case of Doppler measurements one has 

6P Y = P e cos * 0 (B>6) 
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where 6P y is the change in momentum of the Y-ray 

P e is the electron momentum 
and <(. 0 is as defined above 
Using Eq. (B.2) with v « c Eq. (B.6) takes the form 

6P Y = mev e cos * 0 {Bm7) 
It can be shown that in relativistic dynamics 



(B.8) 



• E 2 = P 2 c2 + moc2 
where E is the relativistic energy 

P is the relativistic momentum 

c is the speed of light 
and m 0 is the particle rest mass 

Now for a photon (y-ray) m 0 = 0 so Eq. (B.8) reduces to 

E Y = Py c (B.9) 
Combining Eqs. (B.6), (B.7), and (B.9) yields 

6E Y = 6 PyC = p e (cos * 0 )c = mgcvetcos * 0 ) ( B .l 0 ) 
where 6E Y is the Doppler energy shift of the observed Y-ray 
6p Y is the change in momentum of the Y-ray 
and p e is the momentum of the annihilating electron 

From Eq. (B.lO) 1t can be seen that the momentum distribution of the annihi- 
lating electrons can be deduced from the measured Doppler energy shifts. 

Combining Eq. (B.5) with (B.lO) yields the following relationship 
between 64. and 6E: 

$E = 255. 56$ I ^ j 

where 6E is in eV 
and 6(j> is in mrad 
Finally, by combining Eq. (B.lO) and the fact that for an electron m 0 c2 = 51 1 
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KeV it can be shown that 

6E T = 511 (cos ♦ 0 )(i) (B.12) 

where 6E Y is in KeV 

For the case of a 1 eV electron (v e = 5.93 x 10? cm/sec) moving in the direc- 
tion of Y-rajr emission, £q. B.12 predicts a 1.01 KeV shift in observed energy. 



APPENDIX C 



The material examined in this work was high purity (99.9995%) aluml 
supplied by the Materials Research Corporation. The chemical analysis, in 



r ppm of 


metallic 


interstitials 


(and some gases) is 


listed bel 


Element 


PPM 


Element 


PPM 


Element 


PPM 


Ag 


<0.1 


L1 


<0.1 


U 




A 1 

A1 


major 


Mg 


1.5 


V 


<0.1 


As 


<0.1 


Mn 


<0.1 


W 


<0.2 


Au 


<0.1 


Mo 


<U.l 


y 




B 




Na 




Zn 


0.4 


U r» 

Dd 




Nb 


<0.1 


Zr 


<0.1 


Be 




N1 


U.l 






Bi 


<0.1 


Pb 


<0.1 


CI 


0.6 


Ca 


<0.1 


Pd 


<oa 


F 




Cd 


<0.1 


Pt 


<0.14 


H 


3 


Co 


<0.1 


Hh 


<0.1 


N 


5 


Cr 


<0.1 


Sb 


<0.1 


0 


20 


Cu 


0.1 


S1 


0.9 


P 


0.3 


Fc 


<0.1 


Sn 


<0.3 


S 




Ga 


<0.1 


Ta 








Ge 


<0.1 


Th 








In 


<0.1 


Ti 


<(U 






K 


<U.l 


Tl 









ow. 



(Ill) 



APPENDIX D 



In this Appendix the utility of Fourier analysis in the solution of the 
convolution integral equation will be demonstrated, m the first section of 
Chapter III it was shown that 

h(y) = f( x)r(y . x)dx ([) i} 

By definition of the Fourier transform 

F(w) = f eiwx f(x)dx (d z) 

Therefore 

F(«)K(w) = eiwx f(x)dx]r ^ e iwt r(t)dt] 

co ay 

= [ m [ m e iwx f(x)eiwt r(t)dxdt (D>4) 
S /. L eiW(X+t,f (x)r(t)dxdt (D . 5) 

Now let x + t = y 
so t = y - x 
and dt = dy 
Hence Eg. (D.5) can be rewritten as 

F(w)R(w) = f f eiwy f(x)r(y . x)(Jxd> (Df>) 

CO oo 

- J m elwyr/^ f(x)r(y . x)(Jx]dy (q ?) 

(112) 



113 

But the Integral within the brackets is by Eq. (U.l) equal to h(y). Therefore 
Eq. (U.7) can be rewritten as 

F(w)R(w) = / ei w yh(y)dy (0.8) 

and using Ecj. (D.2) on the right hand side 

F(w)R(w) = H(w) (D.9) 
This demonstrates that in Fourier space convolution is simply a complex multi- 
plication. 



APPENDIX E 



The purpose of this Appendix is to demonstrate some of the mathematical 
difficulties associated with the incorrectly posed convolution integral 
equation. The primary source of the problem is that the left hand side of the 
equation 

^ b 

n U) * / f(y)r(x - y)dy ( E .n 

is only known to within some error e. That is, if h(x) is the "true" function 
then 

max jh*(x) - h(x)| <■ c (E ?) 

a<x<b 

Any function f^y) which when convoluted with the known resolution function 
r(x) yields a function h^x) that satisfies Eq. (E.2) must be considered a 
solution of Eq. (E.l). 

Assume there are two solutions to Eq. (E.l) which are f^x) and f 2 (x) = 
f](x) + c sin(2imw/(b-a). By definition 

* b 

"l(x) = / fj(y)r(x - y)dy (£.3) 

and 

* b 

M*) 3 / Cfi(y) + C sin (2*ny/(b - a))]r(x - y)dy (£.4) 

d 

Subtracting Eq. (E.3) from Eq. (E.4) yields 
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* * b 

h 2 (x) - h^x) = J r(x - y)C sin(2»ny/(b - a))dy (E.5) 
a 

Now inteyrating by parts with 

u = r(x - y) dv = C sin (f^Jdy 

yields 

/* r(x - y)C sln^ldy = r(x - y>[^i|-^ «- <g*ff>< 

Now integrating the last term on the right hand side by parts with 
du = . { 2l^_) sin (jaiL, dy v . r(x . y) 

Yields 

It <~$*i) d(r V y " * ■ Cr(x - y) cos (ffiLjjJ 

- £ - r(x - y)<|^> stn (fSSL, dy (E.7) 

Combining Eqs. (E.6) and (E.7) yields 

/J f(x . y) C sin (fffllL, dx ») |>(x - y) cos (E.8) 

and combining Eqs. (E.8) and (E.£>) yields the result 

h*(x) - h*(x) *itfy=r±) r(x - y) cos (f=*y)]J (E.9) 



So it can be seen that for any value of C one can select a value of N large 
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enough so that the difference between h^x) and h^x) is smaller than any 
previously determined number. Therefore, both h|(x) and h^x) satisfy Eq. 
(£.2) and f,(x) as well as f 2 (x) must be solutions to Eq. (E.l). The important 
point is that one may select C and N in such a way as to allow f,(x) and f 2 (x) 
to differ from each other by as much as one likes. Clearly the problem as 
posed does not have a unique solution and a different solution technique is 
required. 



1 



APPENDIX F 



ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

C THIS IS A LISTING OF THE FORTRAN 77 VERSION OF THE PROGRAM C 

C DECONFFT WHICH IS DESIGNED TO DECONVOLUTE DOPPLER BROADENED C 

C POSITRON ANNIHILATION SPECTRA. FOR A DESCRIPTION OF TBE C 

C MATHEMATICS INVOLVED SEE J. P. SCHAFFER, E.J. SHAOGHNESSY, C 

C AND P.L. JONES, JOURNAL OP NUCLEAR INSTRUMENTS AND METHODS C 

C IN PHYSICS RESEARCH, B5, 74, 1964. C 

C C 
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 



INTEGER PCSP 
CHARACTER* 2 PCSZ,PCB 
CHARACTER* 80 XCH,YCH,TCH 



CHARACTER* 7 NAME, NAME E 

DIMENSION I A (2 04 8), IB(2048), DOPE (12) ,XP(74) ,YP(74) 
DOUBLE PRECISION A(2048) , B(2048) ,C(2048) ,FR(2048) ,HC(2048) 
DOUBLE PRECISION D(2048) , E (2048) ,G (2) ,S (2048) ,S1 (2048) 
DOUBLE PRECISION Fl (2048) ,PI, V,W,W1 
COMPLEX*16 R(4096),H{4096) ,F(4096) ,T,A1,H1(4096) 
COMPLEX*16 A2,R1(4096) ,T1 
PI=3.141592653589793D0 

PRINT* , 'THIS PROGRAM DEALS WITH THE CONVOLUTION" 
PRINT*, »H(X)»F(X) *R(X) , WHERE R IS THE* 
PRINT* , * RESOLUTION FUNCTION, H IS THE' 
PRINT* , ' BROADENED AND OBSERVED FUNCTION, AND* 
PRINT* , 0 F IS THE TRUE AND DESIRED FUNCTION.' 
PRINT* , * ' 

C THE RESOLUTION FUNCTION IS USUALLY OBTAINED FROM 
C THE Sr(8S) GAMMA RAY (51< KeV) . 



C 

2 PRINT* , ' THE PROGRAM REQUIRES N DATA POINTS* 
PRINT*, 'WHERE N*2**M, AND M IS AK INTEGER.' 
PRINT*,' ' 

C THESE DATA POINTS WILL BE H(X) AND R(X) AS COLLECTED 
C BY THE MCA. 

PRINT*, 'ENTER THE VALUE OP M TO BE USED 
PRINT*, 'PROGRAM. (M=ll YIELDS N=2048) ' 



c 
c 



IN THIS' 



PRINT*, ' 
READ(*, ' (12) ' ) M 
N=2**M 
ND2=N/2 
ND4-N/4 

PRINT*, 'ENTER THE RESOLUTION FUNCTION FILENAME' 

PRINT*,' ' 

READ (*, ' (A7) ') NAME 

OPEN (90,FILE=HAME,PAD='yES') 

READ (90,1) (IA(I), 1=1,2048) 

FORMAT - (5X,.5I7 ) 

CLOSE (90) 
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SxS?;! 8 "™.™ 8 EXPEMMENTAL FUNCTION FILENAME 1 

R£AD(*|'(A7) ') NAME 

OPEN (90, PILE= NAME, PAD c • YES' ) 

READ (90,1) (IB(I), 1=1,2048) 

CLOSE (90) r*»«a/ 

NIN-2048 

n ™* SUBROUTINE DELSKIP DELETES THE LINE 
C SKIP IN THE RAW DATA SUPPLIED BY THE MCA 
CALL DELSKIP (IA,A,NIN) 
CALL DELSKIP (IB,B,NIN) 



C 
C 
C 
C 



C 



THE FOLLOWING SUBROUTINE CALCULATES THE CENT PAT 
f? A !SSJf CATION BASED °» THE c¥otI R of 5flf 
11 CHANNEL REGION WHICH CONTAINS THE MAXIMUM 
TOTAL NUMBER OF COUNTS ™* iMUM 



CALL CENTER (A, IMAXA) 
CALL CENTER (B,IMAXB) 



C SSIS'Mr DEALS „™ THE 

Km*! W)"^ 10 correct? • 

IF (NAME -EQ. 'y«) GO TO 3 
CALL NOBACK (N,B,IMAXB) 
GO TO 200 

3 PRINT*,* DO YOU WISH TO BACKGROUND CORRECT dctht" 

SaSl-'SBrigSS",^".- ^TRAL FILTERING • 

READ (*,'(A1) ') NAME 

IF (NAME .EQ. »J') GO TO 5 

PRINT*, 'THE NEXT SEVERAL QUESTIONS REFER TO tbpi 
PRINT*, 'RESOLUTION FUNCTION. • THE 
CALL BACK (A, IMAXA, N) 

PRINT*, 'THE NEXT SEVERAL QUESTIONS REFER TO TUP- 
PRINT*, 'EXPERIMENTAL FUNCTION.' 
CALL BACK (B, IMAXB, N) 
GO TO 377 

5 ocSI r ! E 2?P TSE RESOLUTION OF THE MCA IN' 
PRINT* , eV/ CHANNEL. (USE A DECIMAL POINT)' 
PRINT*, '(SUGGESTION: IF THE SYSTEM^ IS SET 

pS*; : IN " s N0W1AL M0DE °" ™- ' 

READ (*,'(F6.2) ') EV 

X«(1000.D0/EV) +0.5D0 

IK=IFIX (X) 

IKLA= ( IKAXA-2 0 ) / IK 

ILA=514-IKLA 

IKHA= (1024- IMAXA) /IK 

IHA=514+IKHA 

lKLB=(IMAXB-20)/IK 
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ILB-511-IKLB 

IKHB= (1024-IMAXB)/IK 

IHB-511+IKBB 

PRINT* , ' THE EXPERIMENTAL SPECTRUM HAS BEEN COLLECTED ' 
PRINT* , * OVER THE FOLLOWING KeV REGION:' 
WRITE (*,10) ILB,IHB 
10 PORMAT (15,15) 

PRINT*, 'WE ARE NOW READY TO PROCEED WITH THE' 
PRINT* BACKGROUND CORRECTION. OVER WHICB' 
PRINT*,' LOW ENERGY REGION WOULD YOU LIKE THE' 
PRINT*, 'LINEAR LEAST SQUARE FIT? IF POSSIBLE' 
PRINT*, 'SELECT 486 TO 497. THE PROPER FORM' 
PRINT*,' OF ENTRY IS TWO 3 DIGIT INTEGERS' 
PRINT*, 'SEPARATED BY TWO BLANKS. * 
PRINT*,' ' 
READ (*,*) J LAB, J LB B 
IS-IMAXB- (511-JLAB) *IK 
IF-IMAXB- (511-JLBB) *IK 
N1=IF-IS 
DO 20 I«1,N1 
J=IS-1+I 

D(I)=J ) 

E(I)=B(J) 
20 CONTINUE 

CALL FIT (N1,D,E,G) 

BBL=G(1) 

XMBL=G(2) 

PRINT*, 'OVER WHICH HIGH ENERGY REGION WOULD' 

PRINT*, 'YOU LIKE THE LEAST SQUARE FIT?' 

PRINT*, 'IF POSSIBLE SELECT 534 542.' 

PRINT* , ' ' 

READ (*,*) JHAB,JHBB 

IS=IMAXB+( JHAB-511) *IK 

IF=IMAXB+( JHBB-511) *IK 

N1*IF-IS 

DO 30 1=1, Nl 

J=IS-1+I 

D(I)=J 

£(I)«B(J) 
30 CONTINUE 

CALL FIT (Nl,D,E,G) 

BBH«=G(1) 

XMBB=G(2) 

PRINT*,' TO WHICH ENERGIES WOULD YOU LIKE TO' 
PRINT* , 1 EXTRAPOLATE TBE LINEAR BACKGROUND' 
PRINT* , ' CORRECTION? (SUGGESTION: 503 518)' 
PRINT* , ' ' 
READ (*,*) IL,IH 

PRINT* , ' ON THE HIGH ENERGY SIDE OF THE PEAK, ' 
PRINT*,'WHAT IS THE HIGHEST ENERGY YOU WISH' 
PRINT*,' TO HAVE CONSIDERED FOR DECONVOLUTION? ' 
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PRINT*, 'SUGGESTION: 540' 

PRINT* , ' 1 

READ (*,*) I END 

IL=IMAXB- (511-IL) *IK 

IB-IMAXB+(IH-511)*IK 

IEND=IMAXB+ ( IEND-511) *IK 

IEB-IEND 

DO 35 I=IH,IEND 

B ( I ) -B ( I ) - ( XMBH*I+BBH) 

IP (B(I) .LT. 0.) B(I)«0.D0 
35 CONTINDE 

XH-XMBH*IH+BBH 

XL-XMBL*IL+BBL 

XM= ( XL-XH ) / ( FLOAT ( IL- IB ) ) 

XB=XH-( ( XL- XH)/( FLOAT (IL-IH)) ) *IH 

DO 40 I=IHAXB, IH 

B(I)=B(I)-(XM*I+XB) 

IP (B(I) .LT. 0.) B(I)*0.D0 
40 CONTINUE 

PRINT* r ' THE RESOLUTION SPECTRUM HAS BEEN COLLECTED' 
PRINT*, 'OVER THE FOLLOWING KeV REGION: ' 
WRITE (*,10) ILA, IHA 

PRINT*,' OVER WHICH LOW ENERGY REGION WOULD TOD LIKE ' 

PRINT*,' THE LINEAR LEAST SQUARE PIT?' 

PRINT*, 'IP POSSIBLE SELECT 489 500.' 

PRINT*,' ' 

READ (*,*) JLAA, JLBA 

IS=IMAXA- (5 14- JLAA) *IK 

IF=IHAXA- (5 14 -JLBA) *IK 

N1=IF-IS 

DO 45 I«1,N1 

J=IS-1+I 

D(I)=J 

E(I)=A(J) 
45 CONTINUE 

CALL FIT (N1,D,E,G) 

BAL=G(1) 

XMAL=G(2) 

PRINT* * 1 OVER WHICB HIGH ENERGY REGION WOULD YOU LIKE ' 

PRINT*, 'THE LINEAR LEAST SQUARE FIT?' 

PRINT*,'IP POSSIBLE SELECT 537 545.' 

PRINT*,' ' 

READ (*,*) JHAA, JHBA 

IS»IMAXA+(JHAA-514) *IK 

IF=IMAXA+(JHBA-514) *IK 

N1=IP-IS 

DO 50 1=1, Nl 

J=IS-1+I 

D(I)=J 

B(I)=A(J) 
50 CONTINUE 



CALL PIT <N1,D,E,G) 

BAH«G(1) 

XMAH=G(2) 

PRINT*, 'TO WHICH ENERGIES WOULD YOO LIKE TO' 
PRINT*, 'EXTRAPOLATE THE LINEAR BACKGROUND ' 
PRINT* , • CORRECTION? (SUGGESTION: 506 521)' 
PRINT*, • • 
READ (*,*) IL, IB 

PRINT*, 'ON THE HIGH ENERGY SIDE OF THE PEAK, 1 
PRINT*, 'WHAT IS THE HIGHEST ENERGY YOO WISH' 
PRINT*, 'TO HAVE CONSIDERED POR DECONVOLUTION" 
PRINT*, 'SUGGESTION: 521' " ' 

PRINT*,' ' 
READ (*,*) I END 
IL-IMAfA- (514-IL) *IK 
IH=IMAiA+(IH-514)*IK ' 
IEND=IMAXA+(IEND-514) *IK 
I EA= I END 

DO 55 I=IH,IEND 

A ( I) «A ( I ) - ( XMAH*I+BBH) 

IP (A(I) .LT. 0.) A(I)=0.D0 
55 CONTINUE 

XB=XMAH*IH+BAH 

XL=XMAL*IL+BAL 

X«=(XL-XH)/(FLOAT(IL-IH) ) 

XB*XB- ( (XL-XH) / ( FLOAT ( IL- IH) ) ) * IH 

DO 60 I=IMAXA,IH 

A(I)=A(I)-(XM*I+XB) 

IF (A(I) .LT. 0.) A(I)=0.D0 
60 CONTINUE 
GO TO 378 

377 CONTINUE 
IEA-IMAXA+109 
IEB=IMAXB+453 

378 CONTINUE 

C 

C THIS IS THE END OP THE BACKGROUND CORRECTION SEQUENCE 
C NEXT WE NEED TO REPLECT THE HIGH ENERGY HALF OF THE " 
C BACKGROUND CORRECTED PEAK ABOUT THE CENTRAL POINT 
C TO FORM A SYMMETRIC PEAK. IN ADDITION A ZEROFILL 
C TO 2048 DATA POINTS IS REQUIRED 



C 



CALL 2EROFILL (A, IMAXA, IEA,N) 

CALL ZEROFILL (B,IMAXB,IEB,N) 

M1=IEA-IMAXA+10 

CALL TAPER (A,M1,N) 

M2=IEB-IMAXB+20 

CALL TAPER (B,M2,N) 

M3=K1 

CALL SMOOTH (A,M3,N) 
M4=M2 



CALL SMOOTH (B,M4,N) 
C THE NEXT TWO DO LOOPS NORMALIZE THE AREA 
C OP TEE RESOLUTION AND BROADENED SPECTRA. 
V=0.D0 
W-O.DO 
DO 62 1=1, N 
V-V+A(I) 
W=W+B(I) 

62 CONTINUE 
DO 63 I-1,N 

A(I)=A(I)/V 
B(I)=B(I)/W 

63 CONTINUE 
DO 68 1-1, N 

S(I)=B(I) 
68 CONTINUE 

DO 72 1=1, ND2 

J-ND2+I 

V=S(I) 

S(I)=S(J) 

S(J)=V 
72 CONTINUE 

PRINT*, •DO YOU WISH TO SAVE THE BACKGROUND ' 
PRINT*, 'CORRECTED SPECTRUM? 1 
READ (*,' (Al) •) NAME 
IF (NAME .EQ. 'N* ) GO TO 76 
PRINT*, 'ENTER THE FILENAME. ' 
READ (*,' (A7) ') NAME 
OPEN (90 , FILE-NAME , PAD= ' YES * ) 
WRITE (90,74) (S(I),I=1,N) 
74 FORMAT (5X,5D22.10) 
CLOSE (90) 
CALL PARAMS (N,S) 

76 CONTINUE 

DO 77 1-1, N 

R(I)=DCMPLX(A(I) ,0-DO) 
H(I)=DCMPLX(B(I) ,0.D0) 

77 CONTINUE 

CALL REORDER (M,H) 
CALL FFT (M,H) 
CALL REORDER (M,R) 
CALL FFT (M,R) 
CALL DECON (M,F,R,H) 

DOPE(1)=150. 
DOPE (2) =150. 
DOPE{3)=950. 
DOPE (4) =650. 
DOPE (5) =0. 
DOPE (6) =0 . 
PCSZ='NS' 
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PCSP-1 

I ^TS? ^ 0NE T ° ™ *» POWER 
READ (*,' (Al) ') NAME 

£ iTZiA ,N,) 60 T0 179 

78 CON?^ H(I) * (CONJG(H(I))) 

p££*;. ,EOT ? R ™ FILENAME P0R ^ HE *.S. OP THE EX. DATA. • 

R£AD (*,'(A7) •) NAME 

OPEN ( 90 , FILEs NAME , P AD= • YES 1 ) 

am ~ 1TE < 90 ' 90 D (I#D(I) # I=1 # ND2) 
901 FORMAT (5(14,022.12)) 

CLOSE (90) 

CALL REORDER (M,F) 

CALL FFT (M,F) 

TaO.DO 

DO 277 1=1, N 
T=T+F(I) 

277 CONTINUE 

DO 278 1=1, N 
F(I)=P(I)/T 

278 CONTINUE 

CALL REORDER (M,F) 
CALL FFT (M,F) 
DO 178 1=1, ND2 

ITS CoS P(1) * (CO,UG(P(I '" 

p5S?.*:.' E,fT f E T8E PILENAME P0R ™ »••• OP THE DECON. DATA. 1 

READ <*,'<A7)') name 

OPEN (90, FILE- NAME, PAD-' YTS") 

Son Ha,' 901 ' «'«».ii.«H 

NPTS«74 

DO 910 1=1,10 

XP ( I) =ALOG10 ( FLOAT ( I) ) 

YP(I)=E(I) 

YP(I)=ALOG10(YP(I)) 
910 CONTINUE 

DO 912 1=15,100,5 
J=8+(I/5) 

XP ( J) =ALOGl 0 ( FLOAT (I) ) 
YP(J)=E(I) 

VP(J)=ALOG10(YP(J)) 
912 CONTIN0E 
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DO 914 1-120,1000,20 
J-23-MI/20) 

XP { J) = ALOG1 0 ( PLOAT ( I ) ) 

YP(J)=E(I) 

YP(J)=ALOG10(YP(J)) 
914 CONTINUE 

XP(74)«ALOG10(1020.) 
YP(74)«-10. 

XCH»'LOG OP THE FOURIER COEFFICIENT NO. ' 
MX- 3 4 

YCH='LOG OP THE POWER' 
MY=16 

PRINT*, 'ENTER THE TITLE OP THIS PLOT. 1 

PRINT*, ' ' 

READ (*,'(A80) ') TCH 

PRINT*, 'ENTER THE NO. OF CHARACTERS IN THE TITLE. ' 
READ (*,' (12) ') NT 
PCH»*1 1 

CALL CKSTDP ( 'TK15 ' ,9600) 
CALL CKER 

CALL CKGRAF (DOPE,NPTS,XP, YP,XCH,MX,YCH, MY,TCH,NT, PCSZ 
C,PCSP,PCH) 
DO 930 1=1,10 
YP(I)=D(I) 
YP(I)=ALOG10(YP(I)) 
930 CONTINUE 

DO 932 1=15,100,5 
J=8+(l/5) 
YP(J)=D(I) 
YP(J)=ALOG10(YP(J)) 
932 CONTINUE 

DO 934 I«120, 1020,20 
J=23+(I/20) 
YP(J)-D(I) 
YP(J) =ALOG10(YP(J)) 
934 CONTINUE 
PCB-'3 ' 

CALL CKPLT (DOPE,NPTS,XP, YP,PCSZ ,PCSP,PCB) 
CALL CKMA (0. ,0. ,1) 
READ (*,*) 
179 CONTINUE 

C 

C AT THIS STAGE THE RESOLUTION AND EXPERIMENTAL 
C SPECTRA HAVE BEEN TRANSFORMED INTO FOURIER SPACE. 
C THIS NEXT SMALL SECTION DETERMINES THE NUMBER OF 
C POSITIVE DECREASING FOURIER COEFFICIENTS IN THE 
C DECONVOLVED RESULT. THIS INFORMATION CAN THEN 
C BE USED IN THE FILTERING ROUTINES IF SO DESIRED. 
C 

DO 80 1=2, N 
J=I-1 
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D(I)=REAL(P(I)) 
D(J) =REAL(P(J)) 
IP (D(I) .LT. 0.0) GO TO 81 
_ _ IF (D(I) .GT. D(J) ) GO TO 81 

80 CONTINUE 

81 IPD-I-1 

WRITE (*,82) IPD 

82 pS?!S* r ' (5X,,TH ? N °' ° F P0SITIVE DECREASING COEPS. IS • 15) 
PRINT*, 'DO YOO WISH TO OSE A COSINE TAPER FILTER' 

READ (*, ' (Al) •) NAME 

IF (NAME .EQ. 'N') GO TO 189 

PRINT*. 'ENTER THE NUMBER OP THE FIRST COEFFICIENT ' 
PRINT*, 'TO BE SET EQOAL TO ZERO. • 
READ (*»'(I3) ') IZ 

R^JS T (*, , ?uff.) ,H fT NDMBER ° P C0EPFICIENTS TO BE TAPERED. • 
IS-IZ-IT 
XS^REAL(IS) 
XZ-REAL(IZ) 
XY=1.0/(X£-XZ) 
DO 181 I=IS,IZ 
X=REAL(I) 

,o, ~ F(I )= p (D*0.5*(l-COS(PI*XY*(X-XZ))) 
181 CONTINOE 

J=IZ+1 
K=(N/2)+l 
DO 183 I«=J,K 

F(I)=DCMPLX (O.DO,0.DO) 
183 CONTINUE 
K=K— 1 

DO 185 1=2, K 
J-N+2-I 
P(J)=P(I) 
185 CONTINUE 

189 PRINT*, 'DO YOU WISH TO USE A BOXCAR FILTER TO* 
PRINT*, 'MODIFY THE FOURIER COEFFICIENTS IN' 

S^™-? 0 ^ 0 ™ S1SB FRE0 - »° 1SE? • 

IF (NAME .EQ. 'N') GO TO 84 

PRINT*,'ENTER THE NUMBER OF COEFFICIENTS TO BE USED: • 

J=J+1 

K=N+2-J 

DO 83 I=J,K 

F(I)=DCMPLX(0.D0,0.D0) 

83 CONTINUE 

84 CONTINUE 

reaS T (*!'?a™) W Ze T0 SEE THE FO0RIER coeffi cients?' 



IP (NAME . EQ. 'N' ) GO TO 85 
PRINT* , ' ENTER THE FILENAME FOR H(X). ' 
PRINT* , ' ' 
READ (*,' (A7) ') NAME 
OPEN (90, FILE-NAME, PAD=' YES') 
WRITE (90,911) (I,B(I), 1-1, H) 
911 FORMAT (2(I4,2D22.10)) 
CLOSE (90) 

OPEN (gCFILE-'FCSF^PAD-'YES') 
WRITE (90,911) (I,P(I), I«1,N) 
CLOSE (90) 

85 CONTINUE 

CALL REORDER (M,F) 

CALL FFT (M,F) 

CALL REORDER (M,H) 

CALL FFT (M,H) 
C THIS NEXT SECTION ALLOWS ONE TO USE THE ITERATION 
C ROUTINE TO IMPROVE F(X). 

PRINT*, 'WOULD YOO LIKE TO USE THE ITERATION* 

PRINT*, 'SUBROUTINE TO TRY TO IMPROVE THE' 

PRINT*, 'ACCURACY OF THE DECONVOLUTED SPECTRUM? 

READ (*,' (Al) ') NAME 

IF (NAME .EQ. 'N') GO TO 86 

CALL ITER (M,R,F,H) 

86 CONTINUE 

C THIS NEXT DO LOOP SHIFTS THE PHASE OF THE 
C EXPERIMENTAL FUNCTION B, AND THE DECONVOLUED 
C FUNCTION F, SO THAT THEIR MAXIMA ARE LOCATED 
C IN THE CENTRAL CHANNELS. 
ND2=N/2 
DO 88 1=1, ND2 
J-ND2+I 
T=F(I) 
Tl«B(I) 
F(I)»F(J) 
B(I)=B(J) 
F(J)«T 
B(J)=T1 
88 CONTINUE 
DO 98 1=1, N 

Fl(I)=R£AL(F(I) ) 

98 CONTINUE 
W-0.D0 

DO 99 1=1, N 
W-W+FKI) 

99 CONTINUE 

DO 100 1=1, N 

F1(I)=F1(I)/W 
100 CONTINUE 

PRINT* , 1 DO YOU WANT A PLOT OF h(X) AND f(X)? ' 
READ (*, 1 (Al) ') NAME 



IP (NAME .EQ. 'N* ) GO TO 105 
NPTS-60 
W=P1(1025) 
DO 102 1=1,60 
J=1024+I 

XP(I)=.064*(FLOAT(I)) 
YP(I)=P1(J)/W 
102 CONTINUE 

XCH-'KeV 
MX«=3 

YCB-'MX) AND f (X) * 



IN THE TITLE. 1 



MY«13 

PRINT* '! EM ? ER TBB TITLB 0F THIS PLOT. ' 
READ (*,• (A80) •) TCH 
PRINT* , * ENTER THE NO. OP CHARACTERS 
READ <*,'(I2)») NT 
PCH-'l • 

CALL CKSTDP ( 'TK15* ,9600) 
CALL CKER 

C^SP??^ (D0PE,HPTS ' 3(P '»rXCH,MX.YOT,MY,TCH,NT,PCSZ 
DO 103 1=1,60 



J=1024+I 

TP(I)=S(J)/W 
103 CONTINUE 



105 



PCH='3 ■ 

^H?:r s '"' YP ' PC "' i ' CSI, ' I,CB) 

READ (*,*) 

PRINT*,* DO YOU WISH TO SAVE THE • 
PRINT* , * DECONVOLUTED SPECTRUM? ' 
READ (*,' (Al) •) NAME 
IP (NAME .EQ. 'N' ) GO TO 101 
PRINT*, 'ENTER THE FILENAME. ' 
READ (*,'(A7)«) NAME 
OPEN (90, FILE=NAME , PAD= * YES ' ) 
WRITE (90,74) (P1(I) ,I=1,N) 
CLOSE (90) 



CLOSE (90) 
101 CONTINUE 



UJWXiWUt 

PRINT*, 'DO YOU WISH TO CHECK YOUR RESULT BY* 

l^2E!'. ,C0NV0LDTING f(X > WITH r < X > AND COMPARING 1 
PRINT*, 'TBE RESULT WITH h(X)?' 

READ (*, , (A1)«) NAME 



IP (NAME . EQ. »N') GO TO 200 
CALL REORDER (M,P) 
CALL PFT (M,F) 
CALL CONV (M,F,R,H1) 
CALL REORDER (M,H1) 
CALL PFT (M,H1) 
DO 106 1=1, N 
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HC(I)-R£AL(H1(I)) 

106 CONTINUE 
V*0.D0 

DO 107 1=1, N 
V=V+HC(I) 

107 CONTINUE 

DO 108 I«1,N 

HC(I)«=HC(I)/V 

108 CONTINUE 

PRINT* , 1 ENTER THE FILENAME FOR THE RESULT ' 

PRINT*,' OP THIS COMPARISON. ' 

PRINT* , ' ' 

READ (*, ' (A7) ') NAME 

OPEN (90, FILE=NAME,PAD= * YES ' ) 

WRITE (90,110) 

110 FORMAT (3X,»lM4X,'h(X)RAW',14X,'h(X) = {F*R) \14X,'f (X) ') 
WRITE (90,120) (I, S(I),BC(I),F1(I), 1=1025, 1300) 

120 PORMAT (1X,I4,2X,3D22.10) 
CLOSE (90) 

PRINT*,' DO YOO WANT A PLOT OF h(X) AND THE ESTIMATE' 
PRINT*, 'OF h(X) OBTAINED FROM f(X)*r(X)? • 
READ (*,' (Al) ') NAME 
IF (NAME . EQ. 'N') GO TO 200 
NPTS=60 
W=HC (1025) 
DO 130 1=1,60 
J=1024+I 

XP(I)=.064*(FLOAT(I)) 
YP(I)«BC(J)/W 
130 CONTINUE 
XCB='KeV 
MX=3 

YCB='h(X) AND h*(X)' 
MY-14 

PRINT* ENTER THE TITLE OF THIS PLOT.' 

PRINT*,' ' 

READ (*,' (A80) ') TCH 

PRINT* ,' ENTER THE NO. OF CHARACTERS IN THE TITLE. ' 
READ (*, ■ (12) ') NT 
PCB«'l ' 

CALL CKSTUP ( 'TK15' ,9600/ 
CALL CKER 

CALL CKGRAF (DOPE,NPTS,XP, YP,XCH,MX, YCH, MY,TCH,NT,PCSZ 
C,PCSP,PCH) 
DO 140 1-1,60 
J=1024+I 
YP(I)=S(J)/W 
140 CONTINUE 
PCH='3 ' 

CALL CKPLT ( DOPE, NPTS, XP , YP, PCS Z , PCS P, PCH) 
CALL CKMA (0. ,0. ,1) 
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READ (*,*) 
200 CONTINUE 

CALL P ARAMS (N r Pl) 

STOP 

END 

cccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE DELSKIP (IC,C,NIN) 

DIMENSION IC(2048) 

DOUBLE PRECISION C(1024) 

J=0 

K«=0 

DO 10 I«1,NIN 
J=J+1 
K-K+l 

IC(I)«IC(K) 

IP (J . EQ. 5) THEN 

J=0 

K»=K+5 

END IF 
10 CONTINUE 
NID2=NIN/2 
DO 20 I=1,NID2 

C(I)=FLOAT(IC(I)) 
20 CONTINUE 
RETURN 
END 

C 

SUBROUTINE CENTER (A,IMAX) 

DOUBLE PRECISION A (1024) ,X,T(21) ,Y 

X=0.D0 

DO 10 1*300,800 

IP (A(I) .GT. X) THEN 
X»A(I) 
J«I-10 
ENDIP 
10 CONTINUE 
K«=J+21 
Y=0.D0 
DO 20 I=J,K 
L»I-5 
M=I+5 
T(I)=0.D0 
DO 15 N»L,M 

T(I)=T(I)+A(N) 
15 CONTINUE 

IF (T(I) .GT. Y) THEN 
Y=T(I) 
IMAX=I 
END IF 
20 CONTINUE 
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RETURN 
END 

C 

SUBROUTINE PIT (N,X,Y,C) 
C T2IS SUBROUTINE IS A SIMPLIFICATION OF 
C A LARGER SUBROUTINE WBICH FITS A CURVE USING 
C UP TO A NINTH ORDER POLYNOMIAL. FOR A 
C REFERENCE SEE JPS DISSERTATION. 

DOUBLE PRECISION X{300) ,Y{300) ,C(2) ,A(3,3) ,XN(300) 
DO 10 I«1,N 
10 XN(I)=1.D0 
DO 30 1-1,2 
A(I,1)=0.D0 
A(I,3)=0.DD 
DO 20 J=1,N 

A(I,1)=A(I,1)+XN(J) 
A(I,3)=A(I,3)+Y(J)*XN(J) 
20 XN(J)=XN(J)*X(J) 
30 CONTINUE 
A(2,2)=0.D0 
DO 40 J=l ,N 

A(2,2)=A(2,2)+XN(J) 
40 XN(J)=XN(J)*X(J) 
A(1,2)=A(2,1) 
DO 50 1-1,2 

SUM=0.D0 

IF (2 .GT. I) GO TO 45 
SUM-SUM+A ( 1 , 1) *A (1 , 2 ) 
A(I,2)«A(I,2)-SUM 
GO TO 50 

45 IF (ABS(A(I,I)) .LT. l.E-10) GO TO 99 

A(I,2)=(A(I,2)-SUM)/A(I,I) 
50 CONTINUE 

DO 90 J=l,2 
90 C(J)«A(J,3) 

C(l)»C(l)/A(l r l) 

SUM«A(2,1)*C(1) 

C(2)«(C(2)-SUM)/A(2,2) 

SUM*A(1.,2)*C(2) 

C(1)=C(1)-SUM 

RETURN 

99 PRINT*, 'TROUBLE IN SUBROUTINE FIT' 
RETURN 
END 

C 

SUBROUTINE ZEROFILL (C,IM,IE f N) 
DOUBLE PRECISION C (2048) ,D (2048) ,T 
M=IE-IM 
MP1=M+1 
DO 10 1=1, N 
D(I)=0.D0 



10 CONTINUE 

DO 20 1=1, MP 1 
J-IM+I-1 
D(I)=C(J) 
20 CONTINUE 
DO 30 I«1,N 
C(I)«=D(I) 
30 CONTINUE 
RETURN 
END 

SUBROUTINE TAPER (C,M,N) 
DOUBLE PRECISION C(2048) ,T 
T»l. D0/1000. DO 
DO 90 1=1,1! 

ir (C(I) .NE. 0.D0) GO TO 80 

J=I-1 

K«I+1 

C(I)*(C(J) +C(K))/2.D0 

IP (C(I) .LT. T) GO TO 95 

80 CONTINUE 

90 CONTINUE 
GO TO 97 

95 IK»I 

DO 96 J=IK,M 
K=J-1 

C(J)=0.95D0*C(K) 

96 CONTINUE 

97 CONTINUE 
ND2=N/2.D0 

DO 100 1=2 , ND2 

J=ND2+I 

K=ND2+2-I 

C(J)=C(K) 
100 CONTINUE 
RETURN 
END 

SUBROUTINE SMOOTH (C,M,N) 
DOUBLE PRECISION C(2048) 
DO 20 I«2,M 
J=I-1 

IF (C(J) .GT. C(I) ) GO TO 10 
K=I+1 

C(I) = (C{J)+C(K))/2.D0 
5 IF (C(J) .GT. C(I)) GO TO 10 
C(I)=0.999*C(I) 
GO TO 5 
10 CONTINUE 
20 CONTINUE 
ND2=N/2 
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DO 30 1=2, ND 2 

J-ND2+I 

K=ND2+2-I 

C(J)-C(K) 
30 CONTINUE 
RETURN 
END 

C 

SOBROOTINE REORDER (M,C) 
C IN ORDER TO USE THE POLLOWING PFT ONE MOST 
C FIRST* REORDER THE DATA. SEE JPS DISSERTATION 
C FOR DETAILS. 

COMPLEX*16 C(4096),T 

N=2**M 

MV2«N/2 

NM1=N-1 

3=1 

DO 7 1=1, NM1 

IF (I .GE. J) GO TO 5 

T=C(J) 

C(J)«C(I) 

C(I)=T 

5 K = MV2 

6 IF (K .GE. J) GO TO 7 
J=J-K 

K=K/2 
GO TO 6 

7 J=J+K 
RETURN 
END 

C 

SOBROOTINE FFT (K, A) 
C FOR DETAILS SEE JPS DISSERTATION 
DOUBLE PRECISION PI 
COMPLEX*16 A(4096) ,U,W,T 
PI=3.141592653589793D0 
N=2**M 
DO 30 L=1.M 

LE=2**L 

LEI-=LE/2 

0= (l.D0,O.D0) 

W=DCMPLX(DCOS (PI/LEI) ,-DSIN (PI/LE I) ) 
DO 20 J=1,LEI 

DO 10 I=J,N r LE 
IP=I+LEI 
T=A(IP) *0 
A(IP)=A(I)-T 
A(I)=A(I)+T 
10 CONTINUE 
U=U*W 



20 CONTINUE 
30 CONTINUE 

RETURN 

END 

C 

SUBROUTINE CONV (M,F,R,H) 
DOUBLE PRECISION X 

C0MPLEX*16 H(4096),P(4096),R(4096) 

Wl=FLOAT(N) 
DO 10 1=1, N 

H(I)»P(I)*R(I) 
10 CONTINUE 
DO 20 1-1, * 

H(I)-H{I)/W1 
20 CONTINUE 
RETURN 
END 

C 

SUBROUTINE DECON (M,F,R,H) 
DOUBLE PRECISION X 

SSES?* 16 P(4096 > »R(4096) ,H{4096) 

C IN FOURIER SPACE DECONVOLUTION IS SIMPLY A 
C COMPLEX DIVISION. J>iaPi,X A 

DO 10 1*1, N 

F(I)=E(I)/R(I) 
10 CONTINUE 

X=FLOAT(N) 
DO 30 1=1, N 

P(D=P(I)/X 
30 CONTINUE 
RETURN 
END 

C 

SUBROUTINE ITER (M,R,F,H) 
DIMENSION A(2048) ,B(2048) 

COMPLEX*16 R(4096),F(4096),H(4096},H1(40961 T 
CHARACTER* 7 NAME »oAl*U»cj ,T 

PRINT*, 'ENTER THE VALUE FOR THE PARAMETER WHICH* 
PRINT* , ' CONTROLS THE ALLOWABLE ERRORBETWEEN THE ' 
PRINT*, 'OBSERVED SPECTRUM AND F*R * 

» R S!'!E S=THE SQRT 0P THE SD « OF THE SQRS OF' 
PRINT* , 'THE DIFFERENCES BETWEEN THE H (I) S AND' 
PRINT*, 'THE F(I)*R(I)S. ' 1 ' 

PRINT*, '(A REAL NUMBER PLEASE.) • 
READ (*, ' (F13.7) ') EPS 

PRINT*, 'ENTER THE MAXIMUM NUMBER OF ITERATIONS' 
PRINT*, ' YOU WISH TO PERMIT. • 
PRINT*, • (AN INTEGER PLEASE) • 
READ (*, • (13) ') I END 
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N«2**M 

ND2-N/2 

ND4-N/4 

DO 30 I-1,IEND 

CALL REORDER (M,P) 

CALL PFT (M,F) 

CALL CONV (M,F,R,B1) 

CALL REORDER (M,P) 

CALL FPT (M,F) 

CALL REORDER (M,B1) 

CALL PFT (M,B1) 

RSQ'0.0 

DO 2 IA=1,N 

A(IA)=REAL(H(IA)) 

B(IA)-REAL(B1(IA)) 
2 CONTINDE 
V-0.0 
W=0.0 

DO 4 IAA=1,N 

V«V+A(IAA) 

W-W+B(IAA) 
4 CONTINUE 

DO 6 IAB*=1,N 

A(IAB)=A(IAB)/V 

B(IAB)=B(IAB)/W 
6 CONTINDE 

DO 10 IAC=1,ND4 

RSQ=RSQ+ (A ( IAC) -B ( IAC) ) * *2 

10 CONTINUE 
RSQ»RSQ**0.5 
WRITE (*,11) RSQ 

11 FORMAT (5X,«RSQ« *,E13.6) 

IF (RSQ .LT. EPS) GO TO 33 

F(1)=F(1) *(A(1)+2*A(2)+2*A(3) ) / (B(l) +2*B (2) +2*B(3) ) 
F(2)=F(2)*(A(1)+2*A(2)+A(3)+A(4))/(B(1)+2*B(2)+B(3)+B(4)) 
DO 20 IB«3,ND4 
J=IB-2 
K=IB+2 
TB-0.0 
TB1-0.0 
DO 15 II=J,K 
TB«TB+A(II) 
TH1=TB1+B(II) 
IS CONTINUE 

IF (TBI .EQ. 0.0) GO TO 20 
P(IB)«F(IB) *TB/TB1 
20 CONTINUE 
DO 25 IC«2,ND2 
J=ND2+IC 
K=ND2+2-IC 
F(J)=F(K) 
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25 CONTINUE 
30 CONTINUE 
33 CONTINUE 
I-I-l 

"JITE (*,35) I 

SSKS (5X ' ,THEN °- OP ITERATIONS IS 
END 



SUBROUTINE P ARAMS (N,A) 

SK'S A(204S) ' X(12) -■«•» -v.w 

K=1024 

00 20 1=1,12 
X(D-0.D0 
DO 10 J=l,5 
K=K+1 

X (D=X(I)+A(K) 
10 CONTINUE 
20 CONTINUE 
V=0.D0 

DO 30 1=1025,1084 
V=V+A(I) 
30 CONTINUE 

DO 40 1=1,12 

40 CONTINUE 

°0 60 1=2,10,2 
M(I)=0.D0 
DO 50 J=1026,1084 
V=PLOAT(J-1025) 
W=FLOAT(I) 
V=V**W 

50 CONTINUE 

tn /vl M (D=M(I)/59.D0 
60 CONTINUE 

M(2) =DSQRT(M(2) ) 
V-M{2) ' 

DO 65 1=4,10,2 
65 CONTIN& I) - M(1,/(V " 1 > 

READ («,'(A7)') name 

60 sr.sstf'iiSi, 1 
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100 FORMAT (9X,'M(10) *\D17.10) 

WRITE (90,110) 
110 FORMAT (5X,'5 CHANNEL GROUPS (CENTER FIRST) * 

WRITE (90,120) (I,X(I),I»1,9) 
120 FORMAT (10X, «X( ',11, •) = \D17.10) 

WRITE (90,130) (I, X(I) ,1=10,12) 
130 FORMAT (9X, 'X( ' ,12, ') = ',D17.10) 

CLOSE (90) 

RETURN 

END 

SUBROUTINE BACK (A, I MAX, N) 
DOUBLE PRECISION A(1024) ,C (2048) ,V 
COMPLEX*16 B(4 096) 
CHARACTER* 7 NAME 
ND2=N/2 
ND4=N/4 
DO 10 I»1,N 
C(I)=0.D0 
10 CONTINUE 
N1=ND2-IMAX 

IF (Nl .GT. ND4) N1=ND4 
DO 20 1=1, Nl 

J=IMAX+I-1 

C(I)=A(J) 
20 CONTINUE 
N2=IMAX-1 
N3=N-N2 
DO 30 J=1,N2 

I=N3+J 

C(I)=A(J) 
30 CONTINUE 
DO 40 1=1, N 

B(I)=DCMPLX(C(I),0.D0) 

40 CONTINUE 
M=ll 

CALL REORDER (M,B) 
CALL FFT (M,B) 

PRINT*,'WOULD YOU LIKE TO SEE THE FIRST FEW* 

PRINT* , * FOURIER COEFFICIENTS? ' 

READ (*, ' (Al) ') NAME 

IF (NAME .EQ. 'N' ) GO TO 43 

PRINT* , * HOW MANY? ' 

READ (*,*) IFC 

WRITE (*,41) (I,B(I) ,1=1, IFC) 

41 FORMAT (I4,2D22.10) 
43 CONTINUE 

PRINT*, 'WOULD YOU LIKE TO SEE THE FIRST FEW' 
PRINT*, 'VALUES OF THE POWER SPECTRUM? ' 
READ (*,' (Al) ') NAME 
IF (NAME . EQ. 'N' ) GO TO 46 



PRINT*, ' HOW MANY? • 
READ (*,*) IPS 
DO 44 I«1,ZPS 

A(Z)=B(I)*CONJG(B(I)) 

44 CONTINOB 

WRZTE (*,45) <I,A(Z),Z=1,ZPS) 

45 FORMAT (Z4,D22.10) 

46 CONTINUE 

PRINT*, 'DO YOU WISB TO DELETE ALL OF THE SINE' 

PRZNT*,'COEFFZCZENTS? ' 

READ (*,'(A1) •) NAME 

ZF (NAME -EQ. 'N* ) GO TO 51 

DO 50 1=1, ND 2 

B(Z)«(B(Z)+CONJG(B(Z)))/2.D0 

50 CONT1NDE 

51 PRINT*, 'DO YOO WZSH TO DELETE SOME OF THE COSINE 
PRINT*, •COEFFICZENTS? ' 

READ (*, ' (Al) ') NAME 

IF (NAME .EQ. 'N') GO TO 53 

PRINT*,' HOW MANY? ' 

READ (*,*) IDC 

DO 52 1=1, IDC 

c, ~ B(I)=(B(Z)-CONJG(B(Z)))/2.D0 

52 CO NT IN DE 

DO 152 1=2, IDC 
J=2050-I 

B(J)=(B(J)-CONJG(B(J) ))/2.D0 
152 CONTINUE 

53 PRINT*, 'INSTEAD OF DELETZNG SOME OF THE COSZNE' 
PRZNT*,'COEFS. YOU MAY PREFER TO JUST TRIM THEM.' 
PRZNT*,'THZS IS DONE BY SETTZNG THE FIRST N TERMS 

2fE!' EQDAL T ° THE N+1 TERM - *» MAY I NVOKE* 
PRINT*, 'THIS OPTION EVEN IF YOU HAVE PREVZOUSLY' 

«HSI!'! DELETED THE fZRST N COSZNE COEFFICIENTS. * 
PRZNT*,'DO YOU WZSH TO DO SO? ' 
READ (*, ' (Al) ') NAME 
ZF (NAME .EQ. 'N') GO TO 200 
PRINT*, 'HOW MANY? • 
READ (*,*) IT 
DO 190 1=1, IT 
J=IT+1 

B(Z)=DCMPLX(REAL(B(J) ) ,O.D0) 
190 CONTINUE 
200 CONTINUE 

PRINT*, 'DO YOU WZSH TO DELETE SOME OF THE SZNE' 

PRINT*, 'COEFFICIENTS? ' 

READ (*,' (Al) ') NAME 

IF (NAME . EQ. 'N' ) GO TO 55 

PRINT*, ' HOW MANY? ' 

READ (*,*) IDS 

DO 54 1=1, IDS 



B(I)=(B(I)+CONJG(B(I)))/2.D0 

54 CONTINUE 

DO 154 1=2, IDS 
J-2050-I 

B(J)-(B(J)+CONJG(B(J)) )/2.D0 
154 CONTINUE 

55 CONTINUE 

PRINT* , ' DO YOO WISB TO SAVE THE FOURIER COEPS. 

PRINT*,' B(X) PRIOR TO THE TAPER IN THE WINGS? 

READ (*,'(A1) *) NAME 

IP (NAME . EQ. *N* ) GO TO 57 

PRINT*, 'ENTER THE FILENAME. ' 

PRINT*,' • 

READ (*,»(A7) •) NAME 

OPEN (90, FILE" NAME, PADs 'YES') 

WRITE (90,56) ( I ,B ( I) , 1=1 ,2048) 

56 FORMAT (2 ( 14 ,2D22 .10) ) 
CLOSE (90) 

57 CALL REORDER (M,B) 
CALL FFT (M,B) 

DO 60 1=1, ND2 

C(I)=REAL(B(I)) 
60 CONTINUE 
V=O.D0 

DO 70 1=359,484 
V=V+C(I) 
70 CONTIBUE 
V-V/126. 

PRINT*, 'DO YOU WANT A D.C. CORRECTION USING A' 
PRINT*, 'LEAST SQUARE FITTING IN THE TAILS? ' 
READ (*,' (Al) ') NAME 
IF (NAME . EQ. 'N') V=0.D0 
DO 80 1=1, ND 2 
C(I)=C(I)-V 

IF (C(I) .LT. 0.D0) C(I)=0.D0 
80 CONTINUE 

DO 90 1=1, ND4 
J=(3*ND4)+I 
A(I)=C(J) 
90 CONTINUE 
N4=ND4+1 
DO 100 I=N4,ND2 
J=I-ND4 
A(I)=C(J) 
100 CONTINUE 
IMAX=513 
RETURN 
END 

SUBROUTINE NOBACK (N,B,IMAX) 
DOUBLE PRECISION B (1024) ,A (2048) ,V 



C0MPLEX*16 C(4096) 
CHARACTER* 7 NAME 
M=ll 
ND2=N/2 
ND4=N/4 
DO 10 1=1, N 
A(I)=0.D0 
10 CONTINDE 
N1=ND2-IMAX 

IF (Nl .GT. ND4) Nl=ND4 
DO 20 1=1, Nl 
J=IMAX+I-1 
A(I)=B(J) 
20 CONTINDE 
N2-IMAX-1 
N3=N-N2 
DO 30 J=1,N2 
I=N3+J 
A(I)=B(J) 
30 CONTINDE 
V=0.DO 
DO 33 1=1, N 
V=V+A(I) 
33 CONTINDE 
DO 36 1=1, N 
A(I)=A(I)/V 
36 CONTINDE 

OPEN (90,FILE='RAWDATA' ,PAD=* YES* J 
WRITE (90,70) (I,A(I) ,1=1, N) 
CLOSE (90) 
DO 40 1=1, N 

C(I)=DCMPLX(A(I) ,0.D0) 
40 CONTINDE 

CALL REORDER (M,C) 
CALL FPT (M,C) 
DO 50 1=1, N 

A(I)=(Xl)*CONJG(C(I) ) 

50 CONTINDE 

PRINT* , ' ENTER THE FILENAME FOR E(X). 

PRINT*,' • 

READ (*,»(A7) ') NAME 

OPEN OO^ILE-NAME^ADa'yES') 

WRITE (90,60) ( I ,C(I) ,1=1, 2048) 

60 FORMAT (2 ( 14 ,2D22 .10) ) 
CLOSE (90) 

PRINT* , ' ENTER THE FLENAME FOR THE P.S. 
PRINT*, ' • 
READ (*,'(A7) ') NAME 
OPEN (90 , FILE=NAME,PAD= 1 YES' ) 
WRITE (90,70) d,A(I) ,1=1,1024) 
70 FORMAT (5(14, D22 .10) ) 



CLOSE (90) 

RETURN 

END 

cccccccccccccccccccccccccccccccccccccccccccccccccccc 



APPENDIX G 



In this Appendix the theoretical relationship between the vacancy 
concentration, C v , and the concentration of dislocation cores, C d , is 
demonstrated, (from Cotterill et ak 1972). The starting equation is 

C v = A exp(-E vf /KT) (6<1) 
where C v Is the vacancy concentration at thermal equilibrium 

E v f is the vacancy formation energy 
and KT has its usual meaning 

Now the concentration of vacancies after the quench. c\ is related to C v by 



v 

C.' ■ nC. 



L v = n S (6.2) 
where n is the quench efficiency which is assumed to be nearly unity 

After the low temperature anneal the concentration of dislocation 
loops, Cg, is given by 

C * s vC v ' (G.3) 

where v is the inverse of the number of vacancies required to form an 
average loop 

The dislocation concentration, C d , in terms of an atomic fraction is then 

C d - C^O/d,, (G>4) 
where D is the average loop diameter 
and d 0 is the nearest neighbor distance. 
It can be shown that 
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v - Zfl d 2 /n0 2 (G.b) 
Combining Equations (G.3), (6.4). and (G.5) yields 

C d = 20 d Q C;/D (G.6) 



APPENDIX H 



In this Appendix the derivation of the defect specific lineshape 
parameter "R« is recreated from the work of Mantl and Trifthauser (1978). The 
following quantities require definition: 

S.W - are the values of the "S" and "W parameters for the sample 
being analyzed. 

S t .W t - are the values of the "S" and V parameters for a sample in 
which all of the positrons are annihilating from the trapped 
state. 

S B .W B - are the values of the "S» and »W" parameters for a sample in 
which all of the positrons are annihilated from the bulk 
(defect free). 

P t - the probability of a positron annihilation from the trapped 
state. 

P B - the probability of a positron annihilating from the bulk, 
u - the trapping constant. 
C - the trap concentration. 
*B - the annihilation rate from the bulk state. 
Now, it is known that 

S = S B P B + S t P t (H#1) 

and 

P B + h = 1 (H.2) 
(143) 



Also F t ' «*C/(A B ♦ pC) 

Substituting Eq. (H.3) Into Eq. (H.l) yields: 

PC 



Similarly 



Rearranging Eq. (H.4) and Eq. (H.5) yields: 

S - S B = (S t - S B )wC/(A f + U C) 
W - W B = (W t - W B )„C/(X f + „C) 

Taking the ratio of Eq. (H.6) to Eq. (H.7) yields: 



R = 



I s 




w 





S t " S B 



W - U 

t B 



Clearly the ratio on the far right has a constant value for 
a given material. 



144 
(H.3) 

(H.4) 



(H.5) 

(H.6) 
(H.7) 

(H.8) 
a given defect in 
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A positron annihilation lifetime spectroscopy study of porous silicon 
using a continuous lifetime fitting algorithm 

P. M. DERLET, T. C. CHOY 

Department of Physics, Monash University, Clayton, Victoria 3168, Australia 



Porous silicon is a highly complex system containing 
disorder that exists over a range of scales. Depend- 
ing on the anodization current and time, and the 
dopant level and type, the resulting porous layer can 
have a porosity of up to 80% with a unidirectional 
columnar-like structure for p" silicon wafers [I]. 
The remaining silicon varies in size with character- 
istic lengths ranging between from 3-4 to 15- 
30 nm. Such small crystalline structures have an 
enlarged bandgap due to the quantum confinement 
effect The role of confinement in the observed 
increase in radiative inter-band transition efficiency 
is however still unclear. How the broken silicon 
bonds are terminated can directly affect the nature of 
the inter-band transitions. For example, hydrogen 
passivation "sweeps" the dangling bond states out of 
the band gap region; reducing the available non- 
radiative recombination paths and thereby increasing 
the efficiency of the radiative transitions. This, 
however, is not the complete picture since the 
oxygen content of the porous region has been found 
in some cases to amplify and change the photo- 
luminescence properties [2]. The electro-chemical 
processes occurring on the surface during anodiza- 
tion are not yet clearly understood and further 
experimental characterization is needed 

Positron annihilation lifetime spectroscopy (PALS) 
provides a unique non-destructive probe into the 
disorder inherent in porous silicon on the micro- 
scopic scale, through the identification of impurity 
concentrations and types; and on the mesoscopic 
scale at which the bulk of the porosity exists through 
the study of positronium lifetimes. Itoh et al. [3.4] 
and Dannefaer ex at. [5] have done such an 
investigation using a discrete lifetime analysis. Itoh 
zi al. obtained a satisfactory three -component fit, 
finding a bulk silicon component (at 230 ps), a 
defect component (at 600 ps), which they interpreted 
as arising from either silicon mono- vacancies and/or 
oxygen divacancies within an Si- O2 layer [6 J, and a 
longer 25 ns lifetime due to ortho-positronium 
annihilating in the porous region. Dannefaer et al. % 
on the other hand, determined lhat a five-component 
fit was better, with additional lifetimes arising from 
para-positronium (at 1 25 ps) and another ortho- 
positronium component at 2.5 ns. They attributed 
the latter to ortho-postronium annihilating in larger 
vacancy defects either within the remaining silicon 
or the oxide layer. In this letter we report on a PALS 
investigation of pGrotis silicon usini* a continuous 
lifetime fining algorithm. 



Our motivation lies in the underlying disadvantage 
in discrete lifetime fitting algorithms where the 
number of components must initially be assumed 
since in general a realistic spectrum does not 
uniquely determine this number. This becomes 
particularly apparent when looking at highly dis- 
ordered systems where, indeed, the notion of a 
discrete spectrum may be invalid. These problems 
can in part be overcome by using a continuous 
lifetime fitting algorithm; for example, performing 
an inverse Laplace transform on the lifetime 
spectrum, thereby obtaining a probability density 
function of characteristic lifetimes. The Fortran 
program CONTIN [7-9] is such an algorithm. 
Admittedly, use of the inverse Laplace transform 
does not take into account the instrument resolution 
function of the apparatus. For the determination of 
sub-nanosecond lifetimes this entails a distinct loss 
of accuracy. For well-separated lifetimes that are 
significantly larger than the typical full width at half 
maximum (FWHM) of the equipment, the effect is 
reduced and overshadowed by the loss of accuracy 
due to the worsening statistics typical for these 
components. To date, our work appears to be the 
only systematic study of porous silicon using a 
continuous lifetime fitting algorithm. This is the first 
step towards achieving the goal of determining 
intrinsic lifetime distribution widths for positrons 
annihilating within the porous silicon layer, and thus 
its electronic structure. 

For the measurement of the annihilation lifetime 
spectra of porous silicon we employed a standard 
PALS rig at 200 ps per channel. Such a large timing 
window is needed to ensure that all lifetime 
components decay completely into the background. 
The positron source was 22 NaC! evaporated onto a 
6/*m thick titanium film. Using our PALS rig at 
50 ps per channel and the discrete lifetime fitting 
program PFPOSFIT [10], we obtained for bulk 
silicon, 221 ± I ps at 97 ± 1% concentration, with 
a source component lifetime of approximately 
429±25ps at 3 ± i% concentration (the iterared 
instrument resolution had an FWHM of 240 ps with 
long and short time tails of 50 and 20 ps. 
respectively. 

The porous silicon samples for this investigation 
were made from 200 urn p~ doped wafers, anodized 
for 30 min at 20 mA cm" 2 . The electrolyte used was 
a 35% ethanoic solution of HF. Under illumination 
by an ultraviolet (UV) lamp this sample glowed a 
Jeep red. The initial thickness of the porous layer 
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was approximately 50-60 urn. A high resolution 
scanning electron microscopy (SEM) analysis per- 
formed approximately 6 months after the PALS 
experiment revealed a layer that was approximately 
10/im thick, indicating that structural damage had 
occurred (Fig. 1). The sample nevertheless still 
glowed red - at a reduced intensity. Fig. 1 reveals the 
large-scale columnar structure typical of p~ doped 
porous silicon samples; the remaining (porous) 
silicon constituting the much thinner connecting 
membranes between the micrometre-sized pores. 

All PALS runs were taken at room temperature 
(24 °C) in an environment where the temperature 
could be held constant to within a degree. For each 
sample, a varying number of spectra were taken (see 
Table I and Fig. 2), each of which had a maximum 
peak height of 10 6 counts. Statistical analysis across 
the spectra for a particular sample indicated that to 
within the noise level, there was no detectable 
difference between individual runs, justifying the use 
of the mean spectrum in the analysis using CONTIN. 
The associated information content analyses [II, 12] 
are also given in Table I. Here, the Shannon number 
indicates that in principle up to eight distinct lifetime 
peaks are resolvable, and the resolution factor 
suggests that the minimum ratio between these 
resolvable lifetimes ranges between 2 and 3. 

Since CONTIN assumes a spectrum constructed 
from exponentials, we must choose an appropriate 
zero-time channel from which the analysis can 
begin. This must be just beyond the peak of the 
spectrum, yet not too close to include the "rounding 
off* effect of the finite instrument resolution. The 
resulting inverse Laplace transform found by CON- 
TIN is not very sensitive to the choice of the zero- 
time channel as far as determining distinct compo- 
nents. It is. however, sensitive to determining their 



intensity, since a shift dt in the zero-time channel 
will atTect the intensity of a component r by the 
approximate factor exp(-d//r). 

Such an analysis for the 200 ps per channel bulk 
silicon wafer spectrum is displayed in Fig. 2 for two 
different zero-time channels; approximately 400 ps 
(solution I) and 800 ps (solution 2) beyond (he 
instrumental zero-time channel. Solution I contains 
a dominant broad component centred at approxi- 
mately 220 ps, which can be attributed to- the bulk 
silicon component. Two long time components can 
be seen, occurring at approximately 1.5 and 45 ns, 
the shorter of which can be associated with the 
common ortho-positronium "pick oflT signal present 
in semiconductors [13J. The extra long 45 ns 
component does not seem to have been encountered 
in the literature. We do not know at this stage the 
origin of this lifetime peak; however, as we shall see, 
it also appears in the porous silicon sample at an 
increased intensity. In solution 2, the bulk silicon 
component becomes broader and less intense, the 
ratio of its peak height with that of solution 1 being 
approximately exp(r3i/220) where <Jf»400ps. The 
ortho-positronium peaks at 1.5 and 50 ns acquire 
more dominance, with their peak heights largely 
unchanged from solution 1. In addition, we observed 
a shorter, approximately 125 ps, component, which is 
most likely to be the accompanying para-positronium 
signal. That this was not resolved in solution I 
suggests that the "rounding off" effect was not 
completely removed at 400 ps beyond the approx- 
imate zero-time channel. 

Fig. 3 displays the CONTIN solution for the 
porous silicon spectra. Again, solution 1 and solution 
2 are for time shifts of approximately 400 and 
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Figur* I High resolution SEM 
jsed in the PALS investigation. 



picture of the porous silicon sample 



Xgur* 2 PALS lifetime spectra of (a) bulk silicon and <bj porous 
siheon taken at 200 ps per channel. 
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Figure J Characteristic lifetime probability density determined using 
CONTIN for the bulk silicon sample taken at 200 ps per channcL (a) 
Solution J: zero-time approximately 400 ps beyond the instrumental 
zero-time channel, (b) Solution 2: zero-time approximately 800 ps 
beyond the instrumental zero-time channel. 



800 ps beyond the instrumental zero-time channel, 
respectively. Inspection of solution 1 indicates five 
distinct components. Component 1 is centred at 
approximately 200 ps and can be attributed to the 
positrons annihilating either within the silicon 
substrate or in the remaining silicon in the porous 
material. Component 2 is centred at approximately 
300 ps and has also been encountered by both Itoh et 
al. and Dannefaer et al. Their measured lifetimes for 
component 2» however, ranged between 400 and 
600 ps. CONTIN analyses of simulated spectra 
obeying Poisson statistics indicate that this discre- 
pancy (and the reduced bulk lifetime) may be due to 
the 200 ps per channel resolution and the neglect of 
the instrument resolution function. Components 3, 4 
and 5, occurring respectively at approximately 2, 10 
and 50 ns. are the long lifetime components due to 
ortho-postxonium. In solution 2, the sub-nanosecond 
components are reduced in intensity, and crystalline 
silicon and defect components (at 300 ps) are no 
longer separately resolvable, merging into a single 
broad peak centred at approximately 250 ps. The 
long time components now become broader and of a 
reduced height. This indicates that the widths we 
have observed are not intrinsic to the material but 
rather an arte tact of the fining process and the signal 
to noise ratio. Marshall [14] has concluded that to 
determine accurately the intrinsic widths of peaks, 
this ratio must be at least 400. Inspection of Fig. 4 
indicates a signal to noise ratio of less than 40 for 
the positronium decay modes. 

The present method of analysis is unsuitable to 
obtain accurate absolute intensity information since 
the zero-time channel cannot be estimated beyond 
the 200 ps per channel limit Furthermore, modifica- 
tions to CONTIN suggested by Gregory [12, 15-17] 
that take into account the finite resolution (and 
iherefore the zero-time channel), by obtaining the 
probability density profile with respect to a reference 
spectrum (in our case, bulk silicon), have proved 
unsuccessful. This is primarily due to the poor 
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Ftgurr 4 Characteristic lifetime probability density determined using 
CONTIN for the porous silicon sample taken at 200 ps per channel, (a) 
Solution I: zero-urne approximately 400 ps beyond the instrumental 
zero-time channel, (b) Solution 2: zero-time approximately 800 ps 
beyond the instrumental zero-time channel. 

definition of the bulk silicon reference spectrum at 
200 ps per channel. 

Using the PALS data collected from sample A, we 
have confirmed and extended the findings of both 
Itoh et al. and Dannefaer et al. In particular, we 
resolved three rather than two ortho-positronium 
components, suggesting that there may be an 
additional intermediator^ scale of porosity in which 
ortho-postronium annihilates. This has a lifetime 
indicative of hydrogenated amorphous silicon [18], 
suggesting ortho-positronium annihilation again 
within the amorphous layer. The integrated areas 
(corrected for zero-time shift) of the positronium 
peaks in solutions I and 2 show that components 3 
and 5 occur in approximately the same proportion; 
Dannefaer et al. % on the other hand, found that the 
—2.5 ns component occurred an order of magnitude 
less than the longest time ortho-positronium compo- 
nent The present analysis also establishes the 
existence of a very weak ortho-positronium compo- 
nent in the pre-anodized wafers at a time scale 
approximately equal to the longest time ortho- 
positronium component seen in porous silicon, using 
theory [19] that has traditionally been employed in 
the calculation of free volumes in polymer materials. 
In porous silicon Itoh et al. associated a void length 
scale of approximately 1 nra to this ortho-positro- 
nium component, suggesting that irregularities of 
ihis rnagnifcide exist on the pre-anodized wafer and 
that these may. in part, be the catalyst for the initial 
pore formation process. 

The results presented here further demonstrate the 
rich variety of disorder inherent in porous silicon 
and indicate the potential of a continuous lifetime 
analysis of PALS spectra to provide a high resolution 
probe into such structures. If another order of 
magnitude increase in signal to noise ratio can* be 
achieved then there exists the possibility of resolving 
intrinsic peak widths for porous silicon samples, 
thereby extending the PALS method to resolving 
pore size distributions using the ortho-positronium 
components, and conceivably giving insight into the 
electronic structure of the confined electrons wiihin 
the remaining silicon crystallites of the porcus layer. 

1951 
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Abstract 

The study of decotivolubon by Van Gttert and Cold iterative algorithms and their use in the processing of experimental 
spectra of Doppler fr mfn t f""g of the annihilation fine in positron annihilation measurement it described. By comparing 
results from both algorithms it was observed that* the Gold algorithm was able to efirmnate linear instability of the 
measuring equipment if one oses the 1274 keV °Na peak, that was measured simultaneously with the annihilation peak, 
for (^convolution of annihilation peak 511 keV. This permitted the measurement of small changes of the annihilation peak 
(eg. ^parameter) with high confiden ce . The dependence of x-ray-like peak parameters on the number of iterations and 
the ability of these algorimms to distinguish a y-ray doublet with different intensities and positions were also studied. 




ELSEVIER 



1* Introduction 

Positron annihilation (PA) as one of the non-destructive nuclear methods is used for the study of structure 
of different materials and it is especially efficient in determining free volumes or defects in bulk material. The 
main characteristics, that are determined by PA, are lifetimes of positrons in matter, Doppler broadening of 
annihilation line (DBAL), and angular correlation of annihilation y-rays. The last two characteristics give us 
information about the momentum distribution of electrons in matter. The more common angular correlation 
method is approximately one order more accurate, but it is more time consuming and requires a more intensive 
positron source. On the other hand DBAL is simple and not time consuming. It is mainly used for determining 
the relative changes in the annihilation peak (e.g* S-parameter, defined as the ratio between area of central part 
of the peak to the full area of the peak) depending on temperature, pressure etc. After deconvolution of the 
annihilation peak using a resolution function of measuring system, one can obtain the momentum distribution 
of electrons in the sample; If we are only interested in relative changes of the S-parameter, deconvolution is 
not necessary. But as we have seen in our work, deconvolution can partly eliminate the spurious effects due to 
instability of the measuring system. Another possibility to eliminate an instability of the system is the use of 
an expensive spectrum stabilizer unit 

Generally, deconvolution algorithms have recently found many applications in various domains of experi- 
mental science. They have been applied to such different problems as improving resolution in spectroscopy and 
measuring the thickness of multilayer structures. There is also the long-standing problem of determination of 
position and intensities in /-ray multiplets ft]. 

2. Theory 

The relationship between a measured value x(t) and the raw result of measurement y(/) can be described 
by a convolution-type integral equation 

OI6S-9O02/97/SJ7O0 Copyright (g> 1991 Htevier Science B.V. AU rigj»» ittmcd 
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x(T)h(t-r) dr. 



(1) 



where A(l) is an impulse response. Knowledge of the instrumental function A(r) is usually required. For a 
discrete system, Bq. ( I ) can be written as 



y(i) = £ h(i - k)x(k) i = 0. 1.. ; . ,2N - 2 , 



(2) 



where N is the number of samples of vectors A, x. 

lite impulse response has a finite length. Therefore we will assume that A(i) = 0 for / < 0 and / ^ N. Then 
Bq. (2) can be written in a matrix form: 



y(0) 
.v(l) 
>(2) 



.y(2AT-2) 



A<0> 

Ml) 
A<2) 



0 
A(0) 
A(l) 



0 
0 
A(0) 



A(W-l) h(N-2) HN-3) 
0 A(tf-l) A(JV-2) 
0 0 h(N-l) 



x<0) 
x(l) 
x<2) 

x(AT-I) 



(3) 



y = Hx. 



(4) 



It means that the columns of H i 
both sides of Bq. (4) by H T gives 

H T y = H T Hx. 



! represented by vectors A mutually shifted by one position. Multiplying 



or 



Ji^Hix, 



(5) 



(6) 



where Hi is Tbeplitt matrix [23]. 

Solution of the linear equation system (6) (vector x), under the condition that die output vector y and the 
matrix of the impulse response H is known, is a problem of decon volution. The output vector of system y is 
affected by noise that accompanies each measurement. The existence of this noise strongly affects the process 
of deconvolutiom, and can tend to difficulties in solving the linear equation system (6). 

The wtove-fonauiaied problem of input recs as trucsoa is as a rate iH-csndidoncd, i.e. die estimates J»(f) of 
x(r) satisfying Bq. (1) are extremely sensitive to errors in the measured data jr(r). It is expressed by the fact, 
thai the matrix Hi is almost singular. The direct inversion of H, for solving x cannot lead to a stable solution. 

Therefore, in order to solve this problem, the method of regularization must be included. This means that 
origin^ problem is replaced by an approximate one which solutions are significantly less sensitive to errors in 
the data y(t). 

The Van Cittert iterative method of deconvoludon is widely applied in different areas, for example in 
spectroscopy or in image processing (2,4). The Van Cittert algorithm of deconvolution is described in detail in 
Ref. [4], so we will describe it only very briefly. 

Its basic form for a general linear discrete system is 



x**'»x*+MJ'-Ax*). 



(7) 
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where A is system matrix, k represents the number of iterations and /* is the relaxation factor. The convergence 
condition of Eq. (7) is that the diagonal elements of the matrix A satisfy 

An> £ A u . i = 0J /V- i. (8) 

It is obvious that such a diagonal element dominance is rare in physical problems. However, the deconvplution 
algorithm in Eq. (7) can be modified in a such way that it will satisfy the conditions of convergence. Hence, 
Eq. (7) becomes 

^•= M jP+(B-AiA)jr*=/tf + Dx*. (9) 
where B is a unit matrix and 

D = E-M (10) 
Under the condition that x9 = p, the successive substitutions give 

Assuming that Ao 9 At are eigenvalues of A, then (I - M>)»(* - Mi) (I are 

eigenvalues of D. If 

Km (l-M)*-0. ' = 0,1 tf-1. (12) 



then 

Kim 0* = [0J , 



and 

lim jr'-A"'*-*. (13) 
From equation (12) this implies that the necessary and sufficient conditions of convergence are 

|l-#Orl<l. '«M < l4 > 

If we define A* and iu conjugate A? as 

A*-»fli+jfti« -*r »<%-j*#. 
then the convergence condition (14) becomes 

MCM*? + *?)-2<i<]<0. 1=0,1 Af- 1- (15) 

Inequality ( 15) gives two bounds for /* 

A=0. M=-2~J_. i-0il.;...*-l. (16) 
af + of 

conditions determine the bounds of the /i coefficient Unfortunately these conditions are not fulfilled for 
all ca ses However* if the system matrix A is positive definite the convergent solution always exists. So we 
settle the algorithm in a such way that the eigenvalues A* will be positive, real numbers. 

Let us return to Eq. (5). Matrix H T H is symmetric, so its eigenvalues are real. The eigenvalues of matrix 
(H T H)(H T ff) are squares of eigenvalues of maSix H T H and therefore must be positive. Eq. (5) becomes 

(H T HH T ) / = (H T HH T H)jr t (17) 



» 
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and the iterative algorithm of dcconvolution becomes 
x ut = x* + m( (H T HH T ) j - (H T HH T H)x*J t 
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or 



(18) 



(19) 



By this we ensure the existence of a common interval of solution for inequality ( 14) and convergence of the 
deconvolution algorithm. Eigenvalues A< are real, positive numbers, so from Eq. (16) for /* we can write 



0<m<t— . 

Amu 

where U is the greatest eigenvalue of H' 

Aum = max(Ao»Ai A»_i). 

Now we determine the maximum eigenvalue A^, For eigenvalues of system ( 17) we can write 

HjcbAjX, t*0,l iV-1. 

Ifxj is the biggest absolute element in x, then from Eq. (22) we get 

N-l 

^H^-Aix, 



(20) 



(21) 



(22) 



(23) 



or 



Then 



*7 



« = 0.1 N-l. 



(24) 



(25) 



In practical cases we do not know the biggest element in x. We determine the value of /W, as the maximum 
value from Ah determined by Eq. (25) eg., from the sum of absolute values of rows in matrix H\ 

This is the base of the Van Gttert algorithm of deconvolution. Now we introduce, in analogy with Eq. (24), 
a local variable relaxation factor 



4 



(26) 



and we use it in Eq. ( 19). For the ith element of vector x** 1 we get 

5X4 L - 



(27) 
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or 

^-UZT—^- (28) 

Eq. (28) is the Gold algorithm of deconvolution [5J. It is an extension of Van Cittert 's iterative method. 
The advantage of Gold's method is that the solution x is positive. 



3. Theoretical simulation 

To test an applicability of Van Citteit's and Gold's methods of deconvolution for /-rays spectra we performed 
a theoretical simulation. We convoluted the experimental resolution function (RF) obtained from a positron 
annihilation experiments with different kinds of Gaussians describing y-ray peaks. The simulated spectra were 
used as an input to both Gold's and Van Gteert's iterative deconvolution procedures. Simulated spectra included 
a noise component due to an experimental resolution function. 

3.1. The dependence of output peak parameters from the deconvolution on the number of iterations 

Because the decon volutions are using the iterative method, k was interesting to study the dependence of the 
output from die deconvolution on the number of iterations (AO. N was varied in the range 200-50000. 2000 
iteration steps lasted 40 min on a FC486, 33 MHz for 200 channels of the spectrum (a version with elimination 
of the multiplication of zero matrix tenns took appro*. 15 min). 

We analyzed the peak which was' a convolution of a Gaussian with FWHM = 13 channels and a resolution 
function with a FWHM = 10A The resulting peak had a FWHM ■ 16.7. which was a value in the range of 
our experimental values for the FWHM of annihilation peaks measured in OBAL. The deconvolved peak was 
fitted and the result should be one Gaussian with a FWHM » 13. 

Fig. 1 shows the dependences of the x* value and FWHM on the number of iterations for the algorithm 
of Vhn Cittert and Gold. As we can see from this figure; the resulting deconvolution fit for the Van Cittert 
algorithm was much better for two Gaussians than for one until 14 000 iterations, then the x 2 became practically 
the same. 

For the Gold algorithm the x 1 of the deconvolution fit was very good for one Gaussian and reached a 
minimum value for 2000 iterations. The increase of x* with the number of iterations (for both Van Cittert and 
Gold methods) is a consequence of the round-off errors which arise during the computations with real numbers 
represented in standard single precision (REAL in Fortran or FLOAT in C) and a large number of these 
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operations It was verified by a calculation (for the Gold method) using the double precision floating-point 
representation of real number*. The results of the double-precision calculation are shown in Fig I. The fast 
convergence and the stability of the solution can be observed. It clearly shows the possible errors arising from 
using a finite computer word length of and a large number of iteration steps. 

For the Van Cittert algorithm, it can be seen that by increasing the number of iterations the first Gaussian 
became more and mote dominant and approached the input Gaussian in convolution. The same dependence can 
be observed in the intensity of Gaussians. 

For the Gold algorithm* the results were much better. From Fig. 1, it is clear that the FWHM from the fit 
of the deconvoluted peak reached the correct value of 13 for the number of iterations larger than 1000. The 
variance of the FWHM from this fit was < 0.03% and of the amplitude < 0.005% for ^ 1000. 

For the Van Gttot algorithm and N > 2000* the variance of the FWHM from the fit was < 0.3% and of the 
amplitude < 0.05%. 

In Fig. 2 the weighted differences of the original input Gaussian and fitted results of the deconvoluted peak 
after N » 2000 iteration steps axe shown. The difference is evident For die Gold deconvolution the variance 
was still less than 7.5% even for the ±25 displaced channels from the center of the peak. For the large number 
of iterations (H > 30000) there was already deformation in deconvoluted peak, which could be seen as sharp 
peaks in the dependence of weighted differences. This deformation is a consequence of round-off errors in the 
computer. 

Fig. 3 shows typical results of the deconvolution procedures. The peat deconvoluted by the Gold's algorithm, 
is almost identical with the original Gaussian after carrying out 2000 iterations. Van Quest's procedure showed 
negative values and there were oscillations on both wings of the peak on the level of 2% of amplitude of 
peak. Because Fig. 3 is in a logarithm scale the negative values are not displayed. Gold's deconvolution did not 
produce negative value* and oscillations were very small (< 0.01%). 

J.Z The deconvolution of Gaussian* with different FWHM 

To determine the influence of the FWHM of the convoluted Gaussian on the result of deconvolution, we 
also made a convolution of the Gaussian with FWHM » 20 with RF. Then the difference in parameters of the 
resulting deconvoluted peak after 2000 iterations was compared with the difference for FWHM » 13. 

The Van Cittert's deconvolution of convolution of the original Gaussian with a FWHM » 20 was in very 
good agreement with the original Gaussian, x 2 from the fit was 0.04 and the FWHM was exactly the same. 

The weighted difference of the fit of deconvoluted and original Gaussian was < 10% for counts > 0.1% 
of amplitude. The situation was worse for the convolution of the Gaussian with FWHM = 13. The FWHM 
of the deconvoluted peak was larger than about 1.8%. But the worst fact was that the x 1 from the lit was 
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very large (Fig. 1) and the deconvolution was better described by two Gaussians (x 2 = 1.8). That means that 
the deconvolution changed the shape of the peak in such a way that the deconvoluted peak could be better 
described by two Gaussians though there should be only one. The weighted difference of the deconvolution 
fit and the original Gaussian was ^ 10% for counts > 3% of the amplitude. The parameters from the fit of 
deconvoluted peak woe in excellent agreement with the original Gaussian. Problem was that for a Gaussian 
with FWHM = 13 there was a better x 2 for two Gaussians as for one. 

For the Gold deconvolution agreement with the original Gaussian was very good also for Gaussian with 
FWHM = 13 as we can see from Fig. 2. What was important was that there are no deformations in deconvoluted 
peak because the ^ for one Gaussian was very small (Fig. 2). 

3.3. The deconvolution of the sum of two Gaussians with different intensity and position 

Tb verify the possibility of distinguishing components of a doublet (sum of two peaks) by deconvolutioa, 
we made convolutions of the sum of two Gaussians with intensity ratios 1:1, 1:0.8, 1:0.5, 1:0.2 and mutual 
displacements of I, 5,. 10, IS. 20 channels. Both input Gaussians had FWHM » 13, the number of. iteration 
steps N was 2000 and the deconvoluted peaks were fitted. We computed the weighted difference of these fitted 
parameters with the parameters of the original input Gaussians. 

For the Van Cittert's deconvolution we got the following results: 

If the peak displacement was 15 and 20 channels there was no problem to fit the deconvolution results. The 
variation from the parameters of the input Gaussians was small. < 2%. For a displacement of about 10 channels, 
the fitting procedure for the ratio of intensities 1:02 failed. The variance was < 5% for ratios l:l P l:0.8 and 
< 19% for 1:0 J. For a displacement of about 5 channels we were not able to get a fit for 1:02 and 1:0.5 
ratios. The variances were < 42% (the largest for the amplitude, for the remaining parameters < 10%). So for 
peak displacements larger than the FWHM of the peaks, the resolution was good also for the ratio of intensity 
of 20%. For a displacement smaller as the FWHM of the peak, the resolution was possible only for intensities 
close to each other. For a small displacement (< 10%), the fit was impossible. 

For the Gold deconvolution we also got a fit for a displacement of 10 channels and ratio 1:0.2 and a 
displacement of 5 channels and ratios 1:0.5, 1:02 with variances < 27%. For a displacement of 1 channel the 
fitting procedure failed. So we can say that Gold's algorithm is much more able to distinguish a doublet and 
give us mora precise results than the Van Gttert's algorithm. 

In Fig. 4 convolutions of Gaussians and the results of the deconvolution procedure described above can be 
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4. Experimental spectra 

All these studies of the convolution of Gaussiaas and following deconvolution were aimed at finding out 
the possibility and limitations of the deconvolution methods and their possible use in the deconvolution of 
experimental y-ray spectra. Because the methods are quite reliable (especially Gold's method) we have used 
Gold's deconvolution method for our experimental spectra. * 

The effect of DBAL is quite smallscale and one needs a very stable measuring system for its measurement. 
We can get information about its stability from the 1274 keV peak, which is measured simultaneously with the 
511 keV annihilation peak if we use °Na as a radioactive source; 

In our experiment we used an HPGe detector with a resolution of 136 keV at 570 keV with a 207 Bi source. 
The energy calibration of the system was 163 eV/channeL The number of counts at the annihilation peak 
was «630000 and the number of counts at the 1274 keV peak was « 1 8000<X In Fig. 5 we can see the 
S-parameter with physical information about the momentum distribution of the electrons and the ^parameter 
of the 1274 keV peak, which carries information about the stability of the measuring system. If there was no 
instability during die measurement, it should be a constant As we can see, the ^parameter of the 1274 keV 
peak was not constant, so the measuring system was not stable. So we do not. know if the structure in the 
5-parameter of 51 1 keV peak was really a consequence of this instability or if it had a real, physical meaning. 
Of course, there is a possibility to try to use the 5-parameter of the 1274 keV peak as a correction on the 
instability. eg. by calculating the ratio R = 5(511 keV)/J(1274 keV). The result of such procedure is shown 
in Fig. 5. As we can see. this simple procedure does not suppress the effect of instability. Because we are 



0.740 
0.730 

1 0.710 

^a.Too 

0.S90 
0.600 
0170 




r 1.0*0 

1.039 
i.0» ? 

* 

■1.000 ^ 

f 

0.M0 Z 



OlSTO 
4.M0 



m ' ' joo' 1 »so" ' 260' ' ' 'iio' JOO 

T**+9*n*w (X) 

fig. 1 S^mvntttn srkhoui ihe decorobdo* procedure tad the rt6o * of tte S-p*w*r it 31 1 k*V to ihe S-ptrmxna m 1774 keV. 



514 



P. Bandfychei aL/Nuct ttffr. orfMttk im Phyt Res. A 3&4( 1997) 506-515 



0.760 
0.7*0 
a 720 
jj 0-700 
1 0.090 



0.6*0 
0.620 



*f S70H.V (-'Bi) 



o o 



A' "ite"ifc"ii6 iio ifc ito i*s iio 



interested only in the relative changes of the ^parameter on temperature, we can make a deconvolution of the 
51! keV peak with the 1274 keV peak and try to suppress the effect of instability. 

Fit 6 shows the result of these decomioluuons (the ^parameter after deconvolution) along with the S- 
parameter fiom deconvolutions of the 511 keV peak with the 570 keV peak ("^i). which was measured 
at the beginning of our experiment. So we can see that the structure at the temperature of 175 K was only 
a consequence of instability of the measuring system, because it disappeared after the deconvolution of the 
511 keV peaks with the 1274 keV peaks. We can obtain physical information from the dependence of the 
5-oaameter on the temperature by using a deconvolution procedure. 

This procedure can also be used for a reliable estimate of the electron momentum distribution, fn such case 
we need a better resolution function (eg. 514 keV from the decay of "Sr) measured simultaneously with the 
annihilation line. 



5. Conclusion 

It turned out that the described deconvolution methods are in principle convenient for processing experimental 
spectra from DBAL measurements. They help to sup press influence of instability of the measuring system. 
Part~ulari, the Gold deconvolution method works very well and it does not deform the deconvolution result 
The reproducibility of the deconvolution results tested on simulated convolutions of Gaussian is very good. 
Van Cittert's algorithm gives the results with oscillations from positive to negative values on the wings of the 
deconvolulcd peaks. Gold's deccavduisea gives only small non-negative oscillations. Unambiguously we can 
say that Gold's algorithm is more reliable with regard to the precision of the deconvolution. The time needed 
for the deconvolution is practically the same for both algorithms. 

So, the use of Gold deconvolution is a cheap and simple way to eliminate the electronic instabilities which 
blur physical effects. 
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